A&A 493, 565-569 (2009)
DOI: 10.1051/0004-6361:200810889

Rotational spectrum of deuterated and 15N ethyl cyanides: CH3CHDCN and CH2DCH2CN and of CH3CH2C15N[*],[*]

L. Margulès1 - R. Motiyenko1 - K. Demyk2 - B. Tercero3 - J. Cernicharo3 - M. Sheng4 - M. Weidmann4 - J. Gripp4 - H. Mäder4 - J. Demaison1


1 - Laboratoire de Physique des Lasers, Atomes et Molécules, UMR CNRS 8523, Université Lille 1, 59655 Villeneuve d'Ascq Cedex, France
2 - Centre d'Étude Spatiale des Rayonnements, UMR CNRS 5187, Université de Toulouse (UPS), 9 avenue du Colonel Roche, BP 4346, 31028 Toulouse, France
3 - Laboratory of Molecular Astrophysics, CAB-CSIC. INTA, Ctra de Torrejón a Ajalvir, km 4, 28850 Torrejón de Ardoz, Madrid, Spain
4 - Universität Kiel, Institut für Physikalische Chemie Olshausenstrasse 40, 24098 Kiel, Germany

Received 1 September 2008 / Accepted 10 October 2008

Abstract
Context. Ethyl cyanide is an abundant molecule in hot molecular clouds. Its rotational spectrum is very dense and several hundreds of rotational transitions within the ground state have been identified in molecular clouds in the 40-900 GHz frequency range. Lines from 13C isotopically substituted ethyl cyanide were identified in Orion.
Aims. To enable the search and the possible detection of other isotopologues of ethyl cyanide in interstellar objects, we have studied the rotational spectrum of deuterated ethyl cyanide: CH2DCH2CN (in-plane and out-of-plane) and CH3CHDCN and the spectrum of 15N substituted ethyl cyanide CH3CH2C15N. Using these experimental data, we have searched for these species in Orion.
Methods. The rotational spectrum of each species in the ground state was measured in the microwave and millimeter-submillimeter wavelength range using a waveguide Fourier Transform spectrometer (8-17 GHz) and a source-modulated spectrometer employing backward-wave oscillators (BWOs) (150-260 and 580-660 GHz). More than 300 lines were identified for each species, for J values in the range 71-80 and $K_{\rm a}$ values in the range 28-31 depending on the isotopologues. The experimental spectra were analyzed using a Watson's Hamiltonian in the A-reduction.
Results. From the fitting procedure, accurate spectroscopic constants were derived for each of the species. These new sets of spectroscopic constants enable us to predict reliably the rotational spectrum (lines frequencies and intensities) in the 4-1000 GHz frequency range and for J and $K_{\rm a}$ up to 80 and 31, respectively. Combined with IRAM 30 m antenna observations of Orion, this experimental study allowed us to detect 15N substituted ethyl cyanide CH3CH2C15N for the first time in Orion. The derived column density and rotational temperature are 1013 cm-2 and 150 K for the plateau and 3 $\times$ 1014 cm-2 and 300 K for the hot core. The deuterated species were searched for but were not detected. The upper limit to the column density of each deuterated isotopologues was 1014 cm-2.

Key words: line: identification - methods: laboratory - molecular data - ISM: molecules - radio lines: ISM - submillimeter

1 Introduction

Ethyl cyanide, CH3CH2CN, is an asymmetric top molecule with a large dipole moment ( $\mu_{\rm a}$ = 3.85 D and $\mu_{\rm b}$ = 1.23 D) that exhibits a dense and intense rotational spectrum. It is present in the densest parts of hot molecular cores, where it is proposed to form on dust grains. Several hundreds of lines of CH3CH2CN in the ground state have been observed towards hot cores such as Orion, Sgr B2 and W51 (Liu et al. 2001; Miao & Snyder 1997) but also toward low mass star-forming regions (Cazaux et al. 2003; Remijan & Hollis 2006). It has a high abundance of the order of 1015-1017 cm-2 depending on the sources (Miao & Snyder 1997; Remijan et al. 2004; Remijan & Hollis 2006). Transitions from vibrationally excited ethyl cyanide have also been observed in Sgr B2 (Mehringer et al. 2004) and in W51 e2 (Demyk et al. 2008). Numerous lines from 13C-substituted ethyl cyanide have been detected in Orion Irc2 (Demyk et al. 2007).

All of these observations show that the unidentified lines observed in spectral surveys of molecular clouds are partly due to transitions from known species in vibrationally excited states or from isotopologues of known species. The most promising carriers of these transitions are the so-called interstellar weeds, i.e. molecules, such as ethyl cyanide, which have a dense and intense rotational spectrum and/or low-frequency vibrational modes, such as methyl formate HCOOCH3, dimethyl-ether CH3OCH3, or methanol CH3OH. With the increase in sensitivity and frequency coverage that will be achieved with instruments such as HIFI onboard the Herschel Space Observatory and ALMA, the identification of these U-lines will become crucial to the search for new molecules but also to obtain important information about the physical and chemical conditions prevailing in the observed sources.

An enormous amount of experimental work must be undertaken to complete the actual knowledge of the rotational spectra of low-frequency vibrational mode and of the isotopologues of abundant interstellar molecules. When they exist, measurements are indeed often limited to low frequencies and therefore cannot be used to predict reliable line frequencies and intensities in the Herschel and ALMA frequency ranges.

In this context, following up our study on 13C-substituted ethyl cyanide (Demyk et al. 2007), we present an experimental study of the rotational spectrum of deuterated ethyl cyanide, CH2DCH2CN in-plane, CH2DCH2CN out-of-plane, CH3CHDCN, and 15N substituted ethyl cyanide, CH3CH2C15N. These species were studied 30 years ago by Mäder et al. (1973) and Heise et al. (1976). For each species, about 30 low-J (J $\leq$ 6) rotational transitions were recorded and assigned in the 4-40 GHz frequency range. The rotational constants and quadrupole coupling constants derived from these studies were used to extend the measurements to higher frequency. The isotopologues synthesis and the experimental setup are described in Sect. 2. The analysis of the measured rotational transitions, the resulting spectroscopic parameters, and the prediction of rotational spectrum for each isotopologue in the 8-1000 GHz range are presented in Sect. 3. The detection of CH3CH2C15N in Orion and the search for the deuterated species is presented in Sect. 4.

   
2 Experimental setup

The deuterated ethyl cyanides, CH2DCH2CN and CH3CHDCN, have been prepared from the corresponding ethyl iodides CH2DCH2I and CH3CHDI (98 atom%, C/D/N isotopes, Pointe Claire, Canada) and normal potassium cyanide. KC15N (98 atom%, Aldrich, Taufkirchen, Germany) and normal ethyl iodide were used for the preparation of the 15N-isotopologue. The potassium cyanide was dispersed in a solution of ethyl iodide in triethylene glycol and the mixture was stirred and heated slowly to 110 $^{\circ}$C. Potassium cyanide and ethyl iodide were used in a molar ratio 1.25: 1.00, and the concentration of the ethyl iodide solution was 4 mol/l. Nitrogen was bubbled through the mixture and the ethyl cyanide was isolated in an ice-cooled trap. The yield of this modified Kolbe reaction (Mäder et al. 1973; Organikum 1977) was 75%. The product was distilled and controlled spectroscopically. The only detectable impurity was a small amount of the corresponding isonitrile.

In Kiel, the measurements in the centimeter-wave range were performed by means of waveguide Fourier-transform microwave-spectroscopy (Sarka et al. 1997). A spectrometer in the range of about 8-17 GHz was used, employing an oversized X-band sample cell with a rectangular waveguide of quadratic cross-section and 12 m length (Krüger et al. 1993). The experiments were carried out at ambient temperature and at gas pressures of ca. 0.1 Pa. Experimental transition frequencies were obtained from an analysis of the frequency-domain signals; these were derived following Fourier transformation of the transient emission signal, using a peak finder routine that determined line-center frequencies to an accuracy typically superior to 20 kHz, depending on the strength of the line. In the case of observed line splittings due to methyl internal rotation (AE-splittings) and/or the 14N-nuclear quadrupole coupling, the experimental peak frequencies were corrected to obtain hypothetical unsplit line frequencies (without nuclear quadrupole hyperfine structure).

The millimeter spectra were recorded in Lille in the spectral range 150-250 and 580-660 GHz. The sources were Russian Istok backward-wave oscillators (BWO). They were phase locked on an harmonic from a HP synthesizer. Up to 250 GHz, the signal from the synthesizer was directly mixed onto a Russian planar Schottky diode with part of the signal from the BWO. From 500 GHz to 650 GHz, an active sextupler from millitech (75-100 GHz) and a Schottky planar diode placed in a parabolic structure (from Virginia Diodes Inc.) optimized for this range were used. The detector was an InSb liquid He-cooled bolometer from QMC. To improve the sensitivity of the spectrometer, the sources were frequency-modulated at 5 kHz. The absorption cell was a stainless steel tube (6 cm diameter, 110 cm long), and the pressure that we used during measurements was 2.6 Pa (26 $\mu$bar). The accuracy of isolated lines was superior to 30 kHz.

Table 1: Spectroscopic measurements of the ground vibrational state of CH3CHDCN, CH2DCH2CN and CH3CH2C15N.

Table 2: Spectroscopic constants of the ground vibrational state of CH3CHDCN, CH2DCH2CN and CH3CH2C15N.

   
3 Spectral analysis and line predictions

Ethyl cyanide and its isotopologues are prolate asymmetric top molecules. The main isotopologue, CH3CH2CN, has a large dipole moment ( $\mu_{\rm a}$ = 3.83 D, $\mu_{\rm b}$ = 1.23 D; Heise et al. 1976), which was used to calculate the lines intensities of the four studied isotopologues, since it was found that the rotation of the principal axes system upon isotopic substitution does not induce significant variation in the dipole moment. The experimental spectrum of each isotopologue is very dense and intense. It contains lines from the main isotopologue, which may be present as a trace in the samples. It also exhibits lines from the first low-frequency vibrationally excited states (the CH3 torsion mode and the CCN bending mode). Consequently, some lines in the measured spectra are blended or distorded and are therefore not used in the analysis. 14N-nuclear quadrupole coupling and the internal rotation of the CH3 group introduce splitting of the lines. However, these effects were not taken into account in the analysis, since these splittings are not observed in the millimeter range, and in the microwave region the unsplit line frequencies are used when splitting is observed.

The spectral analysis is a step-by-step process with permanent interaction between measurements and theory. First of all, we used the spectroscopic parameters derived from previous experimental studies at low frequency (4-40 GHz) to provide a first prediction of the line positions at low frequency and for low J and $K_{\rm a}$ value for each species. We used the experimental work from Heise et al. (1976) and Mäder et al. (1976) for CH3CH2C15N and CH2DCH2CN, respectively. The quartic distortional constants missing ( $\delta_{J}$ and  $\delta_{K}$) were fixed to the values of the normal species. In the absence of any experimental data, the rotational constants of CH3CHDCN, were calculated using ab initio calculation at the level B3LYP/6-31G* with Gaussian 03 (Frisch et al. 2004). Based on these predictions, new lines were identified in new measurements performed in the 8-17 GHz range. All of the measured lines were then fitted to derive a new and more accurate set of spectroscopic constants. These constants were then used to derive a new prediction at higher frequency and for a higher value of J and $K_{\rm a}$. New lines were measured and added to the fit step by step until the fit of the measured lines and the precision of the derived spectroscopic parameters were sufficiently good for a reliable prediction to be made for the desired frequency range and value of quantum numbers. For the fitting procedure of the measured lines and for the predictions, we have used a Watson's Hamiltonian using A-reduction in $I^{\rm r}$ representation (Watson 1977). The S reduction was also attempted but without significant improvement (as for the 13C species, Demyk et al. 2007). The molecular parameters were determined by fitting the experimental frequencies using the iteratively re-weighted least squares fitting method (Bakri et al. 2002; Hamilton 1992). The objective of this method was to derive suitable weights using the residuals of a previous iteration. It had the advantage of being more robust than the standard least squares fitting methods and automatically rejecting most of the misassigned lines.

The number of measured lines, the maximum value of J and $K_{\rm a}$, and the standard deviation of the fits are presented in Table 1 for the four studied species: CH3CH2C15N, CH2DCH2CN in-plane, CH2DCH2CN out-of-plane, and CH3CHDCN. For the CH3CH2C15N and CH3CHDCN species, it is possible to reduce the standard deviation of the fit slightly by using higher order centrifugal distortion constants; but they are only marginally determined and they worsen the precision of the predictions. They were therefore no longer considered. The spectroscopic parameters derived from the fitting procedure are presented in Table 2. The list of measured lines are accessible in the online section as Tables 3-6 for CH3CH2C15N, CH3CHDCN, CH2DCH2CN in-plane, and CH2DCH2CN out-of-plane, respectively. For each measured line, the tables indicate its assignment (quantum numbers), the measured frequency, the difference between the observed and calculated frequency, the line strength, the dipole component, and the energy of the lower level.


  \begin{figure}
\par\includegraphics[width=16cm,clip]{0889fig1.eps}
\end{figure} Figure 1: 15N-ethyl cyanide isotopologues detection in Orion. The spectra are in units of main beam temperature ( $T_{\rm mb}$). The histogram spectra offset with respect to each other are the observations (black curve) and the model (smooth red line).
Open with DEXTER

The frequency range of the measurements (8-660 GHz) and the value of the quantum numbers of the identified lines (J $\le$ 80 and $K_{\rm a}$ $\le$ 31) are suitable for astronomical studies. Ethyl cyanide is present in hot cores and has a rather high temperature, in the 100-300 K range, as its isotopologues should also have. For such temperatures (100 and 300 K), the most intense rotational transitions of ethyl cyanide occur around 238 and 407 GHz, respectively, corresponding to J values of 25 and 45, respectively. The set of spectroscopic parameters derived for each species thus allows us to predict reliably the line frequencies and the band intensity of the transitions in the spectral range suitable for interstellar detection. For each species we have calculated a prediction of the rotational spectrum in the 8-1000 GHz range for J $\le$ 80 and $K_{\rm a}$ $\le$ 31. A short sample of the predictions is shown in the online section for each isotopologues (Tables 7-10), the entire tables are available in electronic form at the CDS. The tables indicate the quantum numbers of the transition, the calculated frequency and uncertainty, the line strength, the dipole component and the energy of the lower level. The calculated error (third column in the tables) is estimated from the accuracy of the spectroscopic parameters derived by the fitting procedure. However, to get a more realistic estimation of the error, it must be multiplied by a factor 3 for the strongest lines to 10, for the weakest lines. The precision on the line frequency is good enough for line identification in interstellar spectra. For the lines that are the most suitable for detection, i.e., the strongest lines having J value up to $\sim$50-60 and a low $K_{\rm a}$ value, the error on the predicted frequencies is a few hundred kHz. The error is larger for the weakest lines and increases as J and $K_{\rm a}$ become larger, i.e., as the frequency increases.

   
4 CH3CH2C15N detection in Orion

4.1 Astronomical observation

The observations were carried out using the IRAM 30 m radio telescope during 2004 September (3 mm and 1.3 mm), 2005 March (2 mm), 2005 April (3 mm and 1.3 mm). We acquired data for the entire spectral range detectable by the 30-m recievers. The four SiS receivers operating at 3, 2 and 1.3 mm were used simultaneously. Each receiver was tuned to a single sideband with image rejections within 20-27 dB (3 mm receivers), 12-16 dB (2 mm receivers), and 13 dB (1.3 mm receivers).

System temperatures were 100-350 K for the 3 mm receivers, 200-500 K for the 2 mm receivers, and 200-800 K for the 1.3 mm receivers, depending on the particular frequency, weather conditions, and source elevation. The intensity scale was calibrated using two absorbers at different temperatures and using the Atmospheric Transmission Model (Cernicharo 1985; Pardo et al. 2001).

Pointing and focus were regularly monitored by observing the nearby quasars 0420-014 and 0528+134. The observations were completed in the balanced wobbler-switching mode with a wobbling frequency of 0.5 Hz and a beam throw in azimuth of $\pm$240''. The backends used were two filter banks with 512 $\times$ 1 MHz channels and a correlator providing two 512 MHz bandwidths and 1.25 MHz resolution. We performed a spectral-line survey, for which the central frequencies were chosen in a systematic way: from 80 GHz to 115.5 GHz for the 3 mm domain; from 130.25 GHz to 176.75 GHz for 2 mm; from 197 to 241 GHz for 1.3 mm (low frequency) and from 241.25 to 281.75 GHz for the 1.3 mm domain (high frequency), in steps of 500 MHz. We pointed toward the (survey) position $\alpha$ = 5$^{\rm h}$35$^{\rm m}$14.5$^{\rm s}$, $\delta$ = -5$^{\circ}$22'30.0'' (J2000.0), corresponding to IRc2. The detailed procedure used for the analysis of the line survey is described in Tercero et al. (in preparation).

4.2 Astronomical modeling

In agreement with previous observations of Orion, four clearly defined kinematic regions with quite different physical and chemical conditions (Blake et al. 1987,1996) are implied by the observed LSR velocities and line widths: (i) the narrow ($\la$5 km s-1 line width) feature at $v_{{\rm LSR}}$ $\simeq$ 9 km s-1, forming a N-S extended ridge or ambient cloud; (ii) a compact and quiescent region, compact ridge, ( $v_{{\rm LSR}}$ $\simeq$ 8 km s-1, $\Delta v$ $\simeq$ 3 km s-1), identified for the first time by Johansson et al. (1984); (iii) the more turbulent and compact plateau ( $v_{{\rm LSR}}$ $\simeq$ 6-10 km s-1, $\Delta v$ $\ga$ 25 km s-1); (iv) the hot core component ( $v_{{\rm LSR}}$ $\simeq$ 3-5 km s-1, $\Delta v$ $\la$ 10-15 km s-1), first observed in ammonia emission (Morris et al. 1980).

In modeling the emission from the 15N isotopologue of ethyl cyanide, we found that, as for the 13C isotopologues (Demyk et al. 2007), a sum of two components is sufficient to reproduce all line intensities and profiles reasonably well: the hot core component and the plateau. We assumed LTE for both components. For the hot core, a column density of 3 $\times$ 1014 cm-2, a line width of 5 km s-1, and a rotational temperature of 300 K were the most suitable parameters for reproducing the bulk of the CH3CH2C15N emission. A broad velocity component is also required to reproduce the observations accurately, corresponding to the plateau for which we obtain a column density of 1 $\times$ 1013 cm-2, a line width of 20 km s-1, and a rotational temperature of 150 K. For the hot core component, we assumed a source diameter of 7'' with uniform brightness temperature and optical depth over its extent at 3'' from the pointed position (the observation were pointed towards IRc2, while the CN bearing species appears to originate in a small region 3'' North). For the plateau component, we assumed a size for the source of 30''.

The lines of CH3CH2C15N are weak and many of them are heavily blended with lines from other species. However, no missing lines were found in the coverage of our line survey of Orion and the lines of CH3CH2C15N free of blending appear at the correct frequencies and with the correct intensities. Although our modeling lacks a robust analysis, the difference between model and observed intensities is always below 20%. Our modeling of the different components is straightforward, although we emphasize that we modeled the main isotope, for which strong lines are observed, with the same model and the lines are well reproduced. Hence, we are confident about the results for the isotope 15N, in particular that the frequency of all detected lines differ from laboratory measurements by less than 0.5 MHz. Hence, our assignment of the lines shown in Fig. 1 and in Table 11 to the 15N isotopologue is fairly secure. The procedure that we used to identify the carriers of the weak lines in Orion, many of which remain unassigned at present, is the most suitable for this source since it permits to confirm whether there are no missing lines of the molecules for which we are looking.

Figure 1 shows selected observed lines of the 15N isotopologue with model results. Table 11 (see online material) indicates the model predictions, observed peak intensities and frequencies, and predicted frequencies from the rotational constants obtained in this paper, for all lines of CH3CH2C15N that were not strongly blended with other lines. The differences between the intensity of the model and the peak intensity of the observed lines were mostly due to the contribution from many other molecular species (the strong overlap with other lines ensures that it is difficult to provide a good baseline for the weak lines of CH3CH2C15N).

Using the column density derived for the 13C isotopologues (1.6 $\times$ 1015 and 6 $\times$ 1013 cm-2 for the hot core and the plateau, respectively; see Demyk et al. 2007), we derive an isotopic ratio 13C/15N of between 5-6, in agreement with the solar isotopic abundance (13C/15N (solar) $\simeq$ 6) and strengthening our identification of 15N-ethyl cyanide in Orion. The detailed modeling of ethyl cyanide including the main isotopologue, the vibrationally excited states and the detected isotopologues, will be published elsewhere (Tercero et al., in preparation).

Deuterated ethyl cyanide has not been detected in this line survey above the line confusion limit. Assuming the same physical conditions as those derived for CH3CH2C15N, we have derived an upper limit to the column densities of CH2DCH2CN (in-plane and out-of-plane) of 1 $\times$ 1014 cm-2 in both species.

Acknowledgements
This work was supported by the Programme National ``Physico Chimie du Milieu Interstellaire'' and by the European Research Training Network ``Molecular Universe'' (MRTN-CT-2004-512302). R.M. has been supported by INTAS young research fellowship. The Kiel authors thank the Land Schleswig-Holstein for financial support. J.C. and B.T. acknowledge funding support from Spanish MEC through grants AYA2006-14876 and ESP2004-665, from PRICIT CM project S-0505/ESP-0237 (ASTROCAM).

References

 

  
Online Material

Table 3: Measured transitions of the ground vibrational state of CH3CH2C15N.

Table 4: Measured transitions of the ground vibrational state of CH3CHDCN.

Table 5: Measured transitions of the ground vibrational state of CH2DCH2CN in-plane.

Table 6: Measured transitions of the ground vibrational state of CH2DCH2CN out-of-plane.

Table 7: Sample table of the predicted transitions of the ground vibrational state of CH3CH2C15N.

Table 8: Sample table of the predicted transitions of the ground vibrational state of CH3CHDCN.

Table 9: Sample table of the predicted transitions of the ground vibrational state of CH2DCH2CN in-plane.

Table 10: Sample table of the predicted transitions of the ground vibrational state of CH2DCH2CN out-of-plane.

Table 11: Observed transitions of CH3CH2C15N (without high blend) in the frequency range of the Orion KL survey. Column 1 indicates the transition, Col. 2 provides the assumed rest frequency, Col. 3 the line strength, Col. 4 the energy of the upper level, Col. 5 the observed (centroid) frequencies assuming that the radial velocities relative to LSR are 5 km s-1, Col. 6 the observed peak line temperature of main beam, and Col. 7 indicates the main beam temperature derived from the model.



Copyright ESO 2009