Benjamin E.
Arenas
ab,
Sébastien
Gruet
abc,
Amanda L.
Steber
*abc,
Barbara M.
Giuliano
d and
Melanie
Schnell
*abc
aMax Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, D-22761 Hamburg, Germany. E-mail: amanda.steber@mpsd.mpg.de; melanie.schnell@mpsd.mpg.de
bCenter for Free-Electron Laser Science, Notkestrasse 85, D-22607 Hamburg, Germany
cThe Hamburg Centre for Ultrafast Imaging, Luruper Chaussee 149, D-22761 Hamburg, Germany
dMax Planck Institute for Extraterrestrial Physics, Giessenbachstrasse 1, D-85748 Garching, Germany
First published on 25th October 2016
We report here further spectroscopic investigation of the astrochemically relevant molecule i-propyl cyanide. We observed and analysed the rotational spectra of the ground state of the molecule and ten vibrationally excited states with energies between 180–500 cm−1. For this, we used a segmented W-band spectrometer (75–110 GHz) and performed the experiments under room temperature conditions. This approach thus provides access to high-resolution, pure rotational data of vibrational modes that occur in the far-infrared fingerprint region, and that can be difficult to access with other techniques. The obtained, extensive data set will support further astronomical searches and identifications, such as in warmer regions of the interstellar space where contributions from vibrationally excited states become increasingly relevant.
Recently, Belloche et al. detected the first branched alkyl molecule, i-propyl cyanide toward Sgr B2, with an abundance 0.4 times that of its straight-chain structural isomer n-propyl cyanide.8 Based on astrochemical models, the authors proposed that both isomers, i-propyl cyanide and n-propyl cyanide, are produced on dust grain ice mantles through the addition of molecular radicals. For the branched i-propyl cyanide species, the addition of a functional group to the non-terminal carbon would be required. This detection was based on an extensive laboratory spectroscopy investigation of i-propyl cyanide in selected regions between 6 and 600 GHz.9
A considerable proportion of complex molecules reside in low-lying vibrationally excited states at the elevated temperatures that were determined for i-propyl cyanide towards Sgr B2 (around 150 K). These vibrationally excited states, along with isotopologues that can often be seen in natural abundance, are associated with the so-called ‘weeds’ in radioastronomy data.10 These weeds can occupy many channels in a very complex, warm, dense astronomical dataset, therefore, it is important that the vibrationally excited states and isotopologues of abundant astronomical species are characterised. In 1974, Durig and Li reported the analysis of the ground state and three energetically low-lying vibrational states of i-propyl cyanide in a room-temperature microwave spectroscopy study (26.5–40 GHz).11
In the present work, we identified and characterised up to ten vibrationally excited states of i-propyl cyanide on our new segmented chirped-pulse W-band instrument, which is presented here. The energies of the vibrationally exited states ranged up to 500 cm−1, reinforcing how millimeter-wave spectroscopy of room temperature samples provides us with detailed molecular parameters of low-lying vibrationally excited states. We also observed the 13C isotopologues of the vibronic ground state of i-propyl cyanide in natural abundance (1.1%) in the millimeter-wave frequency range. This allowed us to precisely determine the ground-state structure of the molecule.
The spectrometer is combined with a stainless-steel single-pass vacuum chamber, which is 67 cm long and sealed with two Teflon windows at each side. The sample of interest is placed in a sample reservoir and can be heated to increase the vapour pressure. In this case, the i-propyl cyanide (99% pure, purchased from Alfa Aesar and used without further purification) did not require additional heating to generate sufficient vapour pressure. A flow of about 4 μbar of i-propyl cyanide was maintained inside the chamber during the experiment. The pulse train was broadcast with a horn antenna and coupled into the vacuum chamber through the first Teflon window. The molecular response, in the form of a free induction decay (FID), was coupled out of the chamber via the second Teflon window and detected by a receiving horn antenna. The pulse train and corresponding FIDs were collected on the real-time digitiser, averaged, and finally Fourier transformed via a fast Fourier transformation from the time domain to the frequency domain. A typical excitation pulse duration is 500 ns, and the FID is collected for 4 μs, as was the case for this experiment. The achieved frequency accuracy of the instrument is approximately 30 kHz with line widths of about 550 kHz. For this experiment, 500000 FIDs were recorded and averaged in HDR mode, which corresponds to a measurement time of 50 minutes.
Because of the fact that i-propyl cyanide is of Cs symmetry and the Ray's asymmetry parameter value is κ = −0.57, i-propyl cyanide is a near-prolate asymmetric molecule. However, to be consistent with previous datasets,9 the data were fitted to a Watson's S reduction Hamiltonian in a Ir representation. The spectral assignment was performed with the AABS package14,15 and fit with Pickett's programs (SPFIT/SPCAT).16 The experimental analysis was supported with harmonic and anharmonic quantum-chemical calculations using the program package Gaussian 09.17 These calculations, performed at the B3LYP/aug-cc-pVTZ level of theory, provided the rotational constants of the individual vibrationally excited states. The overall dipole moment of i-propyl cyanide is large, with μa = 4.02 D and μc = 0.62 D. For facilitating the assignment of the rich experimental spectrum, the calculated rotational constants for the vibrationally excited states were shifted by the difference we obtained between the calculated and the experimental rotational constants for the ground state (see Table S1 of the ESI† that provides a summary of the quantum-chemical results and the discrepancies between experiment and calculations).
For the rotational spectra of the non-perturbed vibrational states, we assigned between 66 to 128 transitions with a signal to noise ratio of 1200:
1 for the strongest transitions and fit them to a rigid-rotor Hamiltonian including centrifugal distortion within experimental accuracy (30 kHz). The transitions involved rotational quantum numbers J ranging between 11 and 17 and Ka values in the range 0 to 12 for R-type transitions, i.e., those following the selection rule ΔJ = +1. Ka is the quantum number associated with the projection of the rotational angular momentum onto the inertial axis a in the limit of the prolate symmetric top. Note that we also recorded weak Q-type transitions (obeying the selection rule ΔJ = 0) for which we observed higher J values. The rotational constants A, B and C, and the quartic centrifugal distortion rotational parameters of each state are well determined (see Table 1). The precision on the DK values is limited by the range of Ka values observable in this spectral region. For the perturbed vibrational states (those that lie close in energy), we only included transitions with low Ka values (up to 7) into the fit, which impacts the accuracy of the DK values. Therefore, in order to determine reliable rotational constants for those states, the DK constant has been kept fixed to its ground state value.
GS | ν 30 | ν 29* | ν 17* | ν 16 | ν 15 | 2ν30 | (ν30 + ν29)* | (ν30 + ν17)* | (ν29 + ν17)* | 2ν29* | 2ν17* | (ν30 + ν16) | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Symmetry | A′′ | A′′ | A′ | A′ | A′ | A′ | A′ | A′′ | A′′ | A′ | A′ | A′′ | |
ΔEa (cm−1) | 0 | 191.2 | 216.7 | 218.9 | 281.7 | 357.8 | 383.2 | 408.8 | 411.2 | 429.6 | 430.3 | 437.3 | 473.7 |
ΔEh (cm−1) | 0 | 189.7 | 220.9 | 218.3 | 281.9 | 353.2 | 379.3 | 410.6 | 408.0 | 439.3 | 441.8 | 436.7 | 471.6 |
A (MHz) | 7940.8746(16) | 7896.7463(15) | 7932.513(51) | 7993.882(26) | 7923.1764(24) | 7928.435(15) | 7853.1241(69) | 7879.481(28) | 7945.628(27) | 7926.817(25) | 7875.091(21) | ||
B (MHz) | 3968.08888(52) | 3980.03092(59) | 3964.9749(93) | 3971.8633(50) | 3970.27231(63) | 3978.7235(14) | 3991.94840(76) | 3976.7219(54) | 3983.3608(56) | 3963.9296(49) | 3981.9721(22) | ||
C (MHz) | 2901.05458(53) | 2899.86118(61) | 2898.2195(46) | 2904.7205(20) | 2902.82671(69) | 2900.5124(14) | 2898.60041(73) | 2896.9396(33) | 2903.4154(23) | 2896.3754(36) | 2901.5706(16) | ||
D J (kHz) | 0.6127(12) | 0.6386(13) | 0.630(16) | 0.6085(79) | 0.6129(14) | 0.6155(16) | 0.6603(15) | 0.6349(75) | 0.6205(82) | 0.5585(60) | 0.6320(34) | ||
D JK (kHz) | 12.1698(43) | 11.4510(21) | 11.38(24) | 13.162(99) | 11.9422(17) | 12.2310(37) | 10.7860(38) | 10.363(50) | 12.519(92) | 11.231(95) | 11.279(13) | ||
D K (kHz) | −5.231(27) | −7.869(10) | [−5.231] | [−5.231] | −5.236(10) | −4.07(11) | −9.356(85) | [−5.231] | [−5.231] | [−5.231] | −7.20(55) | ||
d 1 (kHz) | −0.24393(11) | −0.25602(18) | −0.2637(72) | −0.2486(39) | −0.25187(18) | −0.25600(83) | −0.27202(36) | −0.3044(60) | −0.2550(48) | −0.2516(46) | −0.2556(22) | ||
d 2 (kHz) | −0.189404(63) | −0.183884(73) | −0.1774(31) | −0.2025(13) | −0.185515(88) | −0.19630(55) | −0.17902(17) | −0.1983(24) | −0.1934(15) | −0.1929(25) | −0.1760(10) | ||
H JK (Hz) | 0.0347(16) | −2.55(77) | −1.55(31) | −1.59(26) | |||||||||
h 2 (Hz) | 0.000931(72) | ||||||||||||
h 3 (Hz) | 0.000326(33) | ||||||||||||
No. lines | 173 | 128 | 45 | 71 | 109 | 88 | 89 | 42 | 68 | 37 | 66 | ||
rms (kHz) | 26.9 | 32.3 | 58.5 | 44.1 | 33.1 | 33.7 | 32.0 | 64.4 | 50.3 | 51.9 | 32.7 | ||
J m/Ka,m | 40/10 | 35/12 | 16/6 | 17/8 | 40/14 | 40/13 | 38/12 | 17/8 | 17/8 | 17/6 | 17/11 |
The additional rotational constants due to 13C isotopic substitution in natural abundance were used to determine the substitution structure of i-propyl cyanide using Kraitchman's equations18 as implemented in the KRA program package.14
![]() | ||
Fig. 1 (a) Broadband rotational spectrum of i-propyl cyanide and (b) zoom of the ![]() ![]() ![]() ![]() ![]() |
Based on the shifted calculated rotational constants (see Section 2 and Table S1 of the ESI†) the assignment of most of the vibrationally excited states was straight forward. We could unambiguously identify ten energetically low-lying vibrationally excited states and fit their rotational transitions to an asymmetric-top Hamiltonian, as described in Section 2. The rotational constants for the ground state and the vibrationally excited states are summarised in Table 1. The calculated vibrational frequencies of the respective states are also included (B3LYP/aug-cc-pVTZ level of theory with and without anharmonic corrections). However, for some pairs of states, such as ν29/ν17 as well as (ν30 + ν29)/(ν30 + ν17), the analysis was more challenging. We found that these vibrationally excited states, which are within only a few cm−1 energetically, perturb each other. Without being able to identify the band centers, this perturbation cannot be analysed in detail. Such a study would require high-resolution far-infrared data, where vibrational transitions are probed, and our anharmonic frequency calculations predict the relevant vibrational transitions to be too weak to be observed. Even though the standard deviation for these states of the corresponding fits are larger than the frequency accuracy of our spectrometer, they are still sufficiently determined to obtain useful molecular parameters to support their radioastronomic detection, which is one of the motivations of the present work.
The results of our quantum-chemical calculations allow us to assign motions for some of the vibrational states. The lowest frequency vibrational state, ν30, corresponds to the C–CN out-of-plane bending motion, and ν17 can be assigned to the associated in-plane motion. The ν16 and ν29 modes concern the symmetric and asymmetric methyl torsional modes, respectively. Finally, the last fundamental mode assigned, ν15, refers to a CC2 deformation involving a bending motion between the two C–CH3 bonds. All the other modes observed in this study are combinations of the previously described motions. The (ν29 + ν17) and 2ν17 modes could not be observed experimentally. We assume that they are too strongly perturbed because they are involved (with 2ν29) in a three party perturbation.
Note that three of the vibrationally excited states (namely ν29, ν17, ν16) were also observed in the previous study by Durig et al. (Table 3 of ref. 11). By comparing the respective rotational constants, we can make the following correlation: the να, νβ and νγ states from ref. 11 correspond to our ν29, ν17, and ν16 states, respectively. Based on transition intensities, the authors of ref. 11 ordered their states according to νβ, να, and νγ, which generally agrees with our findings. However, they did not report on the lowest vibrational state, i.e., ν30 in Table 1. This missing state ν30 thus led to some confusion in the initial assignments.
The authors of ref. 11 also derived an approximate molecular structure for i-propyl cyanide. The length of the C–CN bond is particularly interesting since it is expected to be rather short as a result of the very strong electron-withdrawing cyano group. In ref. 11, the bond length r(C–CN) was determined to be 1.501 Å and the angle between the CN bond and the CCC plane was estimated to be σ = 53.8°. In the present work, we recorded the rotational spectra of all three singly substituted 13C isotopologues of the vibronic ground state of i-propyl cyanide. These additional sets of rotational constants allowed us to determine its substitution structure (rs) of the carbon backbone using the Kraitchman's equations (see also Table S2 of the ESI†). This method exploits mono-isotopic substitution information and provides a straightforward way to use the isotopic changes in the moments of inertia to determine the atom coordinates in the principal axis frame.
Fig. 2 shows a comparison between the calculated structure of i-propyl cyanide (B3LYP/aug-cc-pVTZ) and the experimentally determined rs structure. As can be seen, the experimental atom positions match those from the calculated structure with small deviations. The substitution structure yields the following structural parameters: r(C–CN) = 1.4600(46) Å, which is 0.04 Å shorter than determined in ref. 11, and r(C–CH3) = 1.5359(42) Å. The angle between the CN bond and the CCC plane was determined to be σ = 51.52(32)°, which is also smaller than in ref. 11. This is in line with the qualitative discussions in ref. 11 that a shorter C–CN bond results in smaller σ angles. The discrepancy between our structural parameters and those from Durig and Li might arise from the various assumptions they had to adopt for their analysis. Their data only allowed them to produce a partial structure, and they presented parameters that were determined with the aid of i-propyl chloride.11 In the present work, the complete carbon backbone structure is determined solely based on the experimental data.
The comparison of the C–CN bond length obtained for i-propyl cyanide can be extended to other cyanide containing molecules. Keeping in mind that the electron-withdrawing cyano group is expected to produce a shortening of that C–CN single bond, Table 2 compares the C–CN bond length in i-propyl cyanide with analogous, short-chain cyanides as well as with average C–C bond lengths for different orbital hybridisation regimes. The C–CN bond in i-propyl cyanide is on the same order of magnitude as an average C(sp3)–C(sp) bond, being only slightly shorter by 0.01 Å. Indeed, the bond is substantially shorter than a typical C–C single bond (C(sp3)–C(sp3)) by 0.07 Å. In addition, the C–CN bond length in i-propyl cyanide is in keeping with those of analogous alkyl cyanides – it amounts to 1.46 Å (when errors are considered), regardless of the adjacent alkyl chains. The shortest C–CN bond is the one in benzonitrile, likely due to electronic effects of the aromatic ring giving it some double bond character. The experimental data are compared to computed bond distances that we performed at the B3LYP/aug-cc-pVTZ level of theory. The similar bond lengths in the alkyl cyanides and the respective shorter value for benzonitrile is well reproduced by the calculations.
Compound | r(C–CN)exp (Å) | r(C–CN)calc (Å) |
---|---|---|
Methyl cyanide | 1.4594(6)19 | 1.455 |
Ethyl cyanide | 1.459(1)20 | 1.460 |
n-Propyl cyanide | 1.4621(43)21 | 1.459 |
i-Propyl cyanide | 1.50111 | 1.466 |
i-Propyl cyanide | 1.4600(46) (this work) | 1.466 |
n-Butyl cyanide | 1.459 | |
t-Butyl cyanide | 1.46(2)22 | 1.472 |
Benzonitrile | 1.4509(6)23 | 1.430 |
C(sp3)–C(sp3) | 1.5324 | |
C(sp3)–C(sp) | 1.4724 |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6cp06297k |
This journal is © the Owner Societies 2017 |