Highlight
Free Access
Issue
A&A
Volume 560, December 2013
Article Number A95
Number of page(s) 10
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201322164
Published online 11 December 2013

© ESO, 2013

1. Introduction

Although interstellar clouds have ionization fractions of only 10-4–10-8, ionic species are very useful to probe physical conditions in such clouds (Larsson et al. 2012). In diffuse clouds (AV < 1), the main ion source is UV photoionization of carbon, while in dense clouds (AV > 1), cosmic-ray ionization of hydrogen is the dominant ionization mechanism (Bergin & Tafalla 2007). Proton transfer reactions of interstellar H with abundant species such as CO and N2 lead to HCO+ and N2H+, which are widely observed in the interstellar medium. Such stable ionic species are useful as tracers of the interaction of interstellar gas with magnetic fields (Houde et al. 2004; Schmid-Burgk et al. 2004), whereas ions which react rapidly with H2 trace other parameters such as the gas density and the ionization rate.

At temperatures ≲ 250 K, the formation of interstellar H2O in the gas phase proceeds through a series of ion-molecule reactions. After charge transfer of H+ or H to O, repeated reactions of O+ with H2 produce OH+, H2O+, and finally H3O+, which upon dissociative recombination with a free electron produces H2O. The OH+ and H2O+ ions are well known from the spectra of comets where they appear as photodissociation products of H2O (Swings & Page 1950; Herzberg & Lew 1974). In the interstellar medium, H2O and H3O+ have been known for decades (Phillips et al. 1992), but observation of the intermediate products OH+ and H2O+ had to await the launch of ESA’s Herschel Space Observatory (Pilbratt et al. 2010). Strong absorption in rotational lines of interstellar OH+ and H2O+ is seen with Herschel on many lines of sight in our Galaxy (Gerin et al. 2010; Ossenkopf et al. 2010) and even some beyond (Weiß et al. 2010). In addition, electronic absorption lines of OH+ have been reported in sensitive near-UV spectra of several diffuse interstellar clouds (Krełowski et al. 2010).

The Herschel data, as well as the single OH+ line observed from the ground toward the Galactic Center source SgrB2 (Wyrowski et al. 2010), imply large column densities of OH+ and H2O+. The hydrogen in the absorbing clouds thus cannot be purely in atomic form, because no OH+ and H2O+ would be produced, nor in purely molecular form, because all OH+ and H2O+ would react into H3O+ and H2O. Using models of UV-irradiated interstellar clouds (PDRs), the observed abundances of OH+ and H2O+ can be used to infer the relative fractions of hydrogen in atomic and molecular forms (Neufeld et al. 2010), which itself traces the ionization rates of the clouds (Hollenbach et al. 2012).

While the interpretation of interstellar HnO+ absorption is reasonably well understood, lines of OH+ and H2O+ have also been observed in emission from the nuclei of several active galaxies, most famously Mrk 231 (Van der Werf et al. 2010). The large dipole moments and small reduced masses of the HnO+ ions imply high line frequencies and large radiative decay rates, so that collisional excitation of their rotational levels requires extremely high densities and line emission is not expected to be observable. Understanding this phenomenon benefits from finding a Galactic source of HnO+ line emission, which can be studied in more detail than extragalactic nuclei.

This paper presents the first observation of OH+ line emission toward a source within our Galaxy: the Orion Bar. Due to its brightness and nearly edge-on geometry, this PDR is well-suited to observe physical and chemical changes in the gas as a function of depth into the cloud, as the intensity of UV irradiation by the Trapezium stars decreases (e.g., Van der Wiel et al. 2009). The Orion Bar is also notable as the only known Galactic source of interstellar HF line emission (Van der Tak et al. 2012). In this case, the proximity of this region (420 pc: Menten et al. 2007; Hirota et al. 2007) allows us to resolve the HnO+ line emission both spatially and spectrally. We use non-local thermodynamic equilibrium (non-LTE) radiative transfer models and PDR thermo-chemical models to interpret our results.

2. Observations

The CO+ peak (αJ2000 = 05h35m20.6s, δJ2000 =  −05°25′14′′) in the Orion Bar (Störzer et al. 1995) has been observed as a spectral scan over the full HIFI range as part of the HEXOS (Herschel observations of EXtra-Ordinary Sources) guaranteed-time key program (Bergin et al. 2010) using the HIFI instrument (De Graauw et al. 2010) of the Herschel Space Observatory (Pilbratt et al. 2010). This paper uses data from HIFI bands 1b (H2O+ 110–101), 3b (OH+ 10–01), 4a (OH+ 12–01 and 11–01), 5a (H2O+ 111–000), and 6b (H3O+). These observations were carried out in 2011 March and April in load chop mode with a redundancy of 4, except that frequency switching was used in band 5. On-source integration times are ≈ 50 s for most spectra except OH+ 10–01 (20 s) and H2O+ 111–000 (190 s). The ObsIDs of the spectra, without the leading 1342, are 215 923 for Band 1b, 216 380 for Band 3b, 218 628 for Band 4a, 194 666 for Band 5a and 218 426 for Band 6b.

Table 1 lists the frequencies of the lines as well as other spectroscopic parameters, which have been taken from the CDMS database (Müller et al. 2001)1. The size of the telescope beam for these observations is 19–23′′ FWHM, corresponding to 9000 AU or 0.04 pc, except for the H2O+ 110–101 line where it is 35′′, and for the H3O+ line where it is 15′′. The WBS (Wide-Band Spectrometer) was used as backend, covering 4 GHz bandwidth in four 1140 MHz subbands at 1.1 MHz resolution. The velocity calibration of HIFI data is accurate to  ~ 0.5 km s-1 or better. The data were reduced with HIPE (Ott 2010) pipeline version 6.0, using the task doDeconvolution for the sideband deconvolution, while further analysis was done in the CLASS package.

Table 1

Observed lines.

In addition to the HIFI spectral scans, a 115′′ × 65′′ area centered on α = 05:35:20.81, δ = –05:25:17.1 with a position angle of 145° was mapped in the OH+ 12 − 01 transition with HIFI, in on-the-fly (OTF) mapping mode with position-switch reference, using a total integration time of 20 min. The noise level of the map is 0.22 K per 0.7 km s-1 channel and its ObsID is 218 216. The fully sampled map was reduced with HIPE 6.0 and exported to CLASS for further analysis. We compare these data to a map of the same area in the CO 10–9 line at ν = 1151.985 GHz (Eup = 304 K), also observed within the HEXOS program under ObsID 217736, with a noise level of 0.56 K per 1.0 km s-1 channel.

3. Results

thumbnail Fig. 1

Line profiles of the OH+ 12–01 (top) and 11–01 (bottom) transitions, observed with HIFI toward the CO+ peak of the Orion Bar, after smoothing to 0.6 km s-1 resolution.

Open with DEXTER

Table 2

Line parameters measured from the HIFI spectra.

3.1. Line profiles

Figure 1 presents the line profiles of the OH+ lines at 971 and 1033 GHz, as observed with HIFI toward the Orion Bar. Although hints of OH+ and H2O+ emission have been seen toward the high-mass protostar W3 IRS5 (Benz et al. 2010), this is the first time that a pure emission profile is seen toward a Galactic source. As the line profiles appear single-peaked at our sensitivity and spectral resolution, we have fitted a Gaussian model to extract the line parameters. Table 2 presents the results of these Gaussian fits, and also reports our upper limits to other lines of OH+, H2O+ and H3O+ from our spectral survey of the Orion Bar. We have also searched for the O+ excited fine structure (2D J = 5/2 → 3/2) line near 593.7 GHz (Blagrave & Martin 2004; Sharpee et al. 2004), leading to an upper limit on  of 30 mK rms per 0.5 km s-1 channel. No continuum emission is detected in the spectra, down to upper limits ranging from = 0.1 K at 607 GHz to 0.3 K at 1655 GHz, which is consistent with the SPIRE measurements of Arab et al. (2012).

The difference between the measured central velocities and FWHM widths of the two detected OH+ lines is mostly due to hyperfine blending of the 12–01 line. We thus adopt the measured parameters of the 11–01 line as the best estimate of the central velocity and FWHM width of the OH+ emitting gas in the Orion Bar. While the central velocity of this line of 10.1 km s-1 is in good agreement with the value of 10.0 ± 0.2 km s-1 from ground-based observations of other molecular species toward this source, its width of 4.3 km s-1 is much larger than the value of 1.7 ± 0.3 km s-1 measured for the dense gas deep inside (AV  ≳  1) the Bar (Hogerheijde et al. 1995; Leurini et al. 2006). On the other hand, its width is smaller than the value of ≈ 5 km s-1 measured in species such as HF which trace the surface (AV  ≲  0.1) of the PDR (Van der Tak et al. 2012), and more similar to the width of the [C II] line (Ossenkopf et al. 2013). Based on the line width, the OH+ emission thus seems to originate from close to the PDR surface (0.3 < AV < 0.5) where CH+ and SH+ peak as well in the chemical model of Nagy et al. (2013). In contrast, the non-hydride reactive molecular ions CO+, SO+ and HOC+ have smaller line widths (2–3 km s-1, Fuente et al. 2003). These species have a different formation channel and thus a different sensitivity to the abundances of atomic H and electrons.

3.2. Spatial distribution

Figure 2 presents our maps of the OH+ 971 GHz line emission, as observed with HIFI toward the Orion Bar, integrated over three velocity ranges, shown in greyscale. The contours show the CO 10–9 line emission, integrated over the same velocity ranges. The OH+ emission around the central velocity (middle panel) is seen to be extended over at least an arc minute (25 000 AU or 0.12 pc) on the sky, and to roughly follow the structure of the Bar seen in CO 10–9 and other molecular tracers (e.g., Van der Wiel et al. 2009). The emission from the Bar is concentrated in two clumps: one peaking near the map center at VLSR  ≈  10 km s-1, and another peaking 20–25′′ to the east near VLSR = 12 km s-1, which is seen most pronounced in the figure’s bottom panel. This second peak lies ≈ 10′′ south-east of the CO 10–9 peak at the same velocity, which coincides with H13CN clumps 2 and 3 from Lis & Schilke (2003). Since OH+ is not expected deeper into the PDR than CO and HCN, we suggest that the location of the second OH+ clump marks a deviation from pure plane-parallel geometry.

The map of low-velocity OH+ in the top panel of Fig. 2 shows a third clump which lies to the north-west of the Bar, and is connected to it by a bridge of fainter OH+ emission. This “perpendicular” emission is also seen in CO 10–9 and other tracers of the Bar surface such as OH 119 μm (Goicoechea et al. 2011), [CII] 158 μm (Ossenkopf et al. 2013), and CH+ 3–2 (Nagy et al. 2013). The feature corresponds to the southern tip of the Orion Ridge facing the Trapezium cluster, which confines the HII region, as visible in large-scale maps of CN line emission (Rodríguez-Franco et al. 1998) and 13CO 3–2 emission (Lis & Schilke 2003; Melnick et al. 2012). From multi-line CN observations, Rodríguez-Franco et al. (2001) derive H2 densities as high as 106 cm-3 for this ridge, which may be favourable to excite OH+ line emission. However, this density is probably an overestimate, as collisions of CN with electrons were not taken into account, which are known to be important for CN (Black & van Dishoeck 1991; Harrison et al. 2013). Furthermore, we note that the OH+ line is not detected toward the Orion S clump further up the Orion Ridge, possibly because emission and absorption from different layers cancel out each other.

thumbnail Fig. 2

Spatial distribution of the OH+ 12–01 emission, integrated between VLSR = +6 and +10 km s-1 (top), between +8 and +12 km s-1 (middle), and between +10 and +14 km s-1 (bottom). Greyscale levels start at 0.6 K km s-1 and increase by 0.2 K km s-1. Contours of CO 10–9 emission are at 5, 15, ... 95% of the peak intensity in the respective velocity channel. Vibrationally excited H2 peaks near the north-west edge of the CO emission (cf. Fig. 1 of Nagy et al. 2013).

Open with DEXTER

3.3. Column densities

The variations in OH+ emission level by factors of 2–3 across the maps in Fig. 2 likely correspond to variations in the total column density or the OH+ abundance with position, although excitation conditions (volume density, kinetic temperature) may also play a role. We focus in the following analysis on the CO+ peak where we have limits on the other HnO+ ions. To estimate the column densities of these ions, we use the line fluxes from Table 2 and apply a main beam efficiency of 76% at 607 GHz, 70% at 1655 GHz, and 74% at our other line frequencies (Roelfsema et al. 2012). The column density N depends on the excitation temperature Tex through where ν is the line frequency, Q is the partition function, gu is the upper state degeneracy, and Aul is the spontaneous decay rate. This equation assumes optically thin emission and negligible background radiation (Tbg ≪ Tex), but does not assume the Rayleigh-Jeans limit ( ≪ kTex) because our measurements are at high frequency. The optically thin assumption is reasonable given the low expected abundances of the HnO+ ions, and neglecting background radiation is justified given the low level of continuum radiation in our observations. Adopting a background brightness temperature of 9 K increases the column density estimates by 5% for Tex = 20 K, and by < 1% for higher values of Tex.

To evaluate the above expression for the column density, we use the spectroscopic parameters of the lines in Table 1. Table 3 presents estimates of the column densities of OH+, H2O+ and H3O+ for values of Tex between 10 and 160 K, which is the expected range for the Orion Bar. If the excitation of HnO+ is close to LTE, Tex would be close to the kinetic temperature of the gas, which ranges from ≈ 85 K for the dense gas (Hogerheijde et al. 1995) to ≈ 150 K near the cloud surface (Batrla & Wilson 2003; Goicoechea et al. 2011). If collisional excitation of the lines cannot compete with their radiative decay, Tex will drop below Tkin, while Tex may exceed Tkin if radiative or chemical pumping plays a role. Section 5 discusses these processes in more detail. For now, the table shows that the derived column densities of HnO+ vary by factors of 2–3 for Tex ≳ 20 K, but that the estimate increases by an order of magnitude if Tex is as low as 10 K.

The OH+ column density in Table 3 for Tex = 10 K is comparable to the values toward other Galactic sources, while the estimates for Tex ≳  20 K are ~ 10 ×  lower. Absorption line studies indicate N(OH+) values between a few 1013 and a few 1014 cm-2, both for the diffuse foreground clouds toward G10.6, W49N, W51 and OMC-2 (Gerin et al. 2010; Neufeld et al. 2010; Indriolo et al. 2012; López-Sepulcre et al. 2013) and the dense gas around the protostars AFGL 2591, W3 IRS5, and Orion-KL (Bruderer et al. 2010; Benz et al. 2010; Gupta et al. 2010). Our OH+/H2O+ ratio of > 40 and OH+/H3O+ ratio of > 15 are larger than in previous observations of diffuse clouds, suggesting an origin of the observed emission in very diffuse low-extinction layers of the PDR.

Table 3

Optically thin estimates of HnO+ column densities toward the CO+ peak of the Orion Bar, in units of 1012 cm-2, as a function of assumed excitation temperature.

thumbnail Fig. 3

Abundances of OH+, H2O+, H3O+, and C+, predicted by the Meudon PDR model as a function of visual extinction for a radiation field of χ = 10   000 in units of the local interstellar radiation field according to Draine (1978) and a pressure of P = 108 K cm-3, as applicable to the Orion Bar.

Open with DEXTER

4. PDR models

To understand the chemistry of HnO+ in more detail, we model the Orion Bar with version 1.4.4 of the Meudon PDR code (Le Petit et al. 2006; Goicoechea & Le Bourlot 2007; González-García et al. 2008). We describe the object with an isobaric model, where the pressure is kept constant and the program solves for the thermal and chemical balance as a function of depth. Isobaric models are the simplest description of gas in steady-state, stationary molecular clouds where gravitation is negligible. The model includes a gas-phase chemical network based on the UMIST database, as well as H2 formation on grains and neutralization of ions on grains and PAHs.

Figure 3 shows the predicted abundances of OH+, H2O+ and H3O+ as a function of visual extinction, for a radiation field of χ = 10   000 Draine (1978) units (χ0 = 2.7 × 10-3 erg s-1 cm-2), a cosmic-ray ionization rate of 2 × 10-16 s-1,and a pressure of P = 108 K cm-3. The radiation field is taken from Goicoechea et al. (2011) and the pressure is chosen to match the intensities of the CH+J = 1–0 to 6–5 transitions observed toward the Orion Bar CO+ peak with Herschel HIFI and PACS (Nagy et al. 2013). For these parameters, the total gas density at the depth where the HnO+ abundances peak (AV ~ 0.3 − 1.0) is in the range between 5.6 × 104 and 1.8 × 105 cm-3. The temperatures at these depths in the model (Table 5) agree well with estimates for the PDR surface from lines of H2(400–700 K, Allers et al. 2005) and H (~550 K, Van der Werf et al. 2013). The column densities of our observed species, integrated up to AV = 1 mag, are N(OH+) = 1.6 × 1013 cm-2, N(H2O+) = 4.4 × 1012 cm-2, and N(H3O+) = 7.5 × 1012 cm-2. When integrating up to AV = 10 mag, the values for OH+ and H2O+ remain the same to 1%, and N(H3O+) increases by just 6% which is insignificant. These predictions are consistent with the estimates in Table 3 for Tex = 10–20 K, except that the observed OH+/H2O+ ratio is larger than in the models. Note that the predictions correspond to a face-on model; to compare to the observations, they must be multiplied by a factor of 1/sin i ≈ 4, where i ≈ 15° for the Orion Bar (Nagy et al. 2013).

Table 4

Rate coefficients of the main formation and destruction paths of HnO+.

Table 5

Rates (in cm-3 s-1) corresponding to the main formation and destruction paths of HnO+ in the model at AV = 0.3, 0.4, and 1.0.

As shown in Fig. 3, HnO+ abundances in the high UV illumination Orion Bar PDR peak near the cloud surface, at depths of AV ~ 0.3−0.4, and decrease rapidly beyond AV ~ 1. Given the clumpy structure of the Orion Bar, this narrow range in AV is spatially rather extended, as indicated by the maps in Fig. 2. To understand the chemistry of HnO+ in the Orion Bar, we study the main formation and destruction channels of these ions at depths of AV = 0.3, 0.4 and 1.0 mag. The corresponding reaction coefficients and rates are listed in Tables 4 and 5. At a depth of AV = 0.3–0.4, the dominant pathways for OH+ formation are: The dominant OH+ formation path is via H2 and O+, which accounts for 70% of the total OH+ production at AV = 0.3 and 85% at AV = 0.4. Photoionization of OH accounts for 18% of OH+ formed at AV = 0.3 and 12% at AV = 0.4. Charge exchange between OH and H+ is mostly significant at AV ~ 0.3, resulting in 12% of OH+ production which decreases to 3% at AV ~ 0.4. This reaction is more important in X-ray dominated regions (XDRs), which is a key difference between such regions and high-illumination PDRs (Van der Werf et al. 2010). At a depth of AV ~ 1, the path via H2 and O+ is still the dominant OH+ formation path (72.1%). However, the role of OH photoionization increases, as it accounts for producing 27% of OH+ at this depth.

Once OH+ is formed, H2O+ and H3O+ formation rapidly follows via similar reactions with H2: At a depth of AV = 0.3, 100% of H2O+ forms from OH+ in the reaction with H2. At a depth of AV = 0.4, as the H3O+ abundance is increased compared to that at AV = 0.3, a secondary reaction apart from the reaction from H2 (83%) produces 17% of the H2O+ formed at this depth: (6)At a depth of AV = 1.0, the role of H2O+ production from H3O+ decreases to 9%, compared to the dominant reaction from H2 (88%). At each of these depths in the model,100% of H3O+ is produced in the reaction of H2O+ with H2.

The destruction of OH+ is mainly via reactions with electrons, photons and H2. The reaction with H2 (Eq. (4)) is the most important H2O+ formation path, as mentioned above. It is also the most important destruction path for OH+ at depths of AV = 0.4 and beyond. At smaller depths, both dissociative recombination with electrons and UV photodissociation contribute significantly to the OH+ destruction: at AV ~    0.3, (7)contributes 52%, and (8)contributes 18%. Destruction of H2O+ occurs mainly via dissociative recombination with electrons at low depths (AV = 0.3–0.4): Reaction 9 accounts for 71% of H2O+ destruction at AV = 0.3, while reaction 10 contributes 20%. At a depth of AV = 0.4, reactions 9 and 10 contribute 44 and 12% (respectively) to the H2O+ destruction. As H2O+ does not have photodissociation channels longward of 13.6 eV, the model does not include this process (Hollenbach et al. 2012). Destruction by H2, which controls H3O+ production is only significant at AV ~ 1 (87%). Therefore, the high electron density and UV radiation field in the Orion Bar may explain the H2O+ and H3O+ non-detections, as a significant fraction of OH+ and H2O+ is directly destroyed by UV photons or by recombination with electrons, limiting the H2O+ and H3O+ formation. Our non-detections of H2O+ and H3O+ are thus consistent with an origin of the OH+ emission in diffuse gas where C+ is abundant.

Alternatively, our observed H2O+/OH+ ratio may be due to a low molecular fraction (H2/H ratio) in the gas probed by our observations. Gerin et al. (2010) derive an analytic expression for the OH+/H2O+ abundance ratio as a function of the gas temperature, electron density, and H2 volume density: As most OH+ forms in the outermost layers of the PDR (AV < 0.4) at very low molecular fractions (0.01%–0.4% for AV = 0.3 and AV = 0.4, respectively), this formula is consistent with the observed abundance ratio of N(OH+)/N(H2O+)  >  2. This suggests that, apart from the effect of the high electron density mentioned above, the low molecular fraction in the surface of the PDR also contributes to the observed OH+/H2O+ column density ratios. We conclude that most of our observed OH+ emission originates at extinctions below AV = 0.4.

Figure 3 shows a second peak in the H3O+ abundance at a depth of AV ~ 7 − 8. This second peak is expected in the abundance of HnO+ in interstellar clouds with a large range of physical conditions (e.g., Hollenbach et al. 2012) and does not significantly contribute to the total H3O+ column density. The sequence of HnO+ formation at this depth is initiated by the cosmic-ray ionization of H2, followed by a reaction of H with H2 producing H. At this depth, OH+ is formed from H, whose abundance is increased due to the lower electron abundance. The second abundance peak disappears toward higher AV as oxygen freezes out on the grain surfaces as water.

thumbnail Fig. 4

Face-on column densities of OH+, H2O+ and H3O+, predicted by the Meudon PDR model as a function of radiation field (in Draine units) for a pressure of P = 108 K cm-3, as applicable to the Orion Bar. The integrated AV of the model is 10 mag.

Open with DEXTER

Figure 4 shows the effect of changing the radiation field on the HnO+ column densities. The model corresponds to P = 108 cm-3 K and the figure considers radiation fields between χ = 5000 and 50 000 Draine units. The predicted column densities of the HnO+ species are seen to peak for χ = 20 000–30 000 and to drop for higher radiation fields. The HnO+ column densities in our model are similar to those by Hollenbach et al. (2012) under Orion Bar conditions (their Figs. 8 and 9).

5. Excitation of OH+

Our estimated column densities of OH+, H2O+ and H3O+ are in reasonable agreement with the predictions from the Meudon PDR model, but the appearance of the lines in emission is unlike other Galactic sources observed in OH+ so far. To understand this behaviour we perform a non-LTE analysis of the excitation of OH+, which takes both reactive and inelastic collisions into account, as well as excitation by the background continuum radiation field.

5.1. Collisional and radiative excitation

The calculations in Appendix A indicate inelastic electron collision rates for the OH+ 971 GHz line of  ≈ 6 × 10-7cm3 s-1, which is ~ 10 ×  higher than the dissociative recombination rate at T ≈ 100 K (Table 4). At this temperature, /kT ≈ 0.5 for this transition, so that by detailed balance, the upward and downward rates differ by only a few percent. For collisions with H2, inelastic collision rates are not known, but it may be reasonable to assume that reactive collisions dominate. In contrast, the reaction of OH+ with H is endothermic, and inelastic collisions with H could influence the excitation of OH+. Test calculations assuming a collisional rate coefficient of 10-10 cm3 s-1 (cf. Andersson et al. 2008 for the case of CO+) for all radiatively allowed transitions indicate that this effect is comparable with electron collisions. Very recent calculations on the OH+-He system confirm this estimate, at least to order of magnitude (F. Lique, priv. comm.). We therefore first calculate the excitation of OH+ assuming that inelastic collisions with H and electrons dominate, and consider the effect of reactive collisions with H2 in Sect. 5.2.

We have used the non-LTE radiative transfer program Radex (Van der Tak et al. 2007) to calculate the excitation of OH+ in the Orion Bar, assuming steady-state conditions and using the inelastic collision rates from Appendix A. Observations of C recombination lines toward the Orion Bar surface indicate hydrogen densities of 5 × 104–2.5 × 105 cm-3(Wyrowski et al. 1997), which for all carbon in C+ and C/H = 1.4 × 10-4 translates into an electron density of n(e) ≈ 10 cm-3, as we have used before for HF (Van der Tak et al. 2012). We use an electron temperature of Te = 300 K, but our results are insensitive to variations in Te between 100 and 1000 K. The adopted line width is the observed 4.3 km s-1 and for the background radiation field we adopt a modified blackbody distribution with a dust temperature of Td = 50 K and a dust emissivity index of β = 1.6, as found by Arab et al. (2012) for the interior of the Bar, so that τd = 0.21 at 971 GHz. The model is insensitive to the details of this radiation field; in particular, the results are unchanged when adopting Td = 70 K and β = 1.2 as found by Arab et al. for the Bar’s surface.

The model predicts an excitation temperature of ≈10 K for the OH+ lines near 1 THz. The excitation is due to the combination of collisions and radiation: a model with the modified blackbody replaced by the 2.73 K cosmic microwave background results in Tex  ≈  7 K. The effect on the emerging line intensities is large, since the Planck function at 1 THz increases by a factor of  ~4 from Tex = 7 K to 9 K.

The calculated excitation temperature of ≈10 K is above the background radiation temperature of 9.2 K, which explains the appearance of the line in emission, but the models require N(OH+)  ≈  5 × 1014 cm-2 to match our observed 971 GHz line intensity. At this high column density, the OH+ lines become optically thick which drives the three N = 1–0 lines to approximately the same peak brightness, inconsistent with our non-detection of the 10 − 01 line of OH+. We therefore regard this model as untenable, conclude that the observed line emission is optically thin, and search for a model with a higher excitation temperature in the next section.

5.2. Effect of reactive collisions

The above calculation indicates that inelastic collisions are sufficient to make the OH+ line appear in emission given the low background intensity of the Orion Bar. A fully self-consistent model should however include formation and destruction terms in the rate equations. We use Eq. (12) from Van der Tak et al. (2007), where the number density of OH+ molecules (in cm-3) is the ratio of its formation rate (in cm-3 s-1) and its destruction rate (in s-1). We adopt a destruction rate of OH+ in the Orion Bar of 10-4 s-1, which is due to reactions with H2, electron recombination, and photodissociation. In contrast, the formation of OH+ through ion-molecule reactions should proceed at a rate at or below the Langevin rate of ≈ 10-9 cm3 s-1 and is limited by the supply of hydrogen or oxygen ions.

Lacking state-to-state formation rates, we approximate the distribution of the newly-formed OH+ over its energy levels by a thermal distribution at a temperature Tf. The value of Tf depends on the dominant formation route and is likely to be a significant fraction of its excess energy. The reactions of O+ with H2 and of O with H have exothermicities of 0.55 and 0.66 eV (Federer et al. 1984; Milligan & McEwan 2000), so Tf should be in the range 2000–3000 K.

For OH+ formation rates between 10-12 and 10-9 cm3 s-1 and formation temperatures of 2000–3000 K, the model predicts excitation temperatures of  ≈ 12 K for the OH+ lines near 1 THz. Compared to models with a negligible formation rate, the line brightness increases by a factor of ≈ 5, almost independent of Tf, for the same column density. The observed intensities of the 12 − 01 and 11 − 01 lines are matched for N(OH+)  ≈  1 × 1014 cm-2, which is considerably lower than for the steady-state excitation model. The difference with the estimates in Table 3 is the presence of background radiation. This model predicts some emission in the 10 − 01 line, but only at the 2σ level of our observations. We conclude that formation pumping plays an important role for OH+ in the Orion Bar.

6. Discussion

6.1. Effect of X-ray ionization

Our derived OH+ column density of ≈1 × 1014 cm-2 is similar to that in diffuse clouds, as observed in absorption toward W49N and other sources (see references in Sect. 3.3), but  ~1.6 ×  higher than the model prediction of 6.4 × 1013 cm-2 in Sect. 4. The calculations in Sect. 4 assume a cosmic-ray ionization rate of 2 × 10-16 s-1, which is a typical value for diffuse interstellar clouds in the Solar neighbourhood (Indriolo & McCall 2012) and an order of magnitude higher than the value for dense clouds (Van der Tak & van Dishoeck 2000). However, the actual ionization rate of the Orion Bar may be atypically high because of its proximity to the Trapezium stars. The effects of X-rays from these stars on the chemistry are similar to those of cosmic rays (Meijerink et al. 2006). In particular, Gupta et al. (2010) estimate an ionization rate of 3 × 10-15 s-1 for the Orion KL region, dominated by X-rays from the star θ1C Ori. Calculations using the Meudon code with an ionization rate of 2 × 10-15 s-1 result in very similar HnO+ column densities as before, though, presumably because of the high gas density in the Orion Bar.

6.2. Possible ion sources

Alternatively, the discrepancy between our observed and modeled column densities of OH+ may arise because the assumption of a stationary medium is inappropriate. Regarding the Orion Bar as a molecular cloud under external illumination and heating, the main source of OH+ may be the reaction of H with O, where the H is due to cosmic-ray ionization as well as leakage of UV photons from the ionized region. However, interaction of the PDR with the neighbouring photoionized nebula may supply O+ ions which produce extra OH+ in their reaction with H2. The supply of O+ ions may be due to leakage of hydrogen- and oxygen-ionizing photons, but also to advection of H+ and O+ ions from the ionized nebula into the mostly neutral PDR, as in some planetary nebulae (Black 1983)2. The concentration of H+ with depth into the PDR is then important, because charge transfer reactions (H+ + O ↔  O+ + H) will rapidly couple the O+/O ratio to the H+/H ratio. Observations of the [Oiii] 88 μm line toward the Orion Bar with Herschel/PACS show that ionized oxygen is widespread in the region (C. Joblin & J. Goicoechea, priv. comm.); the emission extends well into the mostly-neutral gas traced by the [Oi] 63 μm line (M. Gerin, priv. comm.), making the O+ + H2 channel a likely source of OH+ at low depths. Detailed comparison of the optical and near-IR (forbidden and/or permitted) lines of O and O+ would be a stronger test of this scenario (Walmsley et al. 2000; Mesa-Delgado et al. 2011). In addition, simulations of the full molecular + atomic + ionized gas in the Orion PDR with the Cloudy program (Ferland et al. 2013) would be useful to constrain the role of possible ion sources in the Orion Bar.

6.3. Radiative pumping

Besides changing the chemistry of OH+, the strong infrared and ultraviolet radiation fields in the Orion Bar may change its excitation. Section 5.1 already showed that far-infrared continuum radiation raises the OH+ excitation temperature significantly above the level due to inelastic collisions alone. In addition, mid-infrared pumping through the v = 1–0 band at 3.38 μm may contribute if the radiative excitation rate of OH+ in the Bar exceeds the collisional excitation rate by electrons, which is n(e)  ×  Clu = 10 cm-3 ×  6 × 10-7 cm3 s-1 = 6 × 10-6 s-1. The radiative rate is BluUrad, the Einstein absorption coefficient times the radiative energy density, which is approximately Avibϵf/(e(/(kTd)) − 1), where ϵ is the dust emissivity, f the dust filling factor, Avib the spontaneous decay rate of 265 s-1, and Td the dust temperature. For the Orion Bar, the filling factor should be close to unity, and for the emissivity we assume 1 at the short wavelength of the fundamental vibrational band of OH+. Equating the radiative rate to the electron collision rate indicates a minimum temperature of ≈ 240 K, which is reasonable for the gas at the Bar’s surface, but too much for the dust, as the PACS and SPIRE data show (Arab et al. 2012). Only PAHs and small grains would reach such high temperatures, but with very low opacities, even though the OH+ vibrational fundamental is close to the PAH 3.3 μm and aliphatic 3.4 μm emission features.

Besides continuum radiation from dust in the Bar itself, pumping by infrared starlight from the Trapezium may influence the excitation of OH+. The brightest of these stars is θ1C Ori, which is 127′′ away from our observing position (Van Leeuwen 2007). If both objects lie at the same distance from the Sun, the stellar continuum flux at the Orion Bar is 2.96 million times stronger than at the Earth. The stellar temperature of 37 000 K implies a radiative intensity of 0.2 Jy/nsr or 2 × 10-15 erg s-1 cm-2 Hz-1 sr-1. Setting this equal to the Planck function at a radiation temperature TR, we obtain TR = 185 K at λ = 3.38   μm. The corresponding pumping rate in the OH+ vibrational fundamental is of order Avib/ { exp(/kTR) − 1 }  = 5.7 × 10-8 s-1, which is much less than the collisional excitation rate by electrons.

While infrared pumping does not seem play a role for OH+ in the Orion Bar, the absorption rate through electronic transitions in the near ultraviolet is rather higher. The radiative intensity of θ1C Ori at the wavelength of the Av = 0–0 band of OH+ at 27 949 cm-1 (3577 Å) is ≈ 1.8 Jy/nsr. The A-value for this band is 8.01 × 105 s-1(de Almeida & Singh 1981), so that the absorption rate in this band alone is ≈ 3.3 × 10-6 s-1, which is only slightly less than the collisional excitation rate.

6.4. Comparison with extragalactic systems

The Orion Bar is the first and so far only position within our Galaxy where lines of OH+ appear purely in emission, although a mix of emission and absorption is seen in W3 IRS5 (Benz et al. 2010) and NGC 3603 (Makai et al., in prep.) and possibly many more star-forming regions. Detections of extragalactic OH+ and H2O+ emission have been made with Herschel-SPIRE by Van der Werf et al. (2010) toward the active nucleus of the galaxy Mrk 231, with Herschel-PACS toward the ultraluminous merger Arp 220 by Rangwala et al. (2011) and with Herschel-SPIRE toward the Seyfert nucleus NGC 1068 by Spinoglio et al. (2012). In addition, detections of extragalactic HnO+ absorption exist toward M 82 using HIFI (Weiß et al. 2010) and SPIRE (Kamenetzky et al. 2012). Recently, PACS observations of excited OH+, H2O+ and H3O+ toward NGC 4418 and Arp. 220 have been discussed by González-Alfonso et al. (2013). The Orion Bar is special in our Galaxy for its large column density of warm, mostly-atomic gas and its weak far-infrared continuum; we suspect that extragalactic nuclei where HnO+ lines appear in emission have similar conditions. In addition, supernova remnants may contribute, as the recent detection of OH+ emission towards the Crab nebula suggests (Barlow et al. 2013).

7. Conclusions

We have presented maps and spectra of OH+ line emission toward the Orion Bar, and limits on lines of H2O+ and H3O+. The OH+ line emission is extended over ~ 1′ (=25 000 AU = 0.12 pc) and traces the Bar itself as well as the southern tip of the Orion Ridge. Analysis of the chemistry and the excitation of OH+ suggests an origin of the emission at a depth of AV = 0.3–0.5, similar to CH+ and SH+. The OH+ column density of ≈ 1.0 × 1014 cm-2, derived using a non-LTE model including both inelastic and reactive collisions and radiative pumping, is similar to that in previous absorption line studies, while our limits on the OH+/H2O+ and OH+/H3O+ ratios are higher than seen before.

Non-LTE models of the excitation of OH+ show that the unusual appearance of the OH+ lines in emission is the combined result of inelastic collisions, far-infrared radiative pumping by dust, and chemical pumping through the O+ + H2 and O + H channels. The same conditions may apply to extragalactic sources of HnO+ line emission. In the future, high-resolution maps of Galactic and extragalactic OH+ line emission with ALMA will shed further light on the chemistry of this reactive ionic species.

Our observed HnO+ column densities are qualitatively reproduced by a model of the Orion Bar (using the Meudon PDR code) using a radiation field of χ = 104 χ0 and a pressure of P = 108 K cm-3 as suggested by previous observations. Analysis of the main formation and destruction paths of the ions indicates that our high OH+/H2O+ and OH+/H3O+ ratios are due to the high UV radiation field and electron density in the Orion Bar. Destruction of OH+ and H2O+ by photodissociation and electron recombination limits the formation of H2O+ and H3O+. In addition, the low molecular fraction at the PDR surface limits the production of H2O+ and H3O+.

Quantitatively, the Meudon PDR models underpredict the absolute OH+ column density by a factor of ~ 1.6. To match the observed line intensity with an OH+ column density similar to that in the PDR model, the electron density would have to be ≈ 100 cm-3, which is much higher than the PDR model predicts at the depth where the HnO+ ions are abundant, as seen in Fig. 10 of Nagy et al. (2013). However, raising the pressure in the PDR model by ≈ 50% would increase the predicted N(OH+) to the value suggested by the non-LTE models. Such an increase is consistent with the CH+ and SH+ observations, and is also suggested by observations of high-J CO lines with PACS (Joblin et al., in prep.). Furthermore, the Meudon PDR model uses a scaling of the average interstellar radiation field, while realistic models should use direct observations of the dominant hot star θ1C Ori for the H-ionizing part of the spectrum and the wavelength range where OH photoionization occurs, which could contribute significantly to the formation of OH+.


2

Very recently, OH+ emission has been detected towards the Helix nebula (Etxaluze et al. Etxaluze et al. (2013)).

Acknowledgments

The authors thank Sébastien Bardeau and Jérôme Pety (IRAM) for help with the data reduction, Franck le Petit (Paris-Meudon) for assistance with the PDR model calculations, Arturo Rodríguez-Franco (Madrid) for sending his CN map in electronic form, Xander Tielens (Leiden) for useful discussions, Inga Kamp (Groningen) for a careful reading of the manuscript, Christine Joblin (Toulouse) & Javier Goecoechea (Madrid) for sharing their PACS observations of the [Oiii] line, and the referee for a useful report. AF acknowledges support by the Agence Nationale de la Recherche (ANR-HYDRIDES), contract ANR-12-BS05-0011-01. HIFI has been designed and built by a consortium of institutes and university departments from across Europe, Canada and the US under the leadership of SRON Netherlands Institute for Space Research, Groningen, The Netherlands with major contributions from Germany, France and the US. Consortium members are: Canada: CSA, U.Waterloo; France: CESR, LAB, LERMA, IRAM; Germany: KOSMA, MPIfR, MPS; Ireland, NUI Maynooth; Italy: ASI, IFSI-INAF, Arcetri-INAF; Netherlands: SRON, TUD; Poland: CAMK, CBK; Spain: Observatorio Astronómico Nacional (IGN), Centro de Astrobiología (CSIC-INTA); Sweden: Chalmers University of Technology – MC2, RSS & GARD, Onsala Space Observatory, Swedish National Space Board, Stockholm University – Stockholm Observatory; Switzerland: ETH Zürich, FHNW; USA: Caltech, JPL, NHSC.

References

Appendix A: Inelastic collision rates for the OH+–e system

The electronic ground state symmetry of the radical OH+ is . Each rotational level N is therefore split by the spin-rotation coupling between N and the electronic spin S = 1 so that each rotational level N has two sub-levels given by j = N ± 1. In addition, owing to the non-zero nuclear spin of the hydrogen atom (I = 1/2), each fine-structure level is further split into 2 hyperfine levels F = j ± 1/2. The rotational constant of OH+ is 492.26 GHz. The fine-structure splitting is of the order of 60–90 GHz while the hyperfine splitting is less than 0.3 GHz. The dipole moment of OH+ is 2.26 D (Werner et al. 1983).

Electron-impact hyperfine excitation rate coefficients for OH+ were computed using a three-step procedure: i) rotational excitation rate coefficients for the dipolar (ΔN = 1) transitions were first computed within the Coulomb-Born approximation; ii) fine-structure excitation rate coefficients were then obtained from the Coulomb-Born rotational rates using the (scaled) infinite-order-sudden (IOS) approximation; and iii) hyperfine excitation rate coefficients were finally obtained using the so-called “statistical” or “proportional” approach. The Coulomb-Born approximation (Chu & Dalgarno 1974) is expected to be accurate for polar molecules with dipoles in excess of  ~2 D because the dipolar cross sections are entirely dominated by long-range effects and cross sections for transitions with ΔN ≥ 2 are significantly smaller (Faure & Tennyson 2001). In practice, Coulomb-Born cross sections were computed for collision energies below 2 eV and rate coefficients were deduced for temperatures ranging from 10 to 2000 K. The IOS approximation was employed to derive the fine-structure rate coefficients in terms of the rotational rates for excitation out of the lowest rotational level N = 0. This IOS formalism was first introduced by Corey & McCourt (1984) for linear molecules with 2S + 1Σ symmetry. As the Coulomb-Born rotational rates do not strictly obey the IOS factorization formulae, however, the IOS fine-structure rate coefficients were scaled, as recommended by Faure & Lique (2012; see their Eqs. (8), (10) and (13), where the quantum number F should be replaced by S). While in principle the IOS approximation should be also applicable to obtain the hyperfine rate coefficients of OH+, there is to our knowledge no available factorization formula for a 3Σ molecule. Hyperfine rate coefficients were therefore obtained from the fine-structure rates by assuming that they are proportional to the degeneracy (2F + 1) of the final hyperfine level. We note that this simple statistical approach does not account for the collisional propensity rule ΔF = Δj. As shown by Faure & Lique (2012), however, at low total optical depth (τ ≤ 10) the statistical approach is applicable because in this regime the relative populations of each hyperfine component are close to the statistical weights.

The above three-step procedure was applied to the first 49 levels of OH+, that is up to the level (N,j,F) = (8,8,17/2) which lies 1689 K above the ground state (0, 1, 3/2), resulting in 176 collisional transitions. A typical accuracy of 30% is expected for these rate coefficients, with the largest rate coefficients being of the order of 3 × 10-6 cm3 s-1. The rates will be posted on the website of the LAMDA database (Schöier et al. 2005)3.

All Tables

Table 1

Observed lines.

Table 2

Line parameters measured from the HIFI spectra.

Table 3

Optically thin estimates of HnO+ column densities toward the CO+ peak of the Orion Bar, in units of 1012 cm-2, as a function of assumed excitation temperature.

Table 4

Rate coefficients of the main formation and destruction paths of HnO+.

Table 5

Rates (in cm-3 s-1) corresponding to the main formation and destruction paths of HnO+ in the model at AV = 0.3, 0.4, and 1.0.

All Figures

thumbnail Fig. 1

Line profiles of the OH+ 12–01 (top) and 11–01 (bottom) transitions, observed with HIFI toward the CO+ peak of the Orion Bar, after smoothing to 0.6 km s-1 resolution.

Open with DEXTER
In the text
thumbnail Fig. 2

Spatial distribution of the OH+ 12–01 emission, integrated between VLSR = +6 and +10 km s-1 (top), between +8 and +12 km s-1 (middle), and between +10 and +14 km s-1 (bottom). Greyscale levels start at 0.6 K km s-1 and increase by 0.2 K km s-1. Contours of CO 10–9 emission are at 5, 15, ... 95% of the peak intensity in the respective velocity channel. Vibrationally excited H2 peaks near the north-west edge of the CO emission (cf. Fig. 1 of Nagy et al. 2013).

Open with DEXTER
In the text
thumbnail Fig. 3

Abundances of OH+, H2O+, H3O+, and C+, predicted by the Meudon PDR model as a function of visual extinction for a radiation field of χ = 10   000 in units of the local interstellar radiation field according to Draine (1978) and a pressure of P = 108 K cm-3, as applicable to the Orion Bar.

Open with DEXTER
In the text
thumbnail Fig. 4

Face-on column densities of OH+, H2O+ and H3O+, predicted by the Meudon PDR model as a function of radiation field (in Draine units) for a pressure of P = 108 K cm-3, as applicable to the Orion Bar. The integrated AV of the model is 10 mag.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.