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.     

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.     

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.     

The high-velocity molecular stream in W51B: a ring structure with compact H ii regions

We present ¹³ CO J  = 1 − 0 line observations of the H  ii region complex W51B located in the high-velocity (HV) stream. These observations reveal a filamentary and clumpy structure in the molecular gas. The mean local standard of rest (LSR) velocity ∼ + 65 km s⁻¹ of the molecular gas in this region is greater than the maximum velocities allowed by kinematic Galactic rotation curves. The size and mass of the molecular clouds are ∼ 48 × 17 pc² and ∼ 2.4 × 10⁵ M⊙ respectively. In a position–velocity diagram, molecular gas in the southern part comprises a redshifted ring structure with v _LSR=+ 60 to +73 km s⁻¹. The velocity gradient of this ring is ∼ 0.5 km s⁻¹ pc⁻¹, and the mass is ∼ 6.2 × 10⁴ M⊙. If we assume that the ring is expanding with a uniform velocity, the expansion velocity, radius and kinetic energy are ∼ 7 km s⁻¹, ∼ 13 pc and ∼ 3.0 × 10 ⁴⁹ erg respectively. The kinetic energy and mass spectrum of the ring could be explained by an expanding cylindrical cloud with a centrally condensed mass distribution. The locations of two compact H  ii regions, G49.0−0.3 and G48.9−0.3, coincide with the two molecular clumps in this ring. We discuss star formation, and the mechanism that produced the ring structure.     

The interstellar medium in the Upper Cepheus-Cassiopeia region

We have studied the kinematics and spatial distribution of the interstellar gas in the sky region  110°≤ l ≤ 135°, 10°≤ b ≤ 20°  , using the extensive Leiden–Dwingeloo Survey of H  i emission and the Columbia Survey of CO emission. The spectra show two main velocity components, namely feature A that has a mean local standard of rest (LSR) velocity of  ∼0  km s⁻¹  and is due to the Lindblad ring of the Gould belt, and feature C that has a mean LSR velocity of  ∼−11  km s⁻¹  and is associated to the local arm or Orion arm. The H  i and CO distributions of feature A in the region trace a large complex of gas and dust known as the Cepheus Flare, which lies at a distance of 300 pc. The spectral line profiles of feature A, which are rather broad and often double-peaked, reveal that the Cepheus Flare forms part of a big expanding shell of interstellar matter that encloses an old supernova remnant associated with a void inside the Cepheus Flare. On the other hand, by analysing the distribution and velocity structure of feature C, we have detected a second large expanding shell in the region, located at a distance of 800 pc in the local arm. This shell surrounds the stellar association Cepheus OB4 and was probably generated by stellar winds and supernovae of Cepheus OB4. The radii, expansion velocities and H  i masses of the two shells are approximately 50 pc, 4  km s⁻¹ and  1.3 × 10⁴ M_⊙  for the Cepheus Flare shell and 100 pc, 4 km s⁻¹ and  9.9 × 10⁴ M_⊙  for the Cepheus OB4 shell. Both shells have similar ages of the order of a few 10⁶ yr.     

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.     

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.     

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

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.     

A Herbig–Haro object associated with GGD30 and its exciting source

GGD30 has been suggested to be either a small reflection nebulosity or a Herbig–Haro (HH) object formed in the outflow from a nearby obscured star. Observations to date have not been able to distinguish between these two scenarios. In addition, there are conflicting proposals for the location of the exciting source for GGD30. To resolve these questions, we have carried out optical spectroscopy and near-infrared ( J , K and 3.6-μm) imaging of GGD30. Taken together, these observations reveal that the bright optical knot in GGD30 must be a HH object, excited by the outflow from an optically obscured pre-main-sequence (PMS) star located ∼3 arcsec to the southwest. Based on mid-infrared fluxes from the Mid-course Space Experiment ( MSX ) satellite, we estimate the luminosity of this PMS star to be  ∼12.5 L_⊙  which suggests it is an intermediate-mass object rather than low-mass as previously proposed. The optical spectroscopy indicates projected velocities of  ∼−270 km s⁻¹  associated with the HH object. The fact that these velocities are blueshifted and relatively high compared to the velocities typical of HH flows suggests that the outflow from the PMS star must be almost aligned with the line of sight. There is an additional low-velocity  (∼−70 km s⁻¹) Hα  component but its origin is not clear.     

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.     

Chimneys in the starburst galaxy M82

We present MERLIN and VLA observations at 1.4 and 5 GHz of the diffuse radio emission in the centre of M82. We detect a large expanding shell of ionized gas surrounding the brightest supernova remnant 41.95+57.5 with a diameter of ∼100 pc and an expansion velocity of ∼100 km s⁻¹. We observe a 50-pc-scale 'blow-out' from this region, in the form of a 'cone' of missing 5-GHz continuum emission, which appears to be an excellent example of a galactic chimney.
On larger radio scales, we observe four chimney structures extending out to the north ∼100–200 pc along the minor axis. One of these features is almost certainly related to the 50-pc-scale blow-out from 41.95+57.5, although this is not the most prominent feature. The other features have also been traced to expulsion of material from the very centre by using an 'unsharp masked' image from 5-GHz VLA B-array observations, with the supernova remnant removed.
We interpret these chimney features as the base of the superwind, which implies that the ejection of material into the halo does not occur smoothly over the starburst region. Instead, very localized venting of hot gas is clearly in evidence.     

Molecular hydrogen emission in Cygnus A

We present J , H and K -band spectroscopy of Cygnus A, spanning 1.0–2.4 μm in the rest-frame and hence several rovibrational H₂, H recombination and [Fe  ii ] emission lines. The lines are spatially extended by up to 6 kpc from the nucleus, but their distinct kinematics indicate that the three groups (H, H₂ and [Fe  ii ]) are not wholly produced in the same gas. The broadest line, [Fe  ii ] λ 1.644, exhibits a non-Gaussian profile with a broad base (FWHM≃1040 km s⁻¹), perhaps because of the interaction with the radio source. Extinctions to the line-emitting regions substantially exceed earlier measurements based on optical H recombination lines.
Hard X-rays from the quasar nucleus are likely to dominate the excitation of the H₂ emission. The results of Maloney, Hollenbach & Tielens are thus used to infer the total mass of gas in H₂ v=1–0 S(1)-emitting clouds as a function of radius, for gas densities of 10³ and 10⁵ cm⁻³, and stopping column densities N _H=10²²–10²⁴ cm⁻². Assuming azimuthal symmetry, at least 2.3×10⁸ M_⊙ of such material is present within 5 kpc of the nucleus, if the line-emitting clouds see an unobscured quasar spectrum. Alternatively, if the bulk of the X-ray absorption to the nucleus inferred by Ueno et al. actually arises in a circumnuclear torus, the implied gas mass rises to ∼10¹⁰ M_⊙. The latter plausibly accounts for 10⁹ yr of mass deposition from the cluster cooling flow, for which within this radius.     

Modelling the bar in the centre of the starburst galaxy M82

We present VLA A-array 21-cm atomic hydrogen (H  i ) absorption observed against the central region of the starburst galaxy M82 with an angular resolution of ∼1.3 arcsec (≃20 pc). These observations, together with MERLIN H  i absorption measurements, are compared with the molecular (CO) and ionized ([Ne  ii ]) gas distributions and are used to constrain the dynamics and structure of the ionized, neutral and molecular gas in this starburst.
A position–velocity diagram of the H  i distribution reveals an unusual 'hole' feature which, when previously observed in CO, has been interpreted as an expanding superbubble contained within a ring of gas in solid body rotation. However, we interpret this feature as a signature of a nearly edge-on barred galaxy. In addition, we note that the CO, H  i and [Ne  ii ] position–velocity diagrams reveal two main velocity gradients, and we interpret these as gas moving on x₁- and x₂-orbits within a bar potential. We find the best fit to the data to be produced using a bar potential with a flat rotation curve velocity v _b=140 km s⁻¹ and a total length of 1 kpc, a non-axisymmetry parameter q =0.9, an angular velocity of the bar Ω_b=217 km s⁻¹ arcsec⁻¹, a core radius R _c=25 pc, an inclination angle i =80° and a projected angle between the bar and the major axis of the galaxy φ '=4°. We also discuss the orientation of the disc and bar in M82.     

Ultra-high-resolution observations of CH in Southern Molecular Cloud envelopes

We present a mini-survey of ultrahigh-resolution spectroscopy (UHRS) of CH towards three southern molecular cloud envelopes. The sightlines are selected to probe physically similar gas in different Galactic environments. With a velocity resolution of ∼0.5 km s⁻¹  ( R =575 000)  these observations resolve most kinematic components of the absorption lines. We do, however, detect one line component in the Lupus region, which is not resolved and for which an upper limit of   b <0.3 km s^-1  is found. We find a correlation between distance of the absorbing gas from the Galactic mid-plane and the fractional abundance of CH. We show that this correlation can be explained as being a result of a fall-off in the ultraviolet radiation field intensity and propose that CH observations in carefully selected sightlines might allow a mapping of the variations in the interstellar radiation field.     

The coronal line regions of planetary nebulae NGC 6302 and 6537: 3–13 μm grating and echelle spectroscopy

We report on advances in the study of the cores of NGC 6302 and 6537 using infrared grating and echelle spectroscopy. In NGC 6302, emission lines from species spanning a large range of ionization potential, and in particular [Si  ix ] 3.934 μm, are interpreted using photoionization models (including cloudy ), which allow us to re-estimate the temperature of the central star to be about 250 000 K. All of the detected lines are consistent with this value, except for [Al  v ] and [Al  vi ]. Aluminium is found to be depleted to one hundredth of the solar abundance, which provides further evidence for some dust being mixed with the highly ionized gas (with photons harder than 154 eV). A similar depletion pattern is observed in NGC 6537. Echelle spectroscopy of IR coronal ions in NGC 6302 reveals a stratified structure in ionization potential, which confirms photoionization to be the dominant ionization mechanism. The lines are narrow (<22 km s⁻¹ FWHM), with no evidence of the broad wings found in optical lines from species with similar ionization potentials, such as [Ne  v ] 3426 Å. We note the absence of a hot bubble, or a wind-blown bipolar cavity filled with a hot plasma, at least on 1 arcsec and 10 km s⁻¹ scales. The systemic heliocentric velocities for NGC 6302 and 6537, measured from the echelle spectra of IR recombination lines, are found to be −34.8±1 km s⁻¹ and −17.8±3 km s⁻¹. We also provide accurate new wavelengths for several of the infrared coronal lines observed with the echelle.     

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.     

Sputtering of grains in C-type shocks

Sputtering yields are reported for the release of Mg, Fe, Si and O under impact of He, C, O, Si and Fe on grain material composed of Mg- and Fe-bearing silicates. The yields were derived using the trim code, which simulates the results of the transport of ions in matter by means of classical Monte Carlo techniques. The energetics of the sputtering process are a key factor in the sputtering calculations, and so detailed determinations have been made of the energy with which atoms are bound to the lattice, using solid-state simulation programs. The probability of ejection of an atom is computed at a given energy, for a number of angles of incidence, and integrated to obtain the mean yield at that energy. These numerical results are then fitted with a simple function of energy for convenience in subsequent applications.
A grid of C-type shock models has been computed, using our new sputtering yields, for pre-shock densities in the range 10⁴ n _H n (H)+2 n (H₂)10⁶ cm⁻³ and shock speeds 20 v _s45 km s⁻¹. Sputtered fractions can be high, exceeding 50 per cent for shock speeds in excess of approximately 40 km s⁻¹. The column densities of Si and SiO were also computed, for comparison with observations.     

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.     

Kinematics of massive stars in the Small Magellanic Cloud

We present radial velocities for 2045 stars in the Small Magellanic Cloud (SMC), obtained from the 2dF survey by Evans et al. The great majority of these stars are of OBA type, tracing the dynamics of the young stellar population. Dividing the sample into ad hoc 'bar' and 'wing' samples (north and south, respectively, of the line:  δ=−77°50'+[4α]'  , where α is the right ascension in minutes of time) we find that the velocities in the SMC bar show a gradient of 26.3 ± 1.6 km s⁻¹ deg⁻¹ at a position angle of 126°± 4°. The derived gradient in the bar is robust to the adopted line of demarcation between the two samples. The largest redshifts are found in the SMC wing, in which the velocity distribution appears distinct from that in the bar, most probably a consequence of the interaction between the Magellanic Clouds that is predicted to have occurred 0.2 Gyr ago. The mean velocity for all stars in the sample is +172.0 ± 0.2 km s⁻¹ (redshifted by ∼20 km s⁻¹ when compared to published results for older populations), with a velocity dispersion of 30 km s⁻¹.