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.     

A multitransition molecular line study of inward motions towards massive star-forming cores

A multitransition 3-mm molecular line single pointing and mapping survey was carried out towards 29 massive star-forming cores in order to search for the signature of inward motions. Up to seven different transitions, optically thick lines HCO⁺(1-0), CS(2-1), HNC(1-0), HCN(1-0) and ¹²CO(1-0), and optically thin lines C¹⁸O(1-0) and ¹³CO(1-0) were observed towards each source. The normalized velocity differences (     ) between the peak velocities of optically thick lines and optically thin line C¹⁸O(1-0) for each source were derived. Prominent inward motions are probably present in either HCO⁺(1-0) or CS(2-1) or HNC(1-0) observations in most sources. Our observations show that there is a significant difference in the incidence of blueshifted asymmetric line profiles between CS(2-1) and HCO⁺(1-0). The HCO⁺(1-0) shows the highest occurrence of obvious asymmetric features, perhaps owing to different optical depth between CS(2-1) and HCO⁺(1-0). HCO⁺(1-0) appears to be the best inward motion tracer. The mapping observations of multiple line transitions enable us to identify six strong infall candidates: G123.07-6.31, W75(OH), S235N, CEP-A, W3(OH) and NGC 7538. The infall signature is extended up to a linear scale  >0.2 pc  .     

Discovery of a dense bipolar outflow from a new class 0 protostar in NGC 2068/LBS 17

We report the discovery of high-velocity dense gas from a bipolar outflow source near NGC 2068 in the L1630 giant molecular cloud. CO and HCO⁺ J =3→2 line wings have a bipolar distribution in the vicinity of LBS 17-H with the flow orientated roughly east–west and perpendicular to the elongation of the submillimetre dust continuum emission. The flow is compact (total extent ∼0.2 pc) and contains of the order of 0.1 M_⊙ of swept-up gas. The high-velocity HCO⁺ emission is distributed over a somewhat smaller area <0.1 pc in extent.
A map of C¹⁸O J =2→1 emission traces the LBS 17 core and follows the ambient HCO⁺ emission reasonably well, with the exception of the direction towards LBS 17-H where there is a significant anticorrelation between the C¹⁸O and HCO⁺. A comparison of beam-matched C¹⁸O and dust-derived H₂ column densities suggests that CO is depleted by up to a factor of ∼50 at this position if the temperature is as low as 9 K, although the difference is substantially reduced if the temperature is as high as 20 K. Chemical models of collapsing clouds can account for this discrepancy in terms of different rates of depletion on to dust grains for CO and HCO⁺.
LBS 17-H has a previously known water maser coincident with it but there are no known near-infrared, IRAS or radio continuum sources associated with this object, leading to the conclusion that it is probably very young. A greybody fit to the continuum data gives a luminosity of only 1.7 L_⊙ and a submillimetre-to-bolometric luminosity ratio of 0.1, comfortably satisfying the criteria for classification as a class 0 protostar candidate.     

HE 1127–1304: the largest radio quasar

By considering the propagation of low-amplitude magnetohydrodynamic waves in partially ionized plasmas, it is shown that the ion-neutral drift (ambipolar diffusion) induced by the waves can have specific effects on the molecular chemistry of cold material. The chemistry occurring in gas swept by Alfvén waves is described and it is shown that this leads to spatial variations in the deuterium fractionation ratios of, for example, HCO⁺ and N₂H⁺, on spatial scales of a few hundredths of a parsec, depending upon the fractional ionization of the ambient medium. The possibility of detecting interstellar Alfvén waves by molecular spectroscopy and their effect of producing small-scale chemical abundance gradients in molecular clouds are briefly discussed.     

A search for CO+ in planetary nebulae

We have carried out a systematic search for the molecular ion CO⁺ in a sample of eight protoplanetary and planetary nebulae in order to determine the origin of the unexpectedly strong HCO⁺ emission previously detected in these sources. An understanding of the HCO⁺ chemistry may provide direct clues for the physical and chemical evolution of planetary nebulae. We find that the integrated intensity of the CO⁺ line may be correlated with that of HCO⁺, suggesting that the reaction of CO⁺ with molecular hydrogen may be an important formation route for HCO⁺ in these planetary nebulae.     

The CO Cameron band emission in the Red Rectangle

The origin of the strong CO Cameron band emission in the Red Rectangle is investigated. From a comparison of laboratory spectroscopic data and astronomical observations in combination with chemical modelling, it is concluded that the Cameron band emission of CO originates from the dissociative recombination of HCO⁺, and its intensity provides a measure of the product of the density of HCO⁺ ions and the electron density, integrated along the line of sight.     

C-type shocks in the interstellar medium: profiles of CH+ and CH absorption lines

We have computed optical absorption-line profiles of CH⁺ and CH, as predicted by a model of a C-type shock propagating in a diffuse interstellar cloud. Both these species are produced in the shock wave in the reaction sequence that is initiated by C⁺(H₂, H)CH⁺. Whilst CH⁺ flows at the ion speed, CH, which forms in the dissociative recombination reaction CH⁺₃(e⁻, H₂)CH, flows at a speed which is intermediate between those of the ions and the neutrals. The predicted velocity shift between the CH⁺ and CH line profiles is found to be no more than approximately 2 km s⁻¹, which is smaller than has previously been assumed. We also investigate OH and HCO⁺, finding that the correlation between their column densities, recently observed in the diffuse interstellar medium, can be reproduced by the model.     

Molecular line mapping of the giant molecular cloud associated with RCW 106 – III. Multimolecular line mapping

We present multimolecular line maps obtained with the Mopra telescope towards the southern giant molecular cloud (GMC) complex G333, associated with the H  ii region RCW 106. We have characterized the GMC by decomposing the 3D data cubes with gaussclumps , and investigated spatial correlations among different molecules with principal component analysis (PCA). We find no correlation between clump size and linewidth, but a strong correlation between emission luminosity and linewidth. PCA classifies molecules into high- and low-density tracers, and reveals that HCO⁺ and N₂H⁺ are anticorrelated.     

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.     

Application of a clumpy core model to NGC 2071

The multi-transitional observations of CS molecules towards the NGC 2071 core have been re-analysed by using a tri-dimensional Monte Carlo radiative transfer code. Better agreement with the observations is made by an introduction of clumpiness to this model than by smoothly varying density to the 1D microturbulent one. The best-fitting model shows that, when a unique density is assumed for clumps, the volume filling factor of the clumps varies as r ⁻² with an average of ∼5 per cent over the entire core, and that the H₂ number density and the CS abundance of the clump relative to H₂ are ∼ 2 × 10⁶ cm⁻³ and ∼ 6 × 10⁻¹⁰, respectively. The radial density gradient ∝ r ⁻² obtained from our clumpy core model is steeper than that (∝ r ^−1.3) obtained from the microturbulent model. Since all clumps are subject to random bulk motions in this 3D clumpy macroturbulent model, synthesized line profiles do not show self-absorption dips even for opaque transitions and the resulting linewidth is in good accordance with the observations.     

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  .     

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

Gas-phase models for the evolved planetary nebulae NGC 6781, M4-9 and NGC 7293

We have studied the chemistry of the molecular gas in evolved planetary nebulae. Three pseudo-time-dependent gas-phase models have been constructed for dense (10⁴–10⁵ cm⁻³) and cool ( T ∼15 K) clumpy envelopes of the evolved nebulae NGC 6781, M4-9 and NGC 7293. The three nebulae are modelled as carbon-rich stars evolved from the asymptotic giant branch to the late planetary nebula phase. The clumpy neutral envelopes are subjected to ultraviolet radiation from the central star and X-rays that enhance the rate of ionization in the clumps. With the ionization rate enhanced by four orders of magnitude over that of the ISM, we find that resultant abundances of the species HCN, HNC, HC₃N and SiC₂ are in good agreement with observations, while those of CN, HCO⁺, CS and SiO are in rough agreement. The results indicate that molecular species such as CH, CH₂, CH₂⁺ , HCl, OH and H₂O are anticipated to be highly abundant in these objects.     

State-to-state rotational transition rates of the HCO+ ion by collisions with helium

The rates of rotational transitions for HCO⁺, the most abundant ion in interstellar space, induced by collision with helium are obtained for temperatures ranging from 10 to 80 K. The calculations are based on a new potential energy surface for the He–HCO⁺ interaction and on a scattering matrix whose accuracy was checked by pressure broadening and shift measurements. The rates     decrease for increasing values of j and  Δ j   , with a temperature trend depending on the energy involved in the transitions: if it is small, the rates are almost constant, while an increase with T is found for other cases. Comparison with previous and less accurate results shows an agreement within 50 per cent. Comparison between state-to-state and pressure broadening cross-sections allows us to discuss importance and influence of elastic and inelastic collisions.     

Electron-impact rotational excitation of linear molecular ions

Molecular R -matrix calculations are performed to give rotational excitation rates for electron collisions with linear molecular ions. Results are presented for CO⁺, HCO⁺, NO⁺ and H₂⁺ up to electron temperatures of 10 000 K. De-excitation rates and critical electron densities are also given. It is shown that the widely used Coulomb–Born approximation is valid for Δ j =1 transitions when the molecular ion has a dipole greater than about 2D, but otherwise is not reliable for studying electron-impact rotational excitation. In particular, transitions with Δ j >1 are found to have appreciable rates and are found to be entirely dominated by short-range effects.     

Cyanopolyynes in hot cores: modelling G305.2+0.2

We present results from a time-dependent gas-phase chemical model of a hot core based on the physical conditions of G305.2+0.2. While the cyanopolyyne HC₃N has been observed in hot cores, the longer chained species, HC₅N, HC₇N and HC₉N, have not been considered as the typical hot-core species. We present results which show that these species can be formed under hot core conditions. We discuss the important chemical reactions in this process and, in particular, show that their abundances are linked to the parent species acetylene which is evaporated from icy grain mantles. The cyanopolyynes show promise as 'chemical clocks' which may aid future observations in determining the age of hot core sources. The abundance of the larger cyanopolyynes increases and decreases over relatively short time-scales,  ∼10^2.5 yr  . We present results from a non-local thermodynamic equilibrium statistical equilibrium excitation model as a series of density, temperature and column density dependent contour plots which show both the line intensities and several line ratios. These aid in the interpretation of spectral-line data, even when there is limited line information available. In particular, non-detections of HC₅N and HC₇N in Walsh et al. are analysed and discussed.     

Deep echelle spectrophotometry of S 311, a Galactic H ii region located outside the solar circle

We present echelle spectrophotometry of the Galactic H  ii region S 311. The data have been taken with the Very Large Telescope Ultraviolet-Visual Echelle Spectrograph in the 3100–10 400 Å range. We have measured the intensities of 263 emission lines; 178 are permitted lines of H⁰, D⁰ (deuterium), He⁰, C⁰, C⁺, N⁰, N⁺, O⁰, O⁺, S⁺, Si⁰, Si⁺, Ar⁰ and Fe⁰; some of them are produced by recombination and others mainly by fluorescence. Physical conditions have been derived using different continuum- and line-intensity ratios. We have derived He⁺, C⁺⁺ and O⁺⁺ ionic abundances from pure recombination lines as well as abundances from collisionally excited lines for a large number of ions of different elements. We have obtained consistent estimations of t ² applying different methods. We have found that the temperature fluctuations paradigm is consistent with the T _e(He  i ) versus T _e(H  i ) relation for H  ii regions, in contrast with what has been found for planetary nebulae. We report the detection of deuterium Balmer lines up to Dδ in the blue wings of the hydrogen lines, whose excitation mechanism seems to be continuum fluorescence.     

Turbulence in Class 0 and I protostellar envelopes

We estimate the levels of turbulence in the envelopes of class 0 and I protostars using a model based on measurements of the peak separation of double-peaked asymmetric line profiles. We use observations of 20 protostars of both class 0 and I taken in the  HCO⁺(J = 3 → 2)  line that show the classic double-peaked profile. We find that some class 0 sources show high levels of turbulence, whilst others demonstrate much lower levels. In class I protostars, we find predominantly low levels of turbulence. The observations are consistent with a scenario in which class 0 protostars form in a variety of environments and subsequently evolve into class I protostars. The data do not appear to be consistent with a recently proposed scenario in which class 0 protostars can only form in extreme environments.     

Australia Telescope Compact Array 1.2-cm observations of the massive star-forming region G305.2+0.2

We report on Australia Telescope Compact Array observations of the massive star-forming region G305.2+0.2 at 1.2 cm. We detected emission in five molecules towards G305A, confirming its hot core nature. We determined a rotational temperature of 26 K for methanol. A non-local thermodynamic equilibrium excitation calculation suggests a kinematic temperature of the order of 200 K. A time-dependent chemical model is also used to model the gas-phase chemistry of the hot core associated with G305A. A comparison with the observations suggest an age of between  2 × 10⁴  and  1.5 × 10⁵ yr  . We also report on a feature to the south-east of G305A which may show weak Class I methanol maser emission in the line at 24.933 GHz. The more evolved source G305B does not show emission in any of the line tracers, but strong Class I methanol maser emission at 24.933 GHz is found 3 arcsec to the east. Radio continuum emission at 18.496 GHz is detected towards two H  ii regions. The implications of the non-detection of radio continuum emission towards G305A and G305B are also discussed.     

Chemical composition of the Orion nebula derived from echelle spectrophotometry

We present echelle spectroscopy in the 3500- to 7060-... range for two positions of the Orion nebula. The data were obtained using the 2.1-m telescope at Observatorio Astronómico Nacional in San Pedro Mártir, Baja California. We have measured the intensities of about 220 emission lines, in particular 81 permitted lines of C⁺, N⁺, N⁺⁺, O⁰, O⁺, Ne⁰, Si⁺, Si⁺⁺ and S⁺, some of them produced by recombination only and others mainly by fluorescence. We have determined electron temperatures, electron densities and ionic abundances using different continuum and line intensity ratios. We derived the He, C and O abundances from recombination lines and find that the C/H and O/H values are very similar to those derived from B stars of the Orion association, and that these nebular values are independent of the temperature structure. We have also derived abundances from collisionally excited lines. These abundances depend on the temperature structure; accurate t ² values have been derived comparing the O II recombination lines with the [O III ] collisionally excited lines. The gaseous abundances of Mg, Si and Fe show significant depletions, implying that a substantial fraction of these atoms is tied up in dust grains. The derived depletions are similar to those found in warm clouds of the Galactic disc, but are not as large as those found in cold clouds. A comparison of the solar and Orion chemical abundances is made.