Free Access
Issue
A&A
Volume 533, September 2011
Article Number L11
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201117753
Published online 07 September 2011

© ESO, 2011

1. Introduction

Molecules with sulfur account for about 10% of the species identified in the interstellar gas and circumstellar envelopes. Among the most prominent of these is sulfur dioxide (SO2). Since its detection in Orion and Sgr B2 more than 35 years ago (Snyder et al. 1975), SO2 has been observed in many astronomical sources including hot cores in giant molecular clouds and low-mass young stellar objects (YSOs; Schilke et al. 2001; Comito et al. 2005; Caux et al. 2011), circumstellar shells of late-type stars (Tenenbaum et al. 2010), cold dense molecular clouds (Turner 1995), and external galaxies (Martín et al. 2003). This simple triatomic molecule is the carrier of a significant fraction of rotational lines in hot sources such as Orion and Sgr B2, which have been observed with ground-based telescopes in the millimeter band or the Herschel space observatory in the THz-band (Crockett et al. 2010; Ceccarelli et al. 2010). Despite the high density of spectral lines and wide range of astronomical sources in which it is found, the formation of SO2 in space and the chemistry of other prominent sulfur bearing molecules in the dense molecular gas such as SO, H2S, and CS is not well understood (Wakelam et al. 2004, and references therein).

In addition to prominent neutral sulfur-bearing molecules, several positive molecular ions containing sulfur are also observed. Of these, SO+ is widespread and conspicuous in over 20 galactic sources both warm and cold (Turner 1994), HCS+ has a similarly wide distribution, and SH+ was recently shown to be ubiquitous in the diffuse interstellar gas (Menten et al. 2011). Although protonated SO2 (HOSO+) has not explicitly been considered in most chemical networks, it is an excellent candidate for astronomical detection for several reasons: (1) it does not react with H2; (2) the proton affinity of SO2 is high (de Petris et al. 2011); (3) it is formed facilely by the reaction of H+3\hbox{$_3^+$} with SO2 (Anicich & Huntress 1986); (4) the dipole moment is large (μa = 1.74 D; Lattanzi et al. 2010); (5) the parent SO2 is very abundant in hot regions; (6) and radio observations of SH+ and SO+ (Menten et al. 2011; Turner 1994) have established that positive ions of sulfur-bearing molecules are surprisingly abundant with respect to their neutral counterparts and other common sulfur species in many molecular sources.

Here we report the millimeter-wave rotational spectrum of cis-HOSO+, the most stable isomeric form of this ion. The present measurements improve the accuracy of the radio spectrum by two orders of magnitude over that obtained from the earlier laboratory measurements of the lowest rotational transitions in the centimeter band (Lattanzi et al. 2010), allowing a deep search for this fundamental positive ion in regions of the interstellar gas where there is an active sulfur chemistry. From the derived spectroscopic constants, the astronomically most interesting lines of HOSO+ up to 700 GHz can be predicted to 4 km s-1 or better in equivalent radial velocity.

2. Experiment and analysis

The rotational spectrum of protonated SO2 (HOSO+) was first observed at centimeter wavelengths in a supersonic molecular beam with a Fourier transform microwave (FTM) spectrometer used to detect many new reactive radicals and ions of astrophysical interest (Lattanzi et al. 2010, and references therein). The HOSO+ ion was produced in a small nozzle by a low-current DC discharge through SO2 heavily diluted in H2. Following the detection of the four lowest a-type transitions, two b-type cross-ladder (Ka = 1 → 0) transitions in the centimeter band were observed in double resonance, thereby confirming the spectroscopic identification and yielding a precise determination of all three rotational constants.

Millimeter-wave rotational lines of the ground state (cis) isomer of HOSO+ were observed in a negative glow discharge cell (7.6 cm diameter, 1.4 m long) similar in design to that originally described by De Lucia et al. (1983). The experimental conditions were optimized by first observing the known positive ions HCO+, SO+, and HCS+ whose rotational lines all exhibited enhancements in intensity of 50 or more upon application of a strong axial magnetic field of 170–200 G. The most intense lines of HOSO+ (e.g., Fig. 1) were observed in a DC discharge (3 mA) through a flowing mixture of Ar, SO2, and H2 with the discharge cell cooled to 245 K. The total pressure with the discharge running (3.5 mTorr) was maintained by a small diffusion pump (800 l/s pumping speed). Cooling of the discharge cell was accomplished by flowing liquid nitrogen through copper tubing in thermal contact with the solenoid, with the space between the glass cell and the aluminum solenoid form filled with aluminum wool.

thumbnail Fig. 1

The J = 150,15−140,14 rotational transition of cis-HOSO+, illustrating the effect of an applied magnetic field on the line intensity in the negative glow discharge. Owing to the modulation and detection scheme employed, the instrumental lineshape is approximately the second derivative of a Lorentzian. Each spectrum is the result of 12 min of integration.

The search for HOSO+ in the millimeter band was guided by the rotational constants that have been derived from the centimeter-wave data and theoretical centrifugal distortion constants (Lattanzi et al. 2010). Owing to the bent heavy atom backbone of this ion (Fig. 2), the A rotational constant is fairly small (44 GHz) and B − C is large (1.8 GHz). As a result, the departure from harmonicity of a-type transitions in the Ka = 0 ladder is large (e.g., 555 MHz at 298.7 GHz). Thus, when the first suspected millimeter-wave rotational line of HOSO+ was observed, the identification could not be confirmed by simply observing other harmonically related lines. Instead, our search relied on frequencies calculated with the experimental and theoretical spectroscopic constants. Recent work on rigid molecules with two or three heavy atoms established that the fourth-order distortion constant DJ calculated at a high level of theory is usually accurate to about 3% or better (McCarthy et al. 2011; Spezzano et al. 2011). We found that subsequent lines of HOSO+ could be predicted to within about 5 MHz near 300 GHz, because the contribution from DJ to transitions in the Ka = 0 ladder is small. The departure from harmonicity can be calculated to fairly high accuracy from previous data.

Five lines between 186 and 347 GHz in the Ka = 0 ladder of the cis isomer of HOSO+ were observed in our negative glow discharge, each exhibiting a small enhancement in intensity when an axial magnetic field was applied (Table 1). The observed enhancement of between a factor of 3−5 with a field strength of  ~170 G was sufficient to allow us to discriminate between ions and neutral species produced in the discharge (Fig. 1). Although the magnitude of the magnetic field enhancement was about 10 times smaller than that of other molecular ions observed with this spectrometer (i.e., HCO+, HOCO+, and HCS+), the effect of the magnetic field was the same for all millimeter-wave lines assigned to HOSO+ and was repeatable. The identification of lines of HOSO+ observed here in the millimeter band was supported by the detection of four lines in the Ka = 1 ladder that were close in frequency to the prediction.

thumbnail Fig. 2

Relative intensities of the a-type spectrum of cis-HOSO+ at 50 K (top) and 200 K (bottom). The intensities were calculated on the assumption that the populations of all rotational levels are described by a single temperature, and only the ground vibrational state is populated. Shown in the insert (top) is the molecular structure from Lattanzi et al. (2010).

The rotational spectrum of cis-HOSO+ was analyzed with the Watson S-reduced Hamiltonian with three rotational and two centrifugal distortion constants (DJ and DK), and the remaining three fourth-order distortion constants constrained to the theoretical values. With seven spectroscopic constants (Table 2), the ten-millimeter-wave and seven-centimeter-wave lines were reproduced with an rms uncertainty (19 kHz) that is well within the measurement uncertainties.

Table 1

Measured rotational transitions of cis-HOSO+.

Table 2

Spectroscopic constants of cis-HOSO+ (in MHz).

3. Discussion

The high observed abundance of the SH+ ion with respect to that of CS, SO, SO2, and H2S (0.5 to  >55; Menten et al. 2011), and the similarly high ratio of SO+ to that of SO and SO2 (~1; Turner 1994) establish that positive ions of sulfur-bearing molecules are conspicuous in the interstellar gas. On the assumption that HOSO+ is formed in space by the reaction of SO2 with H3+\hbox{${\rm H}_3^+$}, the best sources for detecting this ion are those with the highest abundance of SO2, such as hot cores in high-mass star-forming regions and low-mass YSOs.

Many rotational transitions of HOSO+ of comparable predicted intensity are accessible throughout the millimeter- and THz-bands, facilitating a deep search with ground-based and space-borne telescopes. Figure 2 shows a simulated spectrum of HOSO+ calculated for a rotational temperature (Trot) of 200 K appropriate for a typical hot core region, on the assumption that all rotational levels are described by a single temperature. The dominant features are strong R-branch (ΔJ = 1) transitions separated by 16 GHz, with transitions in the Ka = 1 and Ka = 2 ladders nearly as intense as those in the Ka = 0 ladder owing to the fairly small value of the A rotational constant. As this plot illustrates, the THz band is the most favorable spectral region for a radioastronomical search in sources with Trot near 200 K. Also shown in Fig. 2 is the spectrum calculated for a warm region with Trot = 50 K, such as that between the hot corino in IRAS 16293-2442 and the cool (~10 K) surrounding dust envelope (Caux et al. 2011). At this intermediate temperature, the peak line-intensity shifts to the millimeter-band, which is readily accessible with ground based telescopes.

The spectroscopic constants of HOSO+ (Table 2) allow an accurate prediction of transitions up to the 600−700 GHz region (Band 9 of ALMA). The sixth-order distortion constant HJ was not determined in the measurements here, but on the assumption that it is similar to that of SO2 (HJ = 1.1 × 10-8 MHz; Müller, & Brünken 2005), we estimate that the uncertainty in the calculated frequencies near 650 GHz owing to neglect of HJ (~4 km s-1) is comparable to the linewidths in hot cores and corinos.

We find no evidence for HOSO+ in published spectral line surveys in Orion where the column density of SO2 is high (5 × 1016 cm-2; Schilke et al. 2001; Comito et al. 2005), or in Sgr B2 where the abundance of SO2 is at least 20 times higher (≥1 × 1018 cm-2; Nummelin et al. 2000). On the assumption that the ion is present in the same region as SO2, the column density of HOSO+ in Orion (≤3 × 1015 cm-2) implies that the upper limit of the HOSO+/SO2 ratio (≤6%) is comparable to the observed 34SO2/SO2 ratio in this source. Interferometric observations of the hot corino in IRAS 16293-2422 with the SMA confirm that the abundances of sulfur-bearing molecules in this source are high as well (2.9 × 1016 cm-2; Chandler et al. 2005). Chandler et al. covered a range of about 3 GHz near 305 GHz, which unfortunately does not encompass a strong transition of HOSO+, so we are unable to determine whether there is any evidence for this ion in the corino. The HOSO+ ion might be detectable in IRAS 16293-2422 in a dedicated search with the SMA if the abundance of HOSO+ is comparable to that of 34SO2, which is observed with a signal-to-noise of 2−3. The combination of high sensitivity and high spatial resolution should greatly advance the study of corinos and the search for this protonated species when ALMA comes online.

The rotational spectrum of the trans isomer of protonated sulfur dioxide, calculated to lie only 2.4 kcal/mol or  ~1200 K higher in energy, has also been measured in the centimeter-band (Lattanzi et al. 2010). If this isomer is produced in our negative glow discharge in comparable abundance to that of the cis isomer, its millimeter-wave rotational lines should also be observable because of its higher dipole moment (μa = 3.16 D versus 1.74 D for cis-HOSO+). Measurements of both the cis and trans isomers in a molecular beam with our FTM spectrometer imply that this might be the case. By analogy with the cis isomer reported here, a laboratory search for the millimeter-wave transitions of trans-HOSO+ should be undertaken now that the rotational constants have been derived from the rotational data, and the leading centrifugal distortion constants were calculated at the same level of theory as that of the cis isomer (Lattanzi et al. 2010). It is anticipated that an extended search in frequency will not be required for the identification of the millimeter-wave lines of trans-HOSO+.

Other small sulfur-bearing molecules of astrophysical interest whose rotational spectra have not been measured in the millimeter and THz band include the HOSO radical and the HSO+ ion. The HOSO radical is an excellent candidate for detection in space. In some recent calculations (Garrido et al. 2011; Ballester et al. 2010), full dimensional quasi-classical trajectory studies of OH+SO reaction were undertaken. They find that (1) the reaction is exothermic, with an enthalpy of formation of −29.4 kcal/mol; (2) there is no barrier in the entrance channel to form the HOSO intermediate; (3) HOSO represents the global minimum on the potential energy surface, lying 71.8 kcal/mol below the OH + SO limit; and (4) HOSO is the only product that is formed in the energy range studied. Recently rotational lines of the radical have been observed in the centimeter-band and the leading distortion constants have been determined by high-level quantum chemical calculations (Lattanzi et al. 2011).

Protonated SO (HSO+) is an especially attractive candidate for detection: it is predicted to have a closed-shell A1\hbox{$^1A^\prime$} electronic ground state (Li et al. 2006); the dipole moment is high (equilibrium values calculated at the CCSD(T)/cc-pwCVQZ level of theory are μa = 2.9 D and μb = 1.1 D; this work); and it has a rich spectrum in the millimeter and THz bands. With only two heavy atoms the lowest rotational transition of HSO+ lies outside the range of our FTM spectrometer, but the same experimental techniques used here for cis-HOSO+ can be applied to the laboratory search for this ion in the millimeter-band. The column density of HSO+ is predicted to be high in the envelopes around massive YSOs (~1015 cm-2; Stäuber et al. 2005), owing to the considerable predicted enhancement in the envelope by X-rays. Therefore, the HSO+ ion might be an important tracer of X-rays in these objects when its rotational spectrum has been measured in the laboratory.

Acknowledgments

The authors wish to thank N. Adams and V. Bierbaum for helpful discussions on ion chemistry; and E. Palmer, D. Kokkin, and O. Martinez for assistance with the ion cell. This work was supported in part by NASA Grant NNX08AE05G. S.T. gratefully acknowledges support by the Deutsche Forschungsgemainschaft (DFG) through grant TH 1301/3-1.

References

  1. Anicich, V. G., & Huntress, Jr., W. T. 1986, ApJS, 62, 553 [NASA ADS] [CrossRef] [Google Scholar]
  2. Ballester, M. Y., Orozco-Gonzalez, Y., Garrido, J. D., & Santos, H. F. D. 2010, 132, 044310 [Google Scholar]
  3. Caux, E., Kahane, C., Castets, A., et al. 2011, A&A, 532, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Ceccarelli, C., Bacmann, A., Boogert, A., et al. 2010, A&A, 521, L22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Chandler, C. J., Brogan, C. L., Shirley, Y. L., & Loinard, L. 2005, ApJ, 632, 371 [NASA ADS] [CrossRef] [Google Scholar]
  6. Comito, C., Schilke, P., Phillips, T. G., et al. 2005, ApJS, 156, 127 [NASA ADS] [CrossRef] [Google Scholar]
  7. Crockett, N. R., Bergin, E. A., Wang, S., et al. 2010, A&A, 521, L21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. De Lucia, F. C., Herbst, E., Plummer, G. M., & Blake, G. A. 1983, J. Chem. Phys., 78, 2312 [NASA ADS] [CrossRef] [Google Scholar]
  9. de Petris, G., Cartoni, A., Rosi, M., et al. 2011, Chem. Phys. Chem., 12, 2 [Google Scholar]
  10. Garrido, J. D., Ballester, M. Y., Orozco-Gonzalez, Y., & Canuto, S. 2011, J. Phys. Chem. A, 115, 1453 [CrossRef] [Google Scholar]
  11. Lattanzi, V., Thaddeus, P., McCarthy, M. C., & Thorwirth, S. 2010, J. Chem. Phys., 133, 194305 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  12. Lattanzi, V., McCarthy, M. C., & Tamassia, F. 2011, in 66th Ohio State University International Symposium on Molecular Spectrosopy (Columbus: Ohio State Univ.), Talk WJ02 [Google Scholar]
  13. Li, B.-T., Wei, Z.-Z., Zhang, H.-X., & Sun, C.-C. 2006, J. Phys. Chem. A, 110, 10643 [CrossRef] [Google Scholar]
  14. Martín, S., Mauersberger, R., Martín-Pintado, J., García-Burillo, S., & Henkel, C. 2003, A&A, 411, L465 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. McCarthy, M. C., Gottlieb, C. A., Thaddeus, P., Thorwirth, S., & Gauss, J. 2011, J. Chem. Phys., 134, 034306 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  16. Menten, K. M., Wyrowski, F., Belloche, A., et al. 2011, A&A, 525, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Müller, H. S. P., & Brünken, S. 2005, J. Mol. Spec., 232, 213 [CrossRef] [Google Scholar]
  18. Nummelin, A., Bergman, P., Hjalmarson, Å., et al. 2000, ApJS, 128, 213 [NASA ADS] [CrossRef] [Google Scholar]
  19. Schilke, P., Benford, D. J., Hunter, T. R., Lis, D. C., & Phillips, T. G. 2001, ApJS, 132, 281 [NASA ADS] [CrossRef] [Google Scholar]
  20. Snyder, L. E., Hollis, J. M., Ulich, B. L., et al. 1975, ApJ, 198, L81 [NASA ADS] [CrossRef] [Google Scholar]
  21. Spezzano, S., Tamassia, F., Thorwirth, S., et al. 2011, ApJS, submitted [Google Scholar]
  22. Stäuber, P., Doty, S. D., van Dishoeck, E. F., & Benz, A. O. 2005, A&A, 440, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Tenenbaum, E. D., Dodd, J. L., Milam, S. N., Woolf, N. J., & Ziurys, L. M. 2010, ApJS, 190, 348 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  24. Turner, B. E. 1994, ApJ, 430, 727 [NASA ADS] [CrossRef] [Google Scholar]
  25. Wakelam, V., Caselli, P., Ceccarelli, C., Herbst, E., & Castets, A. 2004, A&A, 422, 159 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Measured rotational transitions of cis-HOSO+.

Table 2

Spectroscopic constants of cis-HOSO+ (in MHz).

All Figures

thumbnail Fig. 1

The J = 150,15−140,14 rotational transition of cis-HOSO+, illustrating the effect of an applied magnetic field on the line intensity in the negative glow discharge. Owing to the modulation and detection scheme employed, the instrumental lineshape is approximately the second derivative of a Lorentzian. Each spectrum is the result of 12 min of integration.

In the text
thumbnail Fig. 2

Relative intensities of the a-type spectrum of cis-HOSO+ at 50 K (top) and 200 K (bottom). The intensities were calculated on the assumption that the populations of all rotational levels are described by a single temperature, and only the ground vibrational state is populated. Shown in the insert (top) is the molecular structure from Lattanzi et al. (2010).

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.