Free Access
Issue
A&A
Volume 569, September 2014
Article Number L5
Number of page(s) 5
Section Letters
DOI https://doi.org/10.1051/0004-6361/201424756
Published online 19 September 2014

© ESO, 2014

1. Introduction

Hydrides, which consist of one heavy atom and one or more H atoms, are of fundamental importance in the interstellar medium (ISM). Unfortunately, their fundamental transitions at submillimeter (sub-mm) wavelengths are often difficult to observe from the ground. The Herschel Space Observatory and the Stratospheric Observatory for Infrared Astronomy (SOFIA) have increased our knowledge about hydrides in space considerably. Recent detections include ArH+ with Herschel (Barlow et al. 2013; Schilke et al. 2014) and SH with SOFIA (Neufeld et al. 2012).

The reactive molecular ion SH+ (sulfanylium) is an interesting probe of energetic processes in the ISM. In particular, SH+ is only detectable in significant amounts if the very high endothermicity (0.86 eV or ~9860 K) of the gas-phase reaction (1)can be overcome. Using the Atacama Pathfinder EXperiment (APEX) 12 m telescope, Menten et al. (2011) discovered recently sub-mm SH+ absorption lines in the envelope of the prolific star-forming region Sagittarius B2(M), Sgr B2(M) for short, a very strong continuum source close to the Galactic center, and in the diffuse interstellar clouds on the line of sight toward this source. The ubiquity of SH+ in diffuse clouds was re-emphasized soon thereafter by Herschel/HIFI observations toward several Galactic sight lines by Godard et al. (2012). These authors proposed that the required energy for reaction (1) arises from turbulent dissipation, shocks, or shears in these very low density clouds. Sub-mm SH+ emission lines were also detected by Herschel in higher density environments such as the high-mass star-forming region W3 IRS5 (Benz et al. 2010) or the Orion Bar photo-dissociation region (PDR; Nagy et al. 2013). These are strongly UV-irradiated environments where the gas attains high temperatures (≲ 1000 K) and where H2 molecules are UV-pumped to vibrationally excited states. In addition, in these PDR environments, reaction (1) becomes exothermic when H2 molecules are in the ν= 2 or higher vibrational levels, thus enhancing the SH+ abundance (e.g., Zanchet et al. 2013). Finally, if doubly ionized sulfur atoms co-exist with H2 molecules in X-ray dissociation regions (XDRs), such as Galactic center clouds (e.g., Godard et al. 2012; Etxaluze et al. 2013), the very exothermic reaction (2)can be a significant source of SH+ (e.g., Abel et al. 2008).

The strongest ground state FS component of SH+(X) at ~526 GHz (NJ = 12−01) cannot be observed from the ground because of the water in Earth’s atmosphere. The ~683 GHz (NJ = 11−01) and ~893 GHz (NJ = 21−11) transitions can be observed from the ground, albeit with some difficulty. Transitions of the NJ = 10−01 FS component at ~346 GHz (Savage et al. 2004) lie in a more accesible frequency window. Indeed, Menten et al. (2011) searched for the SH+~346, ~683, and ~893 GHz absorption lines toward Sgr B2(M). Unfortunately, the ~346 GHz lines are heavily blended with absorption components of the nearby CO J = 3−2 line caused by Milky Way spiral arm molecular clouds. Stäuber et al. (2007) detected an emission feature at 345 930 MHz in the massive protostar AFGL 2591, but it may be caused by 34SO2, a common species in these sources, so they could not draw firm conclusions about SH+. Leurini et al. (2006) presented first the results of a 1 mm line survey toward the Orion Bar with the APEX 12 m telescope. They did not report the detection of SH+. However, they observed the denser HCN clump (Lis & Schilke 2003) and not the most exposed UV-irradiated PDR gas layers where SH+ is expected.

In this letter we present the first ground-based detection of SH+ lines toward the Orion Bar CO+ peak (Stoerzer et al. 1995; Goicoechea et al. 2011, and references therein). We show that the previously reported SH+NJ = 10−01 line frequencies at ~346 GHz are not correct and use the frequencies derived from our observations to compute improved SH+ spectroscopic parameters and line frequencies.

thumbnail Fig. 1

From top to bottom: IRAM 30 m, ACA, and ALMA spectra toward the Orion Bar between 345.8 GHz and 345.9 GHz. Blue arrows show the predicted frequencies of the SH+, NJ = 10−01, F = 1/2−1/2 and F = 1/2−3/2 lines before this work (based on experiments by Savage et al. 2004). Red arrows (prediction*) show the expected line frequencies for a Hamiltonian fit that excludes the 345 929.8 MHz line observed by Savage et al. (2004) and attributed to SH+; see Sect. 4.

Open with DEXTER

2. Observations

2.1. Single-dish IRAM 30 m observations

IRAM 30 m telescope observations were taken in 2013 using the EMIR330 receiver and the FFTS backend in the wide mode that covers a 16 GHz instantaneous bandwidth per polarization at a channel spacing of 195 kHz. They are part of a complete millimeter line survey toward the Orion Bar dissociation front (Cuadrado et al. 2014) and include specific deep integration searches for SH+NJ = 10−01 lines. The target position is at α2000 = 05h35m20.8s, δ2000 = −05°25′17.0′′, close to the so-called CO+ peak (Stoerzer et al. 1995). The observing procedure was position switching (PSW) with the reference position located at a (600′′, 0′′) offset to avoid the extended molecular emission from the Orion Molecular Cloud. The antenna temperature, , was converted to the main beam temperature, , using an antenna efficiency of 42 %. At this frequency, the IRAM 30 m telescope provides an angular resolution of 7′′. Figure 1 (top panel) shows the resulting spectrum between 345.8 GHz and 345.9 GHz after baseline subtraction. The integration time was 2 h, and the rms noise was ~30 mK per channel. The angular resolution and sensitivity of these observations are a factor of ~2.5 better than those of Leurini et al. (2006). Atmospheric calibration of the IRAM 30 m data was carried out with the ATM program (Cernicharo 1985, IRAM internal report; Pardo et al. 2001).

2.2. Interferometric ALMA and ACA observations

Atacama Large Millimeter Array (ALMA) observations of the Orion Bar were carried out in 2013 and 2014 during ALMA Cycle 1 observations. They belong to the 2012.1.00352.S project “The fundamental structure of molecular cloud edges: from clumps to photoevaporation” (P.I.: Goicoechea). They consist of mosaicking observations covering a ~40′′ × 40′′ field centered at α2000 = 05h35m20.6s, δ2000 = −05°25′20.0′′.

The ALMA array observations used twenty-seven 12 m antennas. The field of view was observed as a mosaic of 27 Nyquist-sampled pointings in the C32-2 configuration, providing a typical resolution of ~1′′. Observations were carried out with the band 7 receiver near 346 GHz using the ALMA correlator backend at a spectral resolution of 488.3 kHz (~0.4 km s-1) over a total bandwidth of 937.5 MHz. Atacama Compact Array (ACA) is used to complement ALMA observations of extended sources by producing short-spacing visibilities filtered out by the ALMA array. At the time of observations, the ACA array used nine 7 m antennas. The ACA mosaic contained ten pointings Nyquist-sampled fields and provided a typical resolution of ~3.6′′. The total observing time was about 6.5 h and 2.2 h for ACA and ALMA observations respectively. The data were calibrated within the CASA software package using the standard algorithms. The calibrated uv tables were then exported to GILDAS1, where the data were deconvolved using the Högbom CLEAN algorithm. All technical details on the data reduction and the complete data set will be presented elsewhere. For the purpose of this work (showing the 345.8345.9 GHz window observed by different instrumentation) we only deconvolved the ACA and ALMA visibilities separately.

3. Observational results

Figure 1 shows, from top to bottom, the IRAM 30 m, ACA, and ALMA spectra toward the Orion Bar dissociation front; we note that ACA and ALMA spectra are averaged over the same region of the dissociation front. The abscissa axis represents the rest frequency (in MHz) using a local standard of rest velocity (vLSR) of 10.50 km s-1. We derived this vLSR value from the accurately known HCO+J = 4−3 line frequency2 observed in the same region of the dissociation front. It agrees, within ~0.2 km s-1, with typical velocity centroids of other reactive molecular ions observed with single-dish telescopes toward the Orion Bar (see, e.g., Fuente et al. 2003; Leurini et al. 2006; Nagy et al. 2013).

We expected to detect the SH+F = 1/2−1/2 and F = 1/2−3/2 HFS lines of the NJ = 10 −01 FS component near 346 GHz on the basis of the previous Herschel/HIFI detection of the higher energy SH+NJ = 12 −01 lines at ~526 GHz toward the Orion Bar (Nagy et al. 2013). Blue arrows in Fig. 1 show the absence of emission at frequencies derived from previous laboratory data of Savage et al. (2004) for all three telescopes. Our spectra, however, show two emission lines shifted by ~15 MHz (~13 km s-1) to higher frequencies with respect to the frequencies indicated by blue arrows, which could not be assigned to other carriers in the MADEX, CDMS (Müller et al. 2001, 2005), and JPL (Pickett et al. 1998) spectroscopic catalogs. Table A.1 shows results from Gaussian fits to the two lines, which were detected with up to 12σ and 7σ significance levels. Averaging the three observations, we obtain 345 944.35 MHz and 345 858.27 MHz for the two lines. Given the complex gas kinematics in Orion, implying an uncertainty in vLSR, we adopted a 1σ systematic uncertainty of 0.2 MHz for the rest frequency determination. We note that the line near 345.9 GHz is only ~61 MHz (~53 km s-1) above the CO J = 3−2 line frequency.

4. New SH+ spectroscopic parameters

The two unpaired electrons (S = 1) of sulfanylium lead to a electronic ground state which causes all rotational levels N to be split into three with J = N, N ± 1, except for the N = 0 level for which only J = 1 exists. The lowest order FS parameters are λ and γ which describe the (electron) spin-spin coupling and the (electron) spin-rotation coupling. Rotational (or centrifugal distortion) and vibrational corrections, the latter for excited vibrational states, may appear in the Hamiltonian. The proton magnetic moment (IH = 1/2) splits every FS level into two HFS levels. The main HFS parameters here are bF and c, the scalar and tensorial electron spin-nuclear spin coupling parameters.

The selection rules ΔF = 0, ± 1 hold strictly. The spin-conserving transitions (ΔF = ΔJ = ΔN) are the strongest ones at higher quantum numbers. Other transitions have non-negligible intensities, in particular at low-N. The N = 1−0 ground state rotational transition is split into three FS components J = 0−1, 2−1, and 1−1 near 346, 526, and 683 GHz, respectively, which are split further into 2, 3, and 4 HFS components. An energy level diagram is shown, e.g., in Menten et al. (2011).

The rotation-vibration spectrum of SH+ was recorded up to ν= 2−1 (Brown et al. 1986) and ν= 4−3 (Civiš et al. 1989). Hovde & Saykally (1987) investigated the rotational spectrum of SH+ in its ground and first excited vibrational state between 0.51 THz and 1.62 THz by laser magnetic resonance (LMR). More recently, Savage et al. (2004) employed source- and velocity-modulation to determine accurate rest frequencies of all three HFS components near 526 GHz of the N = 1−0 transition of SH+ as well as one of two HFS components near 346 GHz. In each of these four studies, spectroscopic parameters were determined only from their own data. In contrast, Brown & Müller (2009) used the field-dependent (Hovde & Saykally 1987) and field-free (Savage et al. 2004) rotational data as well as the rovibrational data of the fundamental (ν= 1−0) band (Brown et al. 1986; Civiš et al. 1989) to determine ground state spectroscopic parameters as well as vibrational corrections if they were needed. One of the present authors (HSPM) used a slightly different approach to create an entry for the CDMS catalog because the SPFIT program (Pickett 1991) does not permit the use of field-dependent data. The LMR transition frequencies extrapolated to zero field by Brown & Müller (2009) were fit together with the field-free rotational data (Savage et al. 2004) and all of the rovibrational data (Brown et al. 1986; Civiš et al. 1989) to determine Dunham-type parameters Pi,j for SH+, with i and j indicating vibrational and rotational corrections, respectively.

Savage et al. (2004) observed all three HFS components of the NJ = 12−01 transition of sulfanylium (around 526.1 GHz) between 0.2 MHz and 0.5 MHz higher than calculated from the parameters of the LMR study (Hovde & Saykally 1987); the deviations are well within the uncertainties of ~1.6 MHz from that study. These rest frequencies improve the accuracies of the major HFS parameters bF and c, as well as the origin of the FS component, which depends on several rotational and fine structure parameters. The F = 1/2−3/2 HFS component was found by Savage et al. (2004) only 1.23 MHz higher than calculated from parameters in Hovde & Saykally (1987), well within the uncertainty of ~18 MHz. The uncertainty is much larger because this FS component was not accessed in the LMR study. The associated F = 1/2−1/2 HFS component, weaker by a factor of two, was not detected despite significant signal averaging. Its frequency, ~86 MHz lower than the F = 1/2−3/2 HFS component, is well constrained by the frequency of this HFS component and from the data near 526.1 GHz.

Table 1

Spectroscopic parametersa (MHz, cm-1) of sulfanylium, SH+.

As outlined in Sect. 3, no emission lines were detected in the Orion Bar at the frequencies expected from Savage et al. (2004). However, the two emission lines, each observed ~15 MHz higher, had the expected 1:2 intensity ratio and spacing required to assign them to the NJ = 12 −01 FS component. In order to test the feasibility of these assignments, we omitted the laboratory transition frequency near 346 GHz from the data set which was used to create the CDMS catalog entry. Interestingly, the rms error of the fit improved significantly from 0.898 to 0.826. More importantly, the two HFS components were now predicted at 345 856.8 MHz and at 345 943.0 MHz, very close to our observed emission lines (Fig. 1) and each with predicted 1σ uncertainties of 3.9 MHz. These uncertainties are considerably smaller than those from Hovde & Saykally (1987), mainly because of the extensive rovibrational transition frequencies used in the present fit and because of the remaining data from Savage et al. (2004). Hence, we used averages over the three independent observations given in Sect. 3 to redetermine the SH+ spectroscopic parameters. The resulting parameters are given in Table 1. The rms error of this fit is 0.820, demonstrating that our transition frequencies are more compatible with the remainder of the laboratory data than the single frequency near 346 GHz determined by Savage et al. (2004). Updated predictions of the rotational spectrum of SH+ will be available in the CDMS catalog3. Data up to 1 THz are given in Table A.2.

Errors in rest frequencies determined in the laboratory or through astronomical observations are known to occur especially in sparse data sets. Recent error correction of data pertaining to astrophysically important hydrides include the correction of the fundamental transition of H2D+ (Asvany et al. 2008) and of the fundamental transitions of CH+ isotopologues (Amano 2010).

The changes in spectroscopic parameters are small for the most part and occur predominantly in rotational and FS parameters. The most significant changes of 3σ occur in Y01Be and Y02 ≈ − De. The largest change in magnitude occurs in λ00, but corresponds only to ~1σ. The strong, spin-conserving transitions are only slightly affected at lower frequencies. In addition, the transitions of the NJ = 11 −01 FS component (~683 GHz) are now predicted 2 MHz lower than in the first SH+ CDMS catalog entry. Changes from the predictions by Brown & Müller (2009) are a bit more complex, but deviate on average by ~2 MHz. Larger deviations of several megahertz may occur for transitions with flip of the electron spin (e.g., at ~893 GHz) or at higher frequencies.

The 345 929.8 MHz line observed in the laboratory, and attributed to SH+, NJ = 10 −01, F = 1/2−3/2 by Savage et al. (2004), is not detected at the sensitivity limit of our observations. Here we demonstrate that the 345 929.8 MHz line is not SH+, and we improve the SH+ spectroscopic parameters with the detection of the 345 944 and 345 858 MHz lines, which arise from the edge of the Orion Bar, consistent with the expected SH+ formation route in PDRs, the reaction of S+ with vibrationally excited H2 (e.g., Nagy et al. 2013; Zanchet et al. 2013). Nevertheless, it may still be useful to reinvestigate the N = 1−0 transition of SH+ in the laboratory using methods such as those by Savage et al. (2004); Asvany et al. (2008); Amano (2010); Brünken et al. (2014), in particular in combination with a study of 34SH+.


2

We used ν(HCO+4−3) = 356 734.2244 (9) MHz. All publicly available rotational lines up to J = 17 and with uncertainties below 50 kHz (Tinti et al. 2007; Buffa et al. 2006; Hirao et al. 2008; Cazzoli et al. 2012) were used to compute HCO+ line frequencies with MADEX (Cernicharo 2012). The residual of the fit is 39 kHz.

Acknowledgments

We thank the Spanish MINECO for funding support under grants CSD2009-00038, AYA2009-07304, and AYA2012-32032. We thank the European Research Council for funding support under ERC-2013-Syg 610256-NANOCOSMOS. H.S.P.M. is supported by the German Bundesministerium für Bildung und Forschung (BMBF) via the ALMA Regional Center (ARC) Node project 5A11PK3 for maintanance and upgrade of the CDMS. J.R.G. thanks the Observatoire de Paris/ENS and IRAM (Grenoble) where part of this work was carried out. S.C. was supported by a FPI-INTA grant.

References

Online material

Appendix A: Complementary tables

Table A.1

Gaussian fits to the F = 0.5−1.5 and 0.5−0.5 lines of the NJ = 10−01 FS component of SH+ observed with three different telescopes toward the Orion Bar (for vLSR = 10.50 km s-1).

Table A.2

Quantum numbers, frequencies (MHz), uncertainties unc. (MHz), Einstein A values (10-4 s-1), upper gup and lower glo state degeneracies, and upper Eup and lower Elo state energies (cm-1) of sulfanylium, SH+, below 1 THz.

All Tables

Table 1

Spectroscopic parametersa (MHz, cm-1) of sulfanylium, SH+.

Table A.1

Gaussian fits to the F = 0.5−1.5 and 0.5−0.5 lines of the NJ = 10−01 FS component of SH+ observed with three different telescopes toward the Orion Bar (for vLSR = 10.50 km s-1).

Table A.2

Quantum numbers, frequencies (MHz), uncertainties unc. (MHz), Einstein A values (10-4 s-1), upper gup and lower glo state degeneracies, and upper Eup and lower Elo state energies (cm-1) of sulfanylium, SH+, below 1 THz.

All Figures

thumbnail Fig. 1

From top to bottom: IRAM 30 m, ACA, and ALMA spectra toward the Orion Bar between 345.8 GHz and 345.9 GHz. Blue arrows show the predicted frequencies of the SH+, NJ = 10−01, F = 1/2−1/2 and F = 1/2−3/2 lines before this work (based on experiments by Savage et al. 2004). Red arrows (prediction*) show the expected line frequencies for a Hamiltonian fit that excludes the 345 929.8 MHz line observed by Savage et al. (2004) and attributed to SH+; see Sect. 4.

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.