# Spectroscopic and Computational Characterization of 2-Aza-1,3-butadiene, a Molecule of Astrochemical Significance

Ningjing Jiang, Mattia Melosso, Luca Bizzocchi, Silvia Alessandrini,
Jean-Claude Guillemin, Luca Dore, Cristina Puzzarini

<figure><img src="image_1.png" /></figure>

## To cite this version:

Ningjing Jiang, Mattia Melosso, Luca Bizzocchi, Silvia Alessandrini, Jean-Claude Guillemin, et al. Spectroscopic and Computational Characterization of 2-Aza-1,3-butadiene, a Molecule of Astrochemical Significance. Journal of Physical Chemistry A, 2022, 126 (11), pp.1881-1888. 10.1021/acs.jpca.2c00831. hal-03632037v2

HAL Id: hal-03632037

https://hal.science/hal-03632037v2

Submitted on 6 Apr 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

# Spectroscopic and Computational Characterization of 2-Aza-1,3-butadiene, a Molecule of Astrochemical Significance

*Published as part of The Journal of Physical Chemistry virtual special issue “10 Years of the ACS PHYS Astrochemistry Subdivision”.*

Ningjing Jiang, Mattia Melosso,* Luca Bizzocchi, Silvia Alessandrini, Jean-Claude Guillemin, Luca Dore, and Cristina Puzzarini*

<figure><img src="image_7.png" /></figure>

Cite This: *J. Phys. Chem. A* 2022, 126, 1881–1888

<figure><img src="image_8.png" /></figure>

Read Online

## ACCESS

<figure><img src="image_9.png" /></figure>

Metrics & More

<figure><img src="image_10.png" /></figure>

Article Recommendations

<figure><img src="image_11.png" /></figure>

Supporting Information

**ABSTRACT:** Being N-substituted unsaturated species, azabutadienes are molecules of potential relevance in astrochemistry, ranging from the interstellar medium to Titan's atmosphere. 2-Azabutadiene and butadiene share a similar conjugated π system, thus allowing investigation of the effects of heteroatom substitution. More interestingly, 2-azabutadiene can be used to proxy the abundance of interstellar butadiene. To enable future astronomical searches, the rotational spectrum of 2-azabutadiene has been investigated up to 330 GHz. The experimental work has been supported and guided by accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the two possible forms, *trans* and *gauche*. The *trans* species, more stable by about 7 kJ/mol than *gauche*-2-azabutadiene, has been experimentally observed, and its rotational and centrifugal distortion constants have been obtained with remarkable accuracy, while theoretical estimates of the spectroscopic parameters are reported for *gauche*-2-azabutadiene.

<figure><img src="image_12.png" /></figure>

## INTRODUCTION

1,3-Butadiene ($CH_2=CH-CH=CH_2$, hereafter denoted as butadiene) is a textbook species for understanding conjugation effects and, being the simplest conjugated diene, represents the prototypical reactant for Diels-Alder reactions. Molecules that are isoelectronic to butadiene and share a similar conjugated π system form an interesting family of species. Indeed, they allow the investigation of the heteroatom effects in the butadiene skeleton.¹ Limiting ourselves to the second row of the periodic table, there are three singly substituted species isoelectronic to butadiene: acrolein ($O=CH-CH=CH=CH_2$, also known as propenal), 1-aza-1,3-butadiene ($NH=CH-CH=CH_2$, hereafter 1-azabutadiene), and 2-aza-1,3-butadiene ($CH_2=N-CH=CH_2$, hereafter 2-azabutadiene). In this context, high-resolution molecular spectroscopy and computational chemistry are among the most powerful tools to elucidate the molecular structure, internal dynamics, and conformational behavior of these species, especially when the two approaches are combined synergistically.²

These butadiene-like species and butadiene itself are not easy to characterize experimentally. For example, due to the lack of a permanent electric dipole moment, the most stable *trans* form of butadiene does not possess any rotational spectrum. Hence, it is not possible to take advantage of the

inherent accuracy of rotational spectroscopy to determine its molecular structure. In such a case, a rigorous geometrical determination could only be achieved by the analysis of ro-vibrational spectra for several isotopic species.³⁻⁷ As far as *gauche*-butadiene is concerned, its higher energy with respect to the *trans* form (~12 kJ/mol) together with its small dipole moment value (0.09 D) have prevented its nonplanarity from being established unambiguously until recently.⁸

Whenever a heteroatom is introduced into the butadiene skeleton, the molecular symmetry is noticeably reduced, the resulting asymmetry then generating both an appreciable electric dipole moment and an increased number of active infrared modes. For these reasons, acrolein, 1-azabutadiene, and 2-azabutadiene represent more suitable spectroscopic targets. For example, owing to its large dipole moment values (about 3 and 2.5 D for the *trans* and *cis* forms, respectively⁹), acrolein has been well characterized from a spectroscopic point

Received: February 3, 2022

Revised: March 5, 2022

Published: March 11, 2022

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>
of view, as demonstrated by the vast literature on this subject.<sup>10–16</sup> On the other hand, spectroscopic study of 1- and 2-azabutadiene has been quite limited to date<sup>17–21</sup> as both species are unstable and require an efficient and selective production method.

Exhaustive characterization of the high-resolution rotational spectra for these butadiene-like species is also crucial for a different aspect: the study of the chemistry occurring in the interstellar medium (ISM). Despite the harsh conditions of molecular clouds, i.e., very low temperatures and number densities, unsaturated carbon chains are often abundant (see, for example, the case of TMC-1<sup>22</sup> or Lupus-1A<sup>23</sup>), thus making butadiene and heterodienes good candidates for astronomical searches. However, because of its centrosymmetric nature, trans-butadiene cannot be detected in space through rotational emission and observation of the more energetic gauche form in cold molecular clouds is very unlikely. The abundance of centrosymmetric species, such as molecular nitrogen, cyanogen, or benzene, has been indirectly estimated through observation of their protonated,<sup>24</sup> isomeric,<sup>25</sup> or functionalized forms.<sup>26</sup> Analogously, we suggest that the presence of interstellar butadiene can be proxied by the detection of its isoelectronic heterodienes. This strategy represents a promising route given the fact that acrolein has been recently discovered in the ISM<sup>27</sup> and the detection of 1-azabutadiene is being reported tentatively.<sup>28</sup> However, to the best of our knowledge, astronomical searches of 2-azabutadiene have never been reported. This is probably due to the fact that in addition to the limited spectral data reported in the literature, no entry for 2-azabutadiene is available in the widely used CDMS database<sup>29,30</sup> or in the JPL catalog.<sup>31</sup>

In this respect, here, we report a refined analysis of the rotational spectrum of 2-azabutadiene based on new rotational transitions recorded between 225 and 330 GHz. The experimental work has been supported by high-level quantum-chemical calculations, which guided us through the spectral interpretation and provided a good reference set of spectroscopic parameters. This manuscript is organized as follows. In the next two sections, quantum-chemical computations and experimental details are discussed. Then, the results on 2-azabutadiene are presented in the fourth section, while a comparison with other dienes is discussed in the last section.

## ■ COMPUTATIONAL DETAILS

The leading terms in rotational spectroscopy are the rotational constants, which can be computationally predicted by summing two contributions: the equilibrium rotational constants ($B_e$), straightforwardly derived from the equilibrium structure, and the vibrational corrections ($\Delta B_{vib}$)

$$B_0^i = B_e^i + \Delta B_{vib}^i = B_e^i - \frac{1}{2} \sum_r \alpha_r^i \quad (1)$$

In the equation above, $i$ denotes the inertial axis ($i = a, b, c$, that is, for instance, $B_0^i = A$), and the sum runs over all vibrational normal modes. Within vibrational perturbation theory to second order (VPT2), the vibrational corrections require computation of the vibration-rotation interaction constants ($\alpha_r^i$).<sup>32</sup> Their evaluation requires an anharmonic force field,<sup>33</sup> which has been calculated at the MP2/cc-pVTZ<sup>34</sup> level, within the frozen-core (fc) approximation. The MP2

acronym stands for Møller-Plesset perturbation theory to the second order.<sup>35</sup>

Despite being important for quantitative predictions,<sup>36,37</sup> the vibrational corrections provide a small contribution to the $B_0$ constants, indeed accounting only for ~1–3% of the total value. This implies that the equilibrium structure should be evaluated with high accuracy. For this purpose, the so-called CBS+CV composite scheme has been employed,<sup>38</sup> which is based on the CCSD(T) method (coupled-cluster (CC) singles and doubles with perturbative treatment of triples<sup>39</sup>) and requires minimization of the following energy gradient

$$\frac{dE_{\text{CBS+CV}}}{dx} = \frac{dE_{\text{HF-SCF}}^{\infty}}{dx} + \frac{d\Delta E_{\text{CCSD(T)}}^{\infty}}{dx} + \frac{d\Delta E_{\text{CV}}}{dx} \quad (2)$$

where $dE_{\text{HF-SCF}}/dx$ and $d\Delta E_{\text{CCSD(T)}}/dx$ are gradients for the extrapolation to the complete basis set (CBS) limit of the HF-SCF energy (exponential extrapolation formula by Feller<sup>40</sup>) and of the CCSD(T) electron correlation contribution ($n^{-3}$ extrapolation expression<sup>41</sup>), respectively. Since the extrapolation to the CBS limit is performed within the frozen core approximation, core-valence (CV) correlation effects are incorporated by adding the corresponding correction, $d\Delta E_{CV}/dx$. This involves the difference of all-electron (ae) and frozen-core CCSD(T) calculations using the same basis set. In the CBS+CV scheme, for the extrapolation to the CBS limit, we used the cc-pVnZ basis sets,<sup>34</sup> with $n = T - 5$ for HF-SCF and $n = T$, Q for CCSD(T). For the CV correction, we resorted to the cc-pCVTZ basis set.<sup>42</sup>

The semirigidity of the molecules also requires one to account for the effect of centrifugal distortion. Different orders of centrifugal distortion terms are possible according to the power of the angular momentum operator considered in the effective rotational Hamiltonian. However, from a computational point of view, only the quartic and sextic centrifugal-distortion constants can be obtained: the former as a byproduct of harmonic force-field computations, and the latter from the anharmonic force field. Therefore, the sextic terms have been computed at the MP2/cc-pVTZ level, while the quartic centrifugal-distortion constants have been purposely derived from a fc-CCSD(T)/cc-pVTZ harmonic force field. This latter calculation also allowed for determining the electric dipole moment ($\mu$) components as well as the electric field gradients at the nitrogen nucleus, which straightforwardly provide the nitrogen quadrupole coupling constants.<sup>33</sup>

All quantum-chemical computations have been carried out with the CFOUR package.<sup>43,44</sup>

## ■ EXPERIMENTAL SECTION

The sample of 2-azabutadiene was synthesized in a two-step procedure, and its rotational spectrum was recorded with a frequency-modulation spectrometer.

The synthetic reaction reported in the literature for obtaining 2-azabutadiene<sup>45,46</sup> has been slightly revised, and the experimental conditions were optimized with the aim of maximizing the signal-to-noise ratio (S/N) of the rotational spectrum and minimizing the number of possible interfering lines. The production of 2-azabutadiene proceeds in two steps, detailed in the following: synthesis of the precursor (step I) and its transformation into the title molecule (step II).

**Step I: Synthesis of the Precursor.** The precursor is 2-[(2-chloroethyl)amino]acetonitrile, which has been synthesized in a Strecker reaction starting from 2-chloroethylamine
<figure><img src="image_15.png" /><figcaption>**Figure 1.** Portion of the rotational spectra around the $J_{K\alpha,K\epsilon} = 27_{0,27} \leftarrow 26_{1,26}$ transition of *trans*-2-azabutadiene recorded while optimizing the experimental conditions. Two lines adjacent to the one of interest belong to *t*-BuOH, formed during the reaction. (A) Several spectra acquired during optimization of the reactant temperatures, whereas (B) spectra recorded keeping 2-[(2-chloroethyl)amino]acetonitrile at 65 °C and *t*-BuOK at room temperature while varying the temperature of the cooling bath.</figcaption></figure>

hydrochloride (0.1 mol, 11.6 g) and sodium cyanide (0.1 mol, 5 g) in methanol (30 mL) and water (20 mL). The reaction mixture was cooled to 0 °C, and formaldehyde (37% in water, 0.1 mol) was added slowly. The solution was allowed to warm up to room temperature and stirred for 1 h. Organic compounds were then extracted with dichloromethane (3 × 50 mL). The organic phase was dried over magnesium sulfate, and the solvents were removed in vacuo. The crude product obtained in a 90% yield (10.6 g) was sufficiently pure to be used directly in the next step.⁴⁵,⁴⁶ Distillation under vacuum (0.1 mbar) gave the pure product in a 43% yield (5.1 g).

2-Chloroethylamine hydrochloride, formaldehyde in water, and sodium cyanide were purchased from Sigma-Aldrich and used without further purification.

**Step II: Synthesis of 2-Azabutadiene.** 2-Azabutadiene was synthesized by vaporizing 2-[(2-chloroethyl)amino]-acetonitrile (2 g) on potassium *tert*-butylate (*t*-BuOK, 30 g) introduced in a horizontal tube and filling its lower half. In this reaction, *tert*-butanol was also produced and selectively trapped in a U-tube immersed in a cooling bath, while a gaseous flow of almost pure 2-azabutadiene was directly introduced into the spectrometer cell.

The best conditions were optimized by monitoring the spectral S/N for the $J_{K\alpha, K\epsilon} = 27_{0,27} \leftarrow 26_{1,26}$ transition of *trans*-2-azabutadiene, predicted on the basis of literature data.²⁰ Figure 1 shows several spectra recorded upon changing the temperature of (i) the precursor species, (ii) *t*-BuOK, and (iii) the cooling bath. Ultimately, the maximum yield of 2-azabutadiene (and consequently the best S/N in the spectrum) was achieved by maintaining the precursor compound between 65 and 90 °C, the potassium *tert*-butylate sample at room temperature, and the cooling bath at -80 °C. The cooling bath was obtained by mixing acetone with liquid nitrogen and allowed for complete removal of the interfering lines belonging to the *t*-BuOH formed during the reaction.

**Frequency-Modulation Millimeter Spectrometer.** The rotational spectrum of 2-azabutadiene has been recorded between 225 and 330 GHz by means of a frequency-modulation millimeter-wave spectrometer in Bologna; a

detailed description of the instrument is given in previous works.⁴⁷,⁴⁸ Briefly, the primary radiation source—whose frequency and phase stability are ensured by a phase-lock loop (PLL)—is a Gunn diode emitting in the W band (75–110 GHz). The Gunn diode is coupled to a passive multiplier (tripler, WR3.4x3 Virginia Diodes Inc.) to produce higher frequencies. The modulated radiation source is fed into a 3 m glass absorption cell where the vapors of 2-azabutadiene are injected at a pressure between 5 and 10 µbar. The output signal is detected by a Schottky barrier diode and demodulated by a lock-in amplifier set at twice the modulation frequency (2f). The uncertainties of the present experiment are around 30 kHz.

■ RESULTS

2-Azabutadiene exists in two conformations, depending on the C=N—C=C dihedral angle ($\phi$): the most stable planar *trans* form ($\phi = 180^{\circ}$) and the more energetic *gauche* isomer ($\phi = 54^{\circ}$). The optimized equilibrium structures are shown in Figure 2 together with the corresponding CBS+CV geometrical parameters (Cartesian coordinates are available in the Supporting Information). Except for the dihedral angle, most of the bond lengths and angles are similar in the two forms, i.e., they agree within few milliAngstroms for the former and ~1° for the latter. The largest differences are observed for the H–C=N angle (3.7°) and the two involved bonds (7 mÅ for C= N and 9 mÅ for C–H). Understanding the reasons at the basis of these differences is beyond the scope of our paper; however, they can be simply rationalized by performing natural bond orbital (NBO) analysis. The results of such analysis are collected in the SI and can be used to understand how the conjugation effects stabilize the *trans* species more than the *gauche* form.

At the CBS+CV level, the equilibrium energy difference ($\Delta E[\text{gauche} - \text{trans}]$) between the two conformers is 6.5 kJ/mol. By incorporating the harmonic zero-point energy (ZPE) correction, computed at the fc-CCSD(T)/cc-pVTZ level of theory, $\Delta E[\text{gauche} - \text{trans}]$ increases to 7.2 kJ/mol. On the basis of this energy difference, the population of *gauche*-2-
<figure><img src="image_16.png" /><figcaption>**Figure 2.** Molecular structures of *trans*- and *gauche*-2-azabutadiene together with the optimized geometrical parameters (CBS+CV level, see text; bond lengths in Angstroms; angles in degrees).</figcaption></figure>

azabutadiene is about 7% of that of the *trans* form at room temperature. Since the computed ($\mu_b^{(gauche)}/\mu_b^{(trans)}$)² ratio is ~0.7 (see Table 1), rotational transitions belonging to *gauche*-2-azabutadiene are expected to be roughly 20 times weaker than those of the *trans* species. Even though the attained S/N of the spectrum is about 150 for *trans*-2-azabutadiene, this is expected to be less than 10 for the *gauche* form. In addition, the *trans* species possesses eight vibrational states lying below

the energy of *gauche*-2-azabutadiene. Therefore, in conclusion, identification of the weak spectrum of the *gauche* form is prevented by the huge number of interfering lines, and hereafter, we will focus on the only conformer securely characterized in this work, i.e., the *trans* species. Nonetheless, the theoretical spectroscopic constants provided in this manuscript (see Table 1) represent a solid base for attempting the identification of *gauche*-2-azabutadiene in different experiments (e.g., FTMW spectroscopy of a jet-cooled sample).

Moving to the fully characterized *trans*-2-azabutadiene species, this is an asymmetric rotor with a strong nearly prolate nature ($\kappa = -0.979$). In its principal inertia system, the molecule lies on the *ab* plane with the *c* axis perpendicular to it. The computed dipole moment components are $\mu_a = 0.47$ D and $\mu_b = 1.46$ D (see Table 1), with $\mu_c$ being null for symmetry reasons (the point group symmetry is $C_s$).

Exploiting the best experimental conditions, almost 600 new rotational transitions, both of *a* and of *b* type, were recorded between 225 and 330 GHz (examples are provided in Figure 3). They probe energy levels with *J* up to 76 and $K_a$ up to 20, such a difference in the upper limits of *J* and $K_a$ being ascribable to the relatively large *A* rotational constant (~47 GHz) with respect to the much smaller *B* and *C* values (~4.4 GHz). Indeed, the population of the rotational energy levels decreases quickly by increasing $K_a$ and slowly by increasing *J*. Our newly observed rotational transitions were combined with previous low-frequency data²⁰ and analyzed in a weighted

<table><caption>**Table 1.** Ground-State Spectroscopic Parameters of *trans*- and *gauche*-2-Azabutadiene (S Reduction, *I* Representation)</caption><thead><tr><th rowspan="2">parameters</th><th rowspan="2">units</th><th colspan="3">trans</th><th colspan="2">gauche</th></tr><tr><th>exp.<sup>a</sup></th><th>theo.<sup>b</sup></th><th>ref 20<sup>c</sup></th><th>theo.<sup>b</sup></th><th>scaled<sup>d</sup></th></tr></thead><tbody><tr><td>A</td><td>MHz</td><td>47 186.0004(7)</td><td>47231.509</td><td>47 186.012(7)</td><td>24 021.128</td><td>23 997.983</td></tr><tr><td>B</td><td>MHz</td><td>4886.5219(1)</td><td>4887.740</td><td>4886.5223(8)</td><td>6071.991</td><td>6070.478</td></tr><tr><td>C</td><td>MHz</td><td>4430.0772(1)</td><td>4430.824</td><td>4430.0778(8)</td><td>5120.829</td><td>5119.966</td></tr><tr><td>D<sub>J</sub></td><td>kHz</td><td>0.97917(8)</td><td>0.960416</td><td>0.980(3)</td><td>5.42141</td><td>5.52727</td></tr><tr><td>D<sub>JK</sub></td><td>kHz</td><td>-8.0641(7)</td><td>-7.38104</td><td>-7.99(9)</td><td>-40.8174</td><td>-44.5947</td></tr><tr><td>D<sub>K</sub></td><td>kHz</td><td>341.95(2)</td><td>309.12</td><td>341.(1)</td><td>207.669</td><td>229.722</td></tr><tr><td>d<sub>1</sub></td><td>kHz</td><td>-0.117415(5)</td><td>-0.120169</td><td>-0.1173(1)</td><td>-1.82008</td><td>-1.77837</td></tr><tr><td>d<sub>2</sub></td><td>kHz</td><td>-0.006949(1)</td><td>-0.006235</td><td>-0.00695(2)</td><td>-0.19994</td><td>-0.22283</td></tr><tr><td>H<sub>J</sub></td><td>mHz</td><td>0.29(1)</td><td>0.217</td><td></td><td>-12.7</td><td>-16.97</td></tr><tr><td>H<sub>JK</sub></td><td>mHz</td><td>-8.0(2)</td><td>-3.856</td><td></td><td>308.2</td><td>639.4</td></tr><tr><td>H<sub>KJ</sub></td><td>Hz</td><td>-0.267(2)</td><td>-0.411</td><td></td><td>-3.736</td><td>-2.431</td></tr><tr><td>H<sub>K</sub></td><td>Hz</td><td>1.3(1)</td><td>5.75</td><td></td><td>14.019</td><td>3.167</td></tr><tr><td>h<sub>1</sub></td><td>mHz</td><td>0.0691(9)</td><td>0.0786</td><td></td><td>-7.35</td><td>-6.462</td></tr><tr><td>h<sub>2</sub></td><td>mHz</td><td>0.0104</td><td>0.0104</td><td></td><td>-1.12</td><td>-1.12</td></tr><tr><td>h<sub>3</sub></td><td>mHz</td><td>0.0016</td><td>0.0016</td><td></td><td>-0.0691</td><td>-0.0691</td></tr><tr><td>1.5 &times; &chi;<sub>aa</sub></td><td>MHz</td><td>1.411</td><td>1.411</td><td>0.8(15)</td><td>1.275</td><td>1.275</td></tr><tr><td>(&chi;<sub>bb</sub> &minus; &chi;<sub>cc</sub>)/4</td><td>MHz</td><td>-2.142(7)</td><td>-2.13</td><td>-2.1(12)</td><td>-0.839</td><td>-0.839</td></tr><tr><td>|&mu;<sub>a</sub>|</td><td>D</td><td>0.44(3)<sup>e</sup></td><td>0.47</td><td></td><td>0.029</td><td>0.027</td></tr><tr><td>|&mu;<sub>b</sub>|</td><td>D</td><td>1.90(7)<sup>e</sup></td><td>1.46</td><td></td><td>1.206</td><td>1.569</td></tr><tr><td>|&mu;<sub>c</sub>|</td><td>D</td><td></td><td></td><td></td><td>0.806</td><td>0.806</td></tr><tr><td>no. of data</td><td></td><td>675</td><td></td><td>82</td><td></td><td></td></tr><tr><td>J<sub>max</sub> K<sub>a max</sub></td><td></td><td>76, 20</td><td></td><td>31, 4</td><td></td><td></td></tr><tr><td>rms error</td><td>kHz</td><td>32.9</td><td></td><td>31.8</td><td></td><td></td></tr><tr><td>&sigma;</td><td></td><td>1.09</td><td></td><td>1.06</td><td></td><td></td></tr></tbody></table>

<sup>a</sup>Values in parentheses are one standard deviation and refer to the last digits. Parameters without uncertainties are held fixed at the corresponding computed value. <sup>b</sup>CBS+CV equilibrium rotational constants augmented by vibrational corrections at the fc-MP2/cc-pVTZ level. Quartic centrifugal distortion and nuclear quadrupole coupling constants computed at the fc-CCSD(T)/cc-pVTZ level. Sextic centrifugal distortion terms at the fc-MP2/cc-pVTZ level. Equilibrium dipole moment components calculated at the fc-CCSD(T)/cc-pVTZ level and augmented by vibrational corrections at the fc-MP2/cc-pVTZ level. <sup>c</sup>Converted from A to S reduction. <sup>d</sup>Obtained using eq 3. See text for details. <sup>e</sup>Determined by Stark effect measurements in ref 20.
<figure><img src="image_17.png" /><figcaption>**Figure 3.** Portion of the rotational spectrum of 2-azabutadiene (black). Red trace depicts *a*-type transitions (A), while green trace represents *b*-type transitions (B).</figcaption></figure>

<table><caption>**Table 2.** Rotational and Quartic Centrifugal Distortion Constants of the *Trans* Species<sup>*a*</sup></caption><thead><tr><th rowspan="2">constant</th><th rowspan="2">units</th><th colspan="2">butadiene<sup>b</sup></th><th colspan="2">1-azabutadiene<sup>c</sup></th><th colspan="2">acrolein<sup>d</sup></th><th colspan="2">2-azabutadiene</th></tr><tr><th>anti</th><th>syn</th><th>anti</th><th>syn</th><th>acrolein</th><th>2-azabutadiene</th><th>anti</th><th>syn</th></tr></thead><tbody><tr><td>A</td><td>MHz</td><td>41 682.6292(11)</td><td>45 773.6451(9)</td><td>43 755.656(2)</td><td>47 353.7034(20)</td><td>47 186.0004(7)</td><td>41 682.6292(11)</td><td>45 773.6451(9)</td><td>43 755.656(2)</td></tr><tr><td>B</td><td>MHz</td><td>4433.50067(18)</td><td>4560.93131(7)</td><td>4564.54014(7)</td><td>4659.48810(18)</td><td>4886.5219(1)</td><td>4433.50067(18)</td><td>4560.93131(7)</td><td>4564.54014(7)</td></tr><tr><td>C</td><td>MHz</td><td>4008.04876(19)</td><td>4148.24852(6)</td><td>4134.46273(6)</td><td>4242.70127(16)</td><td>4430.0772(1)</td><td>4008.04876(19)</td><td>4148.24852(6)</td><td>4134.46273(6)</td></tr><tr><td>D<sub>J</sub></td><td>kHz</td><td>0.861735(27)</td><td>0.90964(3)</td><td>0.95037(4)</td><td>1.028148(56)</td><td>0.97917(8)</td><td>0.861735(27)</td><td>0.90964(3)</td><td>0.95037(4)</td></tr><tr><td>D<sub>JK</sub></td><td>kHz</td><td>-7.15724(48)</td><td>-7.531(1)</td><td>-7.589(1)</td><td>-8.70159(74)</td><td>-8.0641(7)</td><td>-7.15724(48)</td><td>-7.531(1)</td><td>-7.589(1)</td></tr><tr><td>D<sub>K</sub></td><td>kHz</td><td>218.6356(33)</td><td>296.80(3)</td><td>257.29(7)</td><td>360.142(30)</td><td>341.95(2)</td><td>218.6356(33)</td><td>296.80(3)</td><td>257.29(7)</td></tr><tr><td>d<sub>1</sub></td><td>kHz</td><td>-0.107889(16)</td><td>-0.10715(1)</td><td>-0.114863(5)</td><td>-0.120228(21)</td><td>-0.117415(5)</td><td>-0.107889(16)</td><td>-0.10715(1)</td><td>-0.114863(5)</td></tr><tr><td>d<sub>2</sub></td><td>kHz</td><td>-0.0056450(72)</td><td>-0.005863(2)</td><td>-0.006237(3)</td><td>-0.0069896(77)</td><td>-0.006949(1)</td><td>-0.0056450(72)</td><td>-0.005863(2)</td><td>-0.006237(3)</td></tr></tbody></table>

<sup>*a*</sup>Values in parentheses are one standard deviation and refer to the last digits. <sup>*b*</sup>From ref 3. <sup>*c*</sup>From ref 28. <sup>*d*</sup>From ref 16. Converted from A to S reduction.

least-squares procedure using the SPFIT program.<sup>49</sup> The optimized energy level positions were calculated employing a Watson S-reduced Hamiltonian with centrifugal distortion terms up to the sixth power of the angular momentum operator.<sup>50</sup> Despite the relatively high *J* values targeted in this work, the hyperfine structure arising from the quadrupolar <sup>14</sup>N nucleus has been partially resolved for some low *K*<sub>a</sub> *b*-type transitions. Therefore, we were able to determine the (*χ*<sub>bb</sub> – *χ*<sub>cc</sub>)/4 parameter, whereas the value of 1.5 × *χ*<sub>aa</sub> has been kept fixed to its computed value. The derived spectroscopic parameters are collected in Table 1, where theoretical estimates are also provided.

From inspection of Table 1, we note that the root-mean-square error of the residuals (32.9 kHz) is in line with the estimated accuracy of both our transition frequencies and those reported in the literature.<sup>20</sup> The semirigid model is appropriate for the description of the rotational energy levels, as demonstrated by the standard deviation of the fit (*σ*), which is 1.09. The rotational constants *A*, *B*, and *C* have been determined with remarkably good accuracy, the associated errors being smaller than 1 kHz. Notably, all quartic centrifugal distortion constants have been obtained with an uncertainty smaller than 0.01%, and most of the sextic terms could be

derived with the exception of two off-diagonal constants (*h*<sub>2</sub> and *h*<sub>3</sub>) that have been kept fixed at the corresponding computed values. The agreement between the experimental and the theoretical parameters is very good. Minor discrepancies are only found for the sextic constants and the value of the *b* component of the dipole moment, experimentally determined by means of Stark effect measurements.<sup>20</sup> The latter discrepancy is unexpected based on the literature available for dipole moment computations.<sup>51</sup> Indeed, even though the use of augmented basis sets is recommended to ensure a flexible description of the outer valence region, the accuracy of the fc-CCSD(T)/cc-pVTZ values is expected to be on the order of the experimental uncertainty.<sup>51</sup> In this respect, the limitations of the experimental determination deserve to be mentioned; in fact, for both *μ*<sub>a</sub> and *μ*<sub>b</sub>, this relied on the observation of only two transitions.

The improvement in the accuracy of our spectroscopic parameters over those derived by Sugie et al.<sup>20</sup> is evident (see Table 1); the uncertainty on the rotational and centrifugal distortion constants has been reduced by 1 and 2 orders of magnitude, respectively. Owing to the higher *J* and *K*<sub>a</sub> values sampled in our measurements and to the enlarged data set available, centrifugal analysis of *trans*-2-azabutadiene has been
<table><caption>Table 3. Rotational and Quartic Centrifugal Distortion Constants of the Cis and Gauche Species<sup>a</sup></caption><thead><tr><th rowspan="2">constant</th><th rowspan="2">units</th><th colspan="2">butadiene<sup>b</sup></th><th rowspan="2">acrolein<sup>c</sup></th><th rowspan="2">2-azabutadiene</th></tr><tr><th>0<sup>+</sup> state</th><th>0<sup>-</sup> state</th></tr></thead><tbody><tr><td>A</td><td>MHz</td><td>21 223.0467(18)</td><td>21 232.3005(80)</td><td>22 831.6513(51)</td><td>24 021.128</td></tr><tr><td>B</td><td>MHz</td><td>5671.39702(56)</td><td>5667.6405(10)</td><td>6241.02408(38)</td><td>6071.991</td></tr><tr><td>C</td><td>MHz</td><td>4577.35202(47)</td><td>4581.14000(70)</td><td>4902.23005(24)</td><td>5120.829</td></tr><tr><td>D<sub>J</sub></td><td>kHz</td><td>4.630(34)</td><td>4.630(34)</td><td>4.893610(178)</td><td>5.42141</td></tr><tr><td>D<sub>JK</sub></td><td>kHz</td><td>-6.99(24)</td><td>97.67(40)</td><td>-27.8640(13)</td><td>-40.8174</td></tr><tr><td>D<sub>K</sub></td><td>kHz</td><td>114.65(32)</td><td>1039.6(76)</td><td>106.63(13)</td><td>207.669</td></tr><tr><td>d<sub>1</sub></td><td>kHz</td><td>-1.2045(112)</td><td>-1.932(28)</td><td>-1.480975(128)</td><td>-1.82008</td></tr><tr><td>d<sub>2</sub></td><td>kHz</td><td>0.635(41)</td><td>-3.57(17)</td><td>-0.109707(78)</td><td>-0.199939</td></tr></tbody></table>

<sup>a</sup>Values in parentheses are one standard deviation and refer to the last digits. <sup>b</sup>From ref 8. <sup>c</sup>From ref 16. Converted from A to S reduction.

expanded and now includes distortion terms up to the sextic constants.

Concerning the *gauche* form of 2-azabutadiene, our ab initio estimate of the spectroscopic parameters should provide reliable spectral predictions with an overall uncertainty on transition frequencies of about 0.1–0.2%. This deviation can be further reduced by applying an empirical, well-tested, scaling procedure based on the ratio between the computed and the experimental constants for the *trans* form

$$p_{\text{scaled}}^{(\text{gauche})} = p_{\text{theo}}^{(\text{gauche})} \times \left( \frac{p_{\text{exp}}^{(\text{trans})}}{p_{\text{theo}}^{(\text{trans})}} \right) \quad (3)$$

where $p^{(X)}$ represents the generic parameters of the X isomer. The scaled parameters are collected in the last column of Table 1, and their accuracy is expected to be increased by 1 order of magnitude with respect to the computed values.

## DISCUSSION AND CONCLUSIONS

Given the similar heavy-atom skeleton and the conjugated π system shared by 2-azabutadiene, butadiene, acrolein, and 1-azabutadiene, it is quite interesting to qualitatively compare the values of the rotational and centrifugal distortion constants available for each species. This comparison is performed separately for the *trans* and *cis/gauche* species, the results being reported in Tables 2 and 3, respectively. As far as the *trans* species are concerned, the similarity is evident. The order of magnitude and the sign of each constant listed in Table 2 are the same for all of the species considered. This indicates a similar mass distribution (i.e., molecular geometry) and stiffness of the chemical bonds due to the analogy between the conjugated systems.

The case of the *gauche/cis* species is slightly different. Here, the similarity is less pronounced due to two facts. First, the higher energy conformer of acrolein displays a planar cis geometry, while both butadiene and 2-azabutadiene possess nonplanar gauche forms. Second, *gauche*-butadiene exhibits tunneling between two equivalent minima, whereas no large amplitude motion occurs in the heterodienes. That being said, it is likely that the differences observed for the two inversion states of butadiene are due to an unsatisfactory treatment of the Coriolis interaction occurring between the 0⁺ and the 0⁻ states. In this framework, the centrifugal distortion constants of acrolein and 2-azabutadiene might be used as a reference to guide the analysis toward a better modeling.

Centrifugal analysis of 2-azabutadiene has another important outcome. The rotational spectrum of the most stable *trans* form now can be predicted with high accuracy in a wide range

of frequencies. This represents a key prerequisite for enabling astronomical searches of this species in dense spectral line surveys of molecular-rich objects, where an unambiguous spectroscopic identification typically requires a match within a few tens of kilohertz.

One isomer of 2-azabutadiene, namely, ethyl cyanide, has been observed in high-mass star-forming regions with remarkable high abundance,<sup>52–55</sup> while another isomer, i.e., 1-azabutadiene, is reported as tentatively detected in G +0.693.<sup>28</sup> These discoveries tell us that molecules showing the same degree of complexity as 2-azabutadiene can be formed efficiently in interstellar conditions, thus making the title heterodiene a good candidate for future detection. Besides giving information about the kinetics and thermodynamics of the C<sub>3</sub>H<sub>5</sub>N isomeric family, 2-azabutadiene can be used— together with the other heterodienes—to proxy the abundance of interstellar butadiene or even ethene. While we do not have yet any evidence that this suggestion might work, indirect support to this idea is provided by the relative abundance of CCH over CN cyclic-hydrocarbon derivatives in TMC-1, which are highly correlated.<sup>56</sup> In work by Cernicharo et al.,<sup>56</sup> it was also suggested that their relative abundance probably reflected the abundance of the CCH and CN radicals. In this respect, we note that 2-azabutadiene can be seen as substituted ethene, where the H<sub>2</sub>CN radical is linked to the latter through the N atom.

Finally, Titan’s chemistry deserves a note. Organic nitriles and, more generally, N-substituted hydrocarbons are largely present in the nitrogen-rich atmosphere Titan. In this view, spectroscopic characterization of 2-azabutadiene might be useful to gain further information on the chemistry occurring in what is considered to resemble the primitive Earth.

## ASSOCIATED CONTENT

### Supporting Information

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpca.2c00831.

Cartesian coordinates of the CBS+CV molecular structure for both *trans* and *gauche* forms; results of the NBO analysis; harmonic and anharmonic vibrational frequencies, IR intensities, and vibration-rotation α constants; ¹H and ¹³C NMR spectra of 2-[(2-chloroethyl)amino]acetonitrile and 2-azabutadiene; list of assigned rotational transitions together with the corresponding residuals from the final fit (PDF)
## AUTHOR INFORMATION

### Corresponding Authors

Mattia Melosso - Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, 40126 Bologna, Italy; Scuola Superiore Meridionale, Università di Napoli Federico II, 80138 Naples, Italy; orcid.org/0000-0002-6492-5921; Email: mattia.melosso2@unibo.it

Cristina Puzzarini - Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, 40126 Bologna, Italy; orcid.org/0000-0002-2395-8532; Email: cristina.puzzarini@unibo.it

### Authors

Ningjing Jiang - Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, 40126 Bologna, Italy

Luca Bizzocchi - Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, 40126 Bologna, Italy

Silvia Alessandrini - Scuola Normale Superiore, S6126 Pisa, Italy

Jean-Claude Guillemin - University of Rennes, Ecole Nationale Supérieure de Chimie de Rennes, CNRS, ISCR-UMR6226, F-35000 Rennes, France; orcid.org/0000-0002-2929-057X

Luca Dore - Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, 40126 Bologna, Italy; orcid.org/0000-0002-1009-7286

Complete contact information is available at: https://pubs.acs.org/10.1021/acs.jpca.2c00831

### Notes

The authors declare no competing financial interest.

## ACKNOWLEDGMENTS

This study was supported by the University of Bologna (RFO funds) as well as by the Italian Space Agency (ASI; "Life in Space" project, N. 2019-3-U.0). N.J. thanks the China Scholarships Council (CSC) for financial support. J.C.G. thanks the Centre National d'Etudes Spatiales (CNES) and the Program "Physique et Chimie du Milieu Interstellaire" (INSU-CNRS) for grants.

## REFERENCES

(1) Wiberg, K. B.; Rablen, P. R.; Baraban, J. H. Butadiene and heterodienes revisited. *J. Org. Chem.* **2018**, *83*, 8473-8482.

(2) Barone, V.; Alessandrini, S.; Biczysko, M.; Cheeseeman, J. R.; Clary, D. C.; McCoy, A. B.; DiRisio, R. J.; Neese, F.; Melosso, M.; Puzzarini, C. Computational molecular spectroscopy. *Nat. Rev. Methods Primers* **2021**, *1*, 38.

(3) Martin-Drumel, M.-A.; Porterfield, J. P.; Goubet, M.; Asselin, P.; Georges, R.; Soulard, P.; Nava, M.; Changala, P. B.; Billinghurst, B.; Pirali, O.; et al. Synchrotron-based High Resolution Far-infrared Spectroscopy of trans-butadiene. *J. Phys. Chem. A* **2020**, *124*, 2427-2435.

(4) Craig, N. C.; Hanson, K. A.; Pierce, R. W.; Saylor, S. D.; Sams, R. L. Rotational analysis of bands in the high-resolution infrared spectra of the three species of butadiene-1,4-d₂; refinement of the assignments of the vibrational fundamentals. *J. Mol. Spectrosc.* **2004**, *228*, 401-413.

(5) Craig, N. C.; Davis, J. L.; Hanson, K. A.; Moore, M. C.; Weidenbaum, K. J.; Lock, M. Analysis of the rotational structure in bands in the high-resolution infrared spectra of butadiene and butadiene-2,3-d₂; refinement in assignments of fundamentals. *J. Mol. Struct.* **2004**, *695-696*, 59-69.

(6) Craig, N. C.; Hanson, K. A.; Moore, M. C.; Sams, R. L. Rotational analysis of several bands in the high-resolution infrared spectrum of butadiene-1-¹³C₁: assignment of vibrational fundamentals. *J. Mol. Struct.* **2005**, *742*, 21-29.

(7) Craig, N. C.; Moore, M. C.; Patchen, A. K.; Sams, R. L. Analysis of rotational structure in the high-resolution infrared spectrum and assignment of vibrational fundamentals of butadiene-2,3-¹³C₂. *J. Mol. Spectrosc.* **2006**, *235*, 181-189.

(8) Baraban, J. H.; Martin-Drumel, M.-A.; Changala, P. B.; Eibenberger, S.; Nava, M.; Patterson, D.; Stanton, J. F.; Ellison, G. B.; McCarthy, M. C. The Molecular Structure of gauche-1,3-Butadiene: Experimental Establishment of Non-planarity. *Angew. Chem.* **2018**, *130*, 1839-1843.

(9) Puzzarini, C.; Penocchio, E.; Biczysko, M.; Barone, V. Molecular structure and spectroscopic signatures of acrolein: theory meets experiment. *J. Phys. Chem. A* **2014**, *118*, 6648-6656.

(10) Fine, J.; Goldstein, J.; Simmons, J. Microwave Spectrum of s-trans Acrolein. *J. Chem. Phys.* **1955**, *23*, 601-601.

(11) Wagner, R.; Fine, J.; Simmons, J.; Goldstein, J. Microwave Spectrum, Structure, and Dipole Moment of s-trans Acrolein. *J. Chem. Phys.* **1957**, *26*, 634-637.

(12) Cherniak, E.; Costain, C. Microwave Spectrum and Molecular Structure of trans-Acrolein. *J. Chem. Phys.* **1966**, *45*, 104-110.

(13) Winnewisser, M.; Winnewisser, G.; Honda, T.; Hirota, E. Ground State Centrifugal Distortion Constants of Trans-Acrolein, CH₂=CH-CHO from the Microwave and Millimeter Wave Rotational Spectra. *Z. Naturforsch. A* **1975**, *30*, 1001-1014.

(14) Blom, C.; Bauder, A. Microwave spectrum, rotational constants and dipole moment of s-cis acrolein. *Chem. Phys. Lett.* **1982**, *88*, 55-58.

(15) Evangelisti, L.; Maris, A.; Grieco, F.; Calabrese, C.; Melandri, S. Millimeter wave free-jet spectrum of acrolein and several isotopologues. *Can. J. Phys.* **2020**, *98*, 555-559.

(16) Daly, A. M.; Bermúdez, C.; Kolesniková, L.; Alonso, J. Comprehensive analysis of prebiotic propenal up to 660 GHz. *Astrophys. J. Suppl. Ser.* **2015**, *218*, 30.

(17) Penn, R. E. Microwave spectrum of 2-propene-1-imine, CH₂=CH-CH=NH. *J. Mol. Spectrosc.* **1978**, *69*, 373-382.

(18) Brown, R. D.; Godfrey, P. D.; Winkler, D. A. Hyperfine interactions in the microwave spectrum of 2-propen-1-Imine (vinyl-imine). *Chem. Phys.* **1981**, *59*, 243-247.

(19) Hamada, Y.; Tsuboi, M.; Matsuzawa, T.; Yamanouchi, K.; Kuchitsu, K.; Koga, Y.; Kondo, S. Pyrolysis of amines: Infrared spectrum of allylimine. *J. Mol. Spectrosc.* **1984**, *105*, 453-464.

(20) Sugie, M.; Takeo, H.; Matsumura, C. A study of the thermal decomposition and dehydrochlorination of N-chloroazetidine. Microwave spectra of N-chloromethylenimine, 1-azetine, and 2-azabutadiene. *J. Am. Chem. Soc.* **1989**, *111*, 906-910.

(21) Amatatsu, Y.; Hamada, Y.; Tsuboi, M. FTIR detection of unstable molecules: infrared spectrum of 2-azabutadiene. *J. Mol. Spectrosc.* **1987**, *123*, 276-285.

(22) Gratier, P.; Majumdar, L.; Ohishi, M.; Roueff, E.; Loison, J.; Hickson, K.; Wakelam, V. A new reference chemical composition for TMC-1. *Astrophys. J. Supp.* **2016**, *225*, 25.

(23) Sakai, N.; Shiino, T.; Hirota, T.; Sakai, T.; Yamamoto, S. Long Carbon-chain molecules and their anions in the starless core, Lupus-1A. *Astrophys. J. Lett.* **2010**, *718*, L49.

(24) Turner, B. U93.174 - A new interstellar line with quadrupole hyperfine splitting. *Astrophys. J.* **1974**, *193*, L83-L87.

(25) Agúndez, M.; Marcelino, N.; Cernicharo, J. Discovery of interstellar isocyanogen (CNCN): Further evidence that dicyanopolylynes are abundant in space. *Astrophys. J. Lett.* **2018**, *861*, L22.

(26) Burkhardt, A. M.; Loomis, R. A.; Shingledecker, C. N.; Lee, K. L. K.; Remijan, A. J.; McCarthy, M. C.; McGuire, B. A. Ubiquitous aromatic carbon chemistry at the earliest stages of star formation. *Nat. Astron.* **2021**, *5*, 181-187.

(27) Manigand, S.; Coutens, A.; Loison, J.-C.; Wakelam, V.; Calcutt, H.; Müller, H.; Jørgensen, J. K.; Taquet, V.; Wampfler, S.; Bourke, T.; et al. The ALMA-PILS survey: first detection of the unsaturated 3-
carbon molecules Propenal (C₂H₃CHO) and Propylene (C₃H₆) towards IRAS 16293-2422 B. *Astron. Astrophys.* **2021**, *645*, A53.

(28) Melosso, M.; Bizzocchi, L.; Jiang, N.; Alberton, D.; Caselli, P.; Rivilla, V. M.; Pietropolli Charmet, A.; Dore, L.; Puzzarini, C. Laboratory Spectroscopy of Allylimine and its Tentative Detection in the Interstellar Medium. *International Symposium on Molecular Spectroscopy (Virtual)*, University of Illinois at Urbana-Champaign, June 21-25, 2021; University of Illinois, 2021.

(29) Endres, C. P.; Schlemmer, S.; Schilke, P.; Stutzki, J.; Müller, H. S. The Cologne Database for Molecular Spectroscopy, CDMS, in the virtual atomic and molecular data centre, VAMDC. *J. Mol. Spectrosc.* **2016**, *327*, 95-104.

(30) Cologne Database for Molecular Spectroscopy (CDMS): https://cdms.astro.uni-koeln.de/.

(31) Jet Propulsion Laboratory (JPL) catalog: https://spec.jpl.nasa.gov/.

(32) Mills, I. M. In *Molecular Spectroscopy: Modern Research*; Rao, K. N., Mathews, C. W., Eds.; Academic Press: New York, 1972; Chapter 3.2, pp 115-140.

(33) Puzzarini, C.; Stanton, J. F.; Gauss, J. Quantum-chemical calculation of spectroscopic parameters for rotational spectroscopy. *Int. Rev. Phys. Chem.* **2010**, *29*, 273-367.

(34) Dunning, T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron Through Neon and Hydrogen. *J. Chem. Phys.* **1989**, *90*, 1007-1023.

(35) Moller, C.; Plesset, M. S. Note on an approximation treatment for many-electron systems. *Phys. Rev.* **1934**, *46*, 618.

(36) Puzzarini, C.; Heckert, J.; Gauss, J. The accuracy of rotational constants predicted by high-level quantum-chemical calculations. I. molecules containing first-row atoms. *J. Chem. Phys.* **2008**, *128*, 194108.

(37) Alessandrini, S.; Gauss, J.; Puzzarini, C. Accuracy of Rotational Parameters Predicted by High-Level Quantum-Chemical Calculations: Case Study of Sulfur-Containing Molecules of Astrochemical Interest. *J. Chem. Theory Comput.* **2018**, *14*, 5360-5371.

(38) Heckert, M.; Kállay, M.; Tew, D. P.; Klopper, W.; Gauss, J. Basis-Set Extrapolation Techniques for the Accurate Calculation of Molecular Equilibrium Geometries Using Coupled-Cluster Theory. *J. Chem. Phys.* **2006**, *125*, 044108.

(39) Raghavachari, K.; Trucks, G. W.; Pople, J. A.; Head-Gordon, M. A fifth-order perturbation comparison of electron correlation theories. *Chem. Phys. Lett.* **1989**, *157*, 479-483.

(40) Feller, D. The use of systematic sequences of wave functions for estimating the complete basis set, full configuration interaction limit in water. *J. Chem. Phys.* **1993**, *98*, 7059-7071.

(41) Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-set convergence of correlated calculations on water. *J. Chem. Phys.* **1997**, *106*, 9639-9646.

(42) Woon, D. E.; Dunning, T. H. Gaussian basis sets for use in correlated molecular calculations. V. Core-valence basis sets for boron through neon. *J. Chem. Phys.* **1995**, *103*, 4572-4585.

(43) Stanton, J. F.; Gauss, J.; Harding, M. E.; Szalay, P. G. CFOUR. A quantum chemical program package. 2016. For the current version, see http://www.cfour.de.

(44) Matthews, D. A.; Cheng, L.; Harding, M. E.; Lipparini, F.; Stopkowicz, S.; Jagau, T.-C.; Szalay, P. G.; Gauss, J.; Stanton, J. F. Coupled-cluster techniques for computational chemistry: The CFOUR program package. *J. Chem. Phys.* **2020**, *152*, 214108.

(45) Guillemin, J.-C.; Denis, J.-M. Flash vacuum thermolysis of α-aminonitriles and subsequent HCN removal on solid base, a 'one line' multistep sequence to reactive N-methyleneamines. *J. Chem. Soc., Chem. Comm.* **1985**, 951-952.

(46) Guillemin, J.-C.; Denis, J.-M. Synthèse d'imines linéaires non-stabilisées par réactions gaz-solide sous vide. *Tetrahedron* **1988**, *44*, 4431-4446.

(47) Melosso, M.; Conversazioni, B.; Esposti, C. D.; Dore, L.; Cane, E.; Tamassia, F.; Bizzocchi, L. The Pure Rotational Spectrum of ¹⁵ND₂ Observed by Millimetre and Submillimetre-Wave Spectroscopy. *J. Quant. Spectrosc. Ra.* **2019**, 222-223, 186-189.

(48) Melosso, M.; Bizzocchi, L.; Tamassia, F.; Degli Esposti, C.; Cané, E.; Dore, L. The rotational spectrum of ¹⁵ND. Isotopic-independent Dunham-type analysis of the imidogen radical. *Phys. Chem. Chem. Phys.* **2019**, *21*, 3564-3573.

(49) Pickett, H. M. The fitting and prediction of vibration-rotation spectra with spin interactions. *J. Mol. Spectrosc.* **1991**, *148*, 371-377.

(50) Watson, J. K. G. In *Vibrational Spectra and Structure*; Durig, J. R., Ed.; Elsevier: New York, 1977; pp 1-89.

(51) Helgaker, T.; Jørgensen, P.; Olsen, J. *Electronic-Structure Theory*; Wiley: Chichester, 2000.

(52) Johnson, D.; Lovas, F.; Gottlieb, C.; Gottlieb, E.; Litvak, M.; Guelin, M.; Thaddeus, P. Detection of interstellar ethyl cyanide. *Astrophys. J.* **1977**, *218*, 370-376.

(53) Margulés, L.; Motiyenko, R.; Demyk, K.; Tercero, B.; Cernicharo, J.; Sheng, M.; Weidmann, M.; Gripp, J.; Mäder, H.; Demaison, J. Rotational spectrum of deuterated and ¹⁵N ethyl cyanides: CH₃CHDCN and CH₂DCH₂CN and of CH₃CH₂C¹⁵N. *Astron. Astrophys.* **2009**, *493*, 565-569.

(54) Margulés, L.; Belloche, A.; Müller, H.; Motiyenko, R.; Guillemin, J.-C.; Garrod, R.; Menten, K. Spectroscopic study and astronomical detection of doubly ¹³C-substituted ethyl cyanide. *Astron. Astrophys.* **2016**, *590*, A93.

(55) Zeng, S.; Jiménez-Serra, I.; Rivilla, V.; Martín, S.; Martín-Pintado, J.; Requena-Torres, M.; Armijos-Abendaño, J.; Riquelme, D.; Aladro, R. Complex organic molecules in the Galactic Centre: the N-bearing family. *Mon. Not. R. Astron. Soc.* **2018**, *478*, 2962-2975.

(56) Cernicharo, J.; Agúndez, M.; Kaiser, R. I.; Cabezas, C.; Tercero, B.; Marcelino, N.; Pardo, J. R.; de Vicente, P. Discovery of two isomers of ethynyl cyclopentadiene in TMC-1: Abundances of CCH and CN derivatives of hydrocarbon cycles. *Astron. Astrophys.* **2021**, *655*, L1.