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 shallow though extensive H2 2.122-μm imaging survey of Taurus–Auriga–Perseus – I. NGC 1333, L1455, L1448 and B1

We discuss wide-field near-infrared (near-IR) imaging of the NGC 1333, L1448, L1455 and B1 star-forming regions in Perseus. The observations have been extracted from a much larger narrow-band imaging survey of the Taurus–Auriga–Perseus complex. These H₂ 2.122-μm observations are complemented by broad-band K imaging, mid-IR imaging and photometry from the Spitzer Space Telescope , and published submillimetre CO   J = 3–2  maps of high-velocity molecular outflows. We detect and label 85 H₂ features and associate these with 26 molecular outflows. Three are parsec-scale flows, with a mean flow lobe length exceeding 11.5 arcmin. 37 (44 per cent) of the detected H₂ features are associated with a known Herbig–Haro object, while 72 (46 per cent) of catalogued HH objects are detected in H₂ emission. Embedded Spitzer sources are identified for all but two of the 26 molecular outflows. These candidate outflow sources all have high near-to-mid-IR spectral indices (mean value of  α∼ 1.4  ) as well as red IRAC 3.6–4.5 μm and IRAC/MIPS 4.5–24.0 μm colours: 80 per cent have [3.6]–[4.5] > 1.0 and [4.5]–[24] > 1.5. These criteria – high α and red [4.5]–[24] and [3.6]–[4.5] colours – are powerful discriminants when searching for molecular outflow sources. However, we find no correlation between α and flow length or opening angle, and the outflows appear randomly orientated in each region. The more massive clouds are associated with a greater number of outflows, which suggests that the star formation efficiency is roughly the same in each region.     

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.     

ISO observations of molecular hydrogen in the DR 21 bipolar outflow

We demonstrate that a wide range of molecular hydrogen excitation can be observed in protostellar outflows at wavelengths in excess of 5 μm. Cold H₂ in DR 21 is detected through the pure rotational transitions in the ground vibrational level (0–0). Hot H₂ is detected in pure rotational transitions within higher vibrational levels (1–1, 1–2, etc.). Although this emission is relatively weak, we have detected two 1–1 lines in the DR 21 outflow with the ISO SWS instrument. We thus investigate molecular excitation over energy levels corresponding to the temperature range 1015–15 722 K, without the uncertainty introduced by differential extinction when employing near-infrared data.
This gas is thermally excited. We uncover a rather low H₂ excitation in the DR 21 West Peak. The line emission cannot be produced from single C-shocks or J-shocks; a range of shock strengths is required. This suggests that bow shocks and/or bow-generated supersonic turbulence is responsible. We are able to distinguish this shock-excited gas from the fluoresced gas detected in the K band, providing support for the dual-excitation model of Fernandes, Brand & Burton.     

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.     

The emission and kinematic structures of the bipolar planetary nebula M 1-8

We have undertaken echelle spectroscopy and narrow-band line imaging of the bipolar planetary nebula M 1-8. This has permitted us to map the outflow in [N  ii ]λλ 6548+6583 Å, Hα, and in the v = 1–0 S(1) transition of H₂ at λ 2.122 μm. It has also permitted us to acquire high-resolution spectra for [N  ii ]λ 6583 Å, Hα and He  ii λ 6560 Å. Our observations support the results of a previous 2MASS analysis by two of the authors (J. P. Phillips and G. Ramos-Larios), and confirm that there is strong H₂ emission outside of the ionized zone, as well as along the major axis of the outflow. Finally, we have investigated the spatial structure of the outflow in low and high excitation lines, and noted evidence for strong ionization stratification within the envelope of the source. We also note that major axis spectra show asymmetries attributable to outflow along the lobes, oriented at an angle i ∼ 35°–40° to the line of sight. Asymmetries along the minor axis, by contrast, appear to be associated with the central collimating disc, and may be interpretable in terms of asymmetries in disc structure, or rotation at an angular velocity of Ω∼ 1.4 10⁻¹² rad s⁻¹. If the disc arises due to common-envelope evolution, then it seems that angular momentum constraints must be relatively tight, and can only be satisfied given fairly extreme physical assumptions (such as low disc mass, high primary star mass, a low distance to the source and so forth).     

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.     

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.     

Near-infrared spectroscopy of (proto)-planetary nebulae: molecular hydrogen excitation as an evolutionary tracer

We present an in-depth analysis of molecular excitation in 11 H₂-bright planetary and protoplanetary nebulae (PN and PPN). From newly acquired K -band observations, we extract a number of spectra at positions across each source. H₂ line intensities are plotted on 'column density ratio' diagrams so that we may examine the excitation in and across each region. To achieve this, we combine the shock models of Smith, Khanzadyan & Davis with the photodissociation region (PDR) models of Black & van Dishoeck to yield a shock-plus-fluorescence fit to each data set.
Although the combined shock + fluorescence model is needed to explain the low- and high-energy H₂ lines in most of the sources observed (fluorescence accounts for much of the emission from the higher-energy H₂ lines), the relative importance of shocks over fluorescence does seem to change with evolutionary status. We find that shock excitation may well be the dominant excitation mechanism in the least evolved PPN (CRL 2688 – in both the bipolar lobes and in the equatorial plane) and in the most evolved PN considered (NGC 7048). Fluorescence, on the other hand, becomes more important at intermediate evolutionary stages (i.e. in 'young' PN), particularly in the inner core regions and along the inner edges of the expanding post-asymptotic giant branch (AGB) envelope. Since H₂ line emission seems to be produced in almost all stages of post-AGB evolution, H₂ excitation may prove to be a useful probe of the evolutionary status of PPN and PN alike. Moreover, shocks may play an important role in the molecular gas excitation in (P)PN, in addition to the low- and/or high-density fluorescence usually attributed to the excitation in these sources.     

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.     

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.     

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.     

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.