Next Article in Journal
Therapeutic Potential of Albumin Nanoparticles Encapsulated Visnagin in MDA-MB-468 Triple-Negative Breast Cancer Cells
Previous Article in Journal
Recent Advances in the Transition-Metal-Free Synthesis of Quinazolines
Previous Article in Special Issue
Experimental, Spectroscopic, and Computational Insights into the Reactivity of “Methanal” with 2-Naphthylamines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers

by
Silvia Alessandrini
1,*,
Mattia Melosso
1,
Víctor M. Rivilla
2,
Luca Bizzocchi
1 and
Cristina Puzzarini
1,*
1
Dipartimento di Chimica “Giacomo Ciamician”, Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy
2
Centro de Astrobiología (CSIC-INTA), Ctra. de Ajalvir km. 4, Torrejón de Ardoz, 28850 Madrid, Spain
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(7), 3226; https://doi.org/10.3390/molecules28073226
Submission received: 11 March 2023 / Revised: 30 March 2023 / Accepted: 31 March 2023 / Published: 4 April 2023
(This article belongs to the Special Issue Advances in Computational Spectroscopy)

Abstract

:
The C3H3NO family of isomers is relevant in astrochemistry, even though its members are still elusive in the interstellar medium. To identify the best candidate for astronomical detection within this family, we developed a new computational protocol based on the minimum-energy principle. This approach aims to identify the most stable isomer of the family and consists of three steps. The first step is an extensive investigation that characterizes the vast number of compounds having the C3H3NO chemical formula, employing density functional theory for this purpose. The second step is an energy refinement, which is used to select isomers and relies on coupled cluster theory. The last step is a structural improvement with a final energy refinement that provides improved energies and a large set of accurate spectroscopic parameters for all isomers lying within 30 kJ mol−1 above the most stable one. According to this protocol, vinylisocyanate is the most stable isomer, followed by oxazole, which is about 5 kJ mol−1 higher in energy. The other stable species are pyruvonitrile, cyanoacetaldehyde, and cyanovinylalcohol. For all of these species, new computed rotational and vibrational spectroscopic data are reported, which complement those already available in the literature or fill current gaps.

Graphical Abstract

1. Introduction

The molecular species observed so far in the interstellar medium (ISM) are generally small and light. In fact, based on the 2021 census of the interstellar and circumstellar molecules [1], and excluding the fullerene species C60, C60+, and C70, the average number of atoms contained in the molecules on this list ranges between 5 and 6. However, it has to be noted that a relevant number of discoveries have occurred since then and updated lists can be found in the Cologne Database for Molecular Spectroscopy [2] or on the Astrochymist website [3]. Despite the limited number of atoms and the severe constraints on chemical reactivity imposed by the ISM conditions (i.e., very low temperature and density, and ionizing radiation), interstellar molecules exhibit a large variety of functional groups, including nitriles [4,5,6], imines [7,8,9,10,11,12,13], amines [14,15], alcohols [16,17,18], amides [19,20,21,22,23], aldehydes [24,25,26,27], and many others. Indeed, several chemical linkages are possible within a specific chemical formula so that a generic isomeric family can be constituted by tens—if not hundreds—of structural isomers, each characterized by different functional groups. Therefore, the recognition of the specific members of an isomeric family that are found, or can be found, in the ISM is typically a non-trivial task.
In this respect, the so-called Minimum Energy Principle (MEP) [28] offers a powerful tool. It states that (i) the greater the thermochemical stability of a structural isomer, the higher the expected abundance in the ISM, and (ii) the ratio between the abundances of two isomers is governed by their relative stabilities. With very few exceptions, such as propadienone (H2CCCO) [29] and acetic acid (CH3COOH) [28], the MEP holds. It also represents a useful guide for selecting molecular species to study in the laboratory and subsequently hunt for in the ISM; it is a solid basis for rationalizing the observed molecular abundances.
When this principle was first formulated, the question of relative stability was addressed by employing quantum-chemical calculations rooted in density functional theory (DFT), using the B3LYP method in conjunction with a basis set of triple- ζ quality to compute both the electronic energy and the zero-point vibrational energy (ZPE) contribution [30]. The rationale behind this choice was the capability of this level of theory to provide accurate molecular structures and energies, and the fact that a comparison with higher-level coupled-cluster calculations for a selection of 14 isomeric families showed overall agreement between the two sets of results. However, according to Reference [31], the uncertainty affecting B3LYP energetics can be as large as 15.5 kJ mol−1 when the (aug)-def2-QZVP basis set is employed. This value is two times larger than that expected in Reference [28], where a smaller basis was employed. Therefore, the B3LYP method is only appropriate when the energy difference among the isomers is quite large. This difference is likely more pronounced for smaller species and decreases as the molecular size increases, as more (and more similar) ways to link the atoms within the molecule become possible. Thus, medium-sized species require accurate methodologies for predicting the correct relative stability among structural isomers, as the energy difference between them can be as small as a few kJ mol−1. Moreover, the evaluation of the ZPE correction can be as important as the electronic energy, and its contribution may be decisive for the identification of the most stable isomeric [32] or conformational [33] form. Therefore, it is evident that a more stringent protocol needs to be formulated for the precise identification of the stability order within a specific isomeric family. Furthermore, the limited sizes of interstellar molecules allow pushing the computational effort to the limit.
In this work, we introduce a multi-step approach for the identification of suitable candidates, within a given isomeric family, for astronomical investigations. This protocol exploits feasible yet accurate calculations for the preliminary search of all possible isomers (and their conformers) and then proceeds to refine their energetics using a higher level of theory, thereby aiming at the selection of the most stable species. For these, even higher-level computations are then performed to further improve the structural determination as well as the energetics, with the final aim of identifying and characterizing the best candidate for astronomical searches. As a test case, we present the study of the C3H3NO isomeric family, which is thus far elusive to astronomical observations, with the aim of identifying which isomer has the highest chance of being detected in the ISM.
Astronomical observations require accurate knowledge of the rotational spectra of the species under investigation. To date, spectroscopic studies have been reported for cyanoacetaldehyde [34], cyanooxirane [35], isoxazole [36,37,38], oxazole [36,39,40], pyruvonitrile [41], and vinylisocyanate [42,43,44,45]. While measurements of the rotational spectra have been extended up to the submillimeter-wave region for some of these species [34,35,41,45], experimental data are limited to the microwave domain (i.e., below 40 GHz) for the aromatic rings, i.e., oxazole and isoxazole, and are even totally missing for other species, such as cyanovinylalcohol. Therefore, in addition to evaluating the energetics of these isomers, we provide accurate estimates of rotational, centrifugal distortion, and nuclear quadrupole coupling constants for the most interesting species. Furthermore, we review the current state of available spectroscopic information to identify possible data gaps that might prevent the detection of C3H3NO species in interstellar space.
The manuscript is organized as follows. In the following section, we present and discuss the results, starting with the preliminary energetic investigation of the C3H3NO family. Next, we report the outcomes of the second step of the protocol, which leads to the selection of the isomers for further consideration. Subsequently, we focus on the third step, which involves further structural determinations and improvements in the energetics. Section 2 concludes with the spectroscopic characterization of the most promising candidates for astronomical investigations. All details of the multi-step protocol introduced here are provided in Section 3 (Materials and Methods). Finally, the concluding remarks are reported.

2. Results and Discussion

In this section, the results obtained for the C3H3NO family of isomers by exploiting the multi-step protocol mentioned in the Introduction (for details, see Section 3) are presented and discussed. Although the family of C3H3NO isomers has already been studied in previous works [46,47], an exhaustive characterization of the entire isomeric family is still missing. Therefore, a systematic search of all possible isomeric forms was carried out using the SciFinder-n application. Bicyclic or carbene-like structures and species with separation of charge were filtered out, as our goal was to identify the most stable isomers. The numbering used in the following is based on the initial discovery of the isomers and does not follow any energetic scale.

2.1. Step 1: Preliminary Investigation

The starting point of our selection procedure consisted of 42 species, including the E and Z forms of structural isomers whenever applicable. In addition, other stable conformers were manually generated, resulting in a total of 67 species considered in this preliminary step. For all of them, the structure was optimized using the double-hybrid rev-DSD-PBEP86 functional [48] in conjunction with the jun-cc-pVTZ basis set [49,50], also including empirical dispersion (overall, this level is briefly denoted as revDSD/junTZ). For further details, please refer to Section 3. The revDSD/junTZ electronic energy, corrected for the harmonic ZPE (hZPE) at the same level of theory (obtained straightforwardly from harmonic force field, HFF, calculations), was used to derive the stability order of the C3H3NO isomers. The list, in decreasing order of stability, is presented in Figure 1, where the arbitrary numbering mentioned above identifies the isomers. The correspondence between the structure and label, together with revDSD/junTZ and revDSD/junTZ + hZPE relative energies, is reported in the Supporting Material (SM).
The most stable isomer at the revDSD/junTZ + hZPE level is the compound 5, i.e., trans-vinylisocyanate, which thus represents our energy reference. The latter is followed by its gauche conformer (5b) and pyruvonitrile (3), located 3.74 and 6.75 kJ mol−1 higher in energy than 5, respectively. These are the only isomers that lie below the relative energy threshold of 10 kJ mol−1.
The 10–50 kJ mol−1 energy range, represented by orange dots in Figure 1, is dominated by four isomers: cyanoacetaldehyde, oxazole, 3-hydroxy-2-propenenitrile, and 2-hydroxy-2-propenenitrile. The trans form of cyanoacetaldehyde (6) is located at 12.4 kJ mol−1, while its gauche form 6b is about 3 kJ mol−1 higher in energy, i.e., at 15.8 kJ mol−1. The relative energy of oxazole, labeled as 1, is 13.6 kJ mol−1, thus lying in between the two previously mentioned species. Moreover, 3-hydroxy-2-propenenitrile, also known as cyanovinylalcohol, presents both Z and E configurations, which in turn can show the trans (t) and gauche (g) forms, depending on the orientation of the -OH moiety with respect to the CC double bond. The revDSD/junTZ + hZPE energy scale predicts Zg-cyanovinylalcohol (9b) lying at 21.0 kJ mol−1, and followed by the Eg form (8b) at 30.3 kJ mol−1. The Et form of cyanovinylalcohol (8) is about 1 kJ mol−1 higher in energy than the Eg species, while the Zt conformer is the highest in energy due to steric effects and it is located at 34.5 kJ mol−1. Two conformers are also present for 2-hydroxy-2-propenenitrile (1-cyanoethanol); they are located at 46.3 kJ mol−1 (g form, 11) and 48.1 kJ mol−1 (t form, 11b).
Moving to higher energies, the four conformers of 3-iminopropenone (13, 13b, 13c, and 13d) are located between 50 and 100 kJ mol−1, which is the range represented by the green dots in Figure 1. These four species differ for (i) the orientation of the imino group, with respect to the cumulenic C=C=O moiety, and (ii) the orientation of the H atom with respect to the C=N bond. The 50–100 kJ mol−1 range also contains the most stable isocyanide compound of the C3H3NO family, i.e., acetyl isocyanide (22), lying at 65.1 kJ mol−1, thus about 58 kJ mol−1 higher in energy than its cyanide counterpart (3). Propiolamide (4) is located at 82 kJ mol−1 and is followed by the 4-member cycle 2(3H)-azetone (18) and 3-amino-1,2-propadienone (33), which lie at 89.5 kJ mol−1 and 90.3 kJ mol−1, respectively.
The light-blue dots of Figure 1 represent the 100–200 kJ mol−1 energy range. The first member of the C3H3NO family belonging to this group, i.e., N-ethenylideneformamide (23), is found at 104.8 kJ mol−1 and it is followed by isoxazole (2), which lies only 3 kJ mol−1 above 23. Other members of the group are (i) the isocyanic counterparts of the species mentioned before, such as isocyanoacetaldehyde (27), isocyanovinylalcohol (41, and 37, which are the Z and E forms, respectively), and 1-isocyanoethenol (12); (ii) 3- or 4-member cycles. In the latter group, the main component is cyanooxirane (7), located at 118.0 kJ mol−1, which is preceded in stability by 2(1H)-azetone (16). The last species that is worth mentioning is vinylcyanide (10), which is the OCN counterpart of 5.
Finally, the 200–450 kJ mol−1 energy range is considered. This includes 25 structures, which are represented by pink dots in Figure 1. Most of the species of this group (∼60 %) are characterized by a NO bond in the N=O or N−OH form. This suggests that nitroso and N-hydroxylamine groups are unlikely to be formed in astronomical environments. Indeed, the MEP suggests that either their formation is not favored or they exhibit an extremely low abundance of such species.
Seeking the most stable isomer within the C3H3NO family, it is necessary to improve the estimates obtained in the first step. From a computational perspective, this involves improving their electronic description and, consequently, the quantum-chemical methodology used to compute their energy. Since such an improvement can be computationally expensive, only the isomers predicted to lie below 100 kJ mol−1 were selected for additional computations, resulting in a total of 20 isomers.

2.2. Step 2: Energy Refinement

In this step, the energy of the selected 20 isomers is computed using a composite scheme based on the CCSD(T) method (coupled cluster singles and doubles approximation augmented by a perturbative treatment of triple excitations) [51]. This composite scheme accounts for extrapolation to the complete basis set (CBS) limit and the core-valence (CV) correlation effects. While the details of this composite scheme are provided in Section 3, it is abbreviated as CBS + CV hereafter. These energy computations were performed on top of the revDSD/junTZ optimized geometries, and always included the hZPE contribution at the same level of theory. A summary of the results is presented in Figure 2, which shows the CBS + CV + hZPE energies and the corresponding isomer labels.
The CBS + CV + hZPE approach still predicts t-vinylisocyanate as the most stable species, which is followed by its g conformer at 3.1 kJ mol−1 and then by oxazole at 3.4 kJ mol−1. Such a small energy difference between oxazole and t-vinylisocyanate compared to the revDSD/junTZ result, which is about 14 kJ mol−1, is remarkable. This different behavior has to be ascribed to the CBS + CV electronic energy, which predicts oxazole as the most stable isomer of the C3H3NO family, while the incorporation of the hZPE inverts the trend between the latter and t-vinylisocyanate. Moving to higher energies, oxazole is followed by pyruvonitrile (3), at 5.0 kJ mol−1, and then by cyanoacetaldehyde (6), at 10.2 kJ mol−1. The remaining isomers exhibit the same energetic order as that obtained from the revDSD/junTZ calculations. Overall, the CBS + CV energies differ in absolute terms from the revDSD counterparts by an average of 4.2 kJ mol−1, with the largest discrepancy observed for oxazole, which is predicted to be 10 kJ mol−1 lower in energy by CBS + CV + hZPE.
Useful information on potential candidates for radioastronomical detection is already obtained in this step of the protocol. The possibility of observing a given molecule in the ISM depends on the intensity of its rotational transitions. While the abundance of the molecule in the environment under consideration clearly plays a role, the rotational line intensity is proportional to the square of the electric dipole moment. Therefore, the value of the latter needs to be considered together with relative energy. For this reason, Figure 3 plots the total electric dipole moment for the 20 most stable isomers, computed at the revDSD/junTZ level, against their relative CBS + CV + hZPE energy. This figure points out that the most stable species of the family (namely, 5, 5b, and 1) have a small dipole moment in comparison with the other members. All forms of cyanovinylalcohol are favored by the large dipole moment, but pyruvonitrile (3) and the trans form of cyanoacetaldehyde (6b) seem to be those that better combine a low relative energy and a large dipole moment. The other species are very high in energy, and this is not counterbalanced by the magnitude of the dipole moment, making their detection unlikely. Interestingly, the gas-phase formation of cyanoacetaldehyde has been suggested as feasible in the ISM via the reaction between the CN radical and oxirane [52], both species being already detected in different interstellar regions (see Reference [52] for details). The same consideration applies to cyanovinylalcohol, which can be easily formed from the reaction between the CN radical and vinylalcohol, the latter being present in the ISM as well [18,53,54]. Pyruvonitrle is instead a co-product of the reaction between acetahldehyde and the CN radical [55].
The number of species of the C3H3NO family that are potentially observable in the ISM is narrowed by employing the more accurate energy list of Figure 2 and the dipole moment considerations from Figure 3. Therefore, only the species within 30 kJ mol−1, with respect to the most stable species (5), which are 10 isomers, have been selected for the next step of the procedure, i.e., the structural improvement.

2.3. Step 3: Structure Improvement and Final Energy Refinement

In this step, the molecular structures of the 10 selected isomers (blue-colored molecules in Figure 2) were refined to further improve the electronic energies and provide accurate estimates of the corresponding rotational parameters, especially for those molecules that could potentially be present in the ISM and have not been studied experimentally via rotational spectroscopy. To achieve this goal, the CBS + CV composite scheme is employed in geometry optimizations, yielding improved equilibrium rotational constants (which only depend on the molecular equilibrium structure) and electronic energies (obtained on top of geometries at the same level of theory). To predict the rotational constants accurately, the vibrational corrections to the equilibrium counterparts have to be computed, requiring anharmonic force field calculations that have been performed at the revDSD/junTZ level of theory (see Section 3 for details). These anharmonic computations also provide better estimates of the ZPE corrections (anharmonic ZPE, aZPE), which further improve the relative energy scale of isomer stability. While this section discusses the new energetic trend, the spectroscopic characterization is addressed in the next section.
Moving on from the CBS + CV + hZPE level (revDSD/junTZ structure) to the CBS + CV + aZPE one (CBS + CV geometry; see Table S2), the relative energy of each isomer remains nearly unchanged; this means that—as expected [56,57]—double-hybrid DFT geometries are suitable for high-level energy computations and the anharmonic contribution to the ZPE correction is small. Indeed, the latter accounts—on average—for less then 2 kJ mol−1. Oxazole is now predicted as the second most stable species, followed by the g-vinylisocyanate, which are at 3.2 kJ mol−1 and 3.3 kJ mol−1, respectively, above t-vinylisocyanate.
Since the last step of our computational protocol leads to oxazole and t-vinylisocyanate as the two most stable isomers lying close in energy, a final refinement is required to confirm (or not) the latter species as the most stable form of the C3H3NO family. For this purpose, the electronic energy was computed using the HEAT-like approach [58,59], which incorporates higher-order terms in the coupled cluster expansion and relativistic effects (see Section 3); moreover, the aZPE correction was improved by means of a hybrid approach, which combines the HFF at the fc-CCSD(T)/cc-pVTZ level with the revDSD/junTZ anharmonic contribution. At the CBS + CV/CBS + CV level (which means that the CBS + CV level has been used for both structure and energy), oxazole is located −4.52 kJ mol−1 below t-vinylisocyanate. The incorporation of the fT contribution (i.e., the full treatment of triple excitations) reduces the energy difference to −2.86 kJ mol−1, thus favoring t-vinylisocyanate by 1.66 kJ mol−1. The effect of the perturbative treatment of quadruple excitations (pQ) is counterbalanced by the DBOC term (which is the diagonal Born–Oppenheimer correction, see Section 3), thus confirming −2.86 kJ mol−1 as the energy gap. The scalar relativistic contribution further reduces the energy difference to −2.73 kJ mol−1. Overall, the HEAT-like approach predicts oxazole as the most stable isomer when considering only the electronic energy. However, the ZPE term has to be taken into account. Its incorporation as hZPE at the fc-CCSD(T)/cc-pVTZ level favors the trans form of vinylisocyanate by 8.4 kJ mol−1, thus locating oxazole at 5.6 kJ mol−1 above the former. The anharmonic correction to the hZPE contribution has a small impact on the overall value, leading to an energy difference of 5.4 kJ mol−1 and establishing t-vinylisocyanate as the most stable isomer of the C3H3NO family.
The energy gap between oxazole and t-vinylisocyanate is strongly affected by the vibrational energy, as observed for the two conformers of cyclopropanecarboxaldehyde [33]; this can be attributed to the high vibrational frequencies of oxazole, a stable aromatic species, in contrast to the flexible t-vinylisocyanate.

2.4. Spectroscopic Characterization

Several C3H3NO isomers have already been characterized by means of rotational spectroscopy. Rotational transitions have been reported for high-energy isomers, such as propiolamide (4) [60] and isoxazole (2) [36,37,38]. Similarly, six isomers among the ten most stables species (according to the present protocol), have also been experimentally observed. These are conformers of (i) vinylisocyanate [42,43,44], (ii) cyanoacetaldehyde [34], (iii) pyruvonitrile [41], and (iv) oxazole [36,39,40,61]. The only species that have not yet been experimentally characterized are the four conformers of cyanovinylalcohol, which are the remaining molecules in the top 10 of the stable isomers. Therefore, our protocol has been employed to derive accurate rotational spectroscopic parameters of the six isomers that have been experimentally studied, which can be used to estimate the uncertainties affecting our predictions for cyanovinylalcohol.
Table 1 reports the CBS + CV equilibrium rotational constants, corrected for vibrational contributions at the revDSD/junTZ level, for the species that have already been experimentally characterized. The A 0 rotational constant is reproduced with a relative error of 0.5%, while B 0 and C 0 have both a mean error of 0.1%. The average error on A 0 significantly reduces to 0.3% if one removes from the statistics the A 0 value for t-NCCH2CHO, which might be affected by experimental misassignments, as suggested in Reference [34]. These conservative estimates for uncertainty can be applied to the rotational constants obtained for the four isomers of cyanovinylalcohol collected in Table 2, where the centrifugal distortion and nuclear quadrupole coupling constants, together with the electric dipole moment components, are also reported.
The values of Table 2 can be compared with those of Reference [34], where the conformers of cyanovinylalcohol were investigated by means of the second-order Møller–Plesset perturbation theory (MP2) in conjunction with the aug-cc-pVTZ basis set. The energy differences reported in Reference [34] are consistent with our computations, which, however, predict that all conformers are planar and reveal the presence of a trans form for (E)-cyanovinylalcohol that was not reported in the previous study. The rotational constants predicted in Reference [34] refer to the equilibrium structure and exhibit deviations with respect to ours up to 100 MHz ( A 0 , (E)-g-cyanovinylalcohol). In most cases, the discrepancies are larger than the vibrational corrections, which should be included whenever aiming for accurate estimates [62,63], thus highlighting the inadequacy of the MP2 method in predicting rotational constants of spectroscopic accuracy. There is no comparison available for centrifugal distortion constants and nuclear quadrupole coupling constants, which are also essential parameters to ensure a correct assignment of rotational transitions.
The spectroscopic constants of the four forms of cyanovinylalcohol were used to simulate their rotational spectra, which are reported in Figure 4. The intensity of the transitions follows the opposite trend of stability, i.e., the most intense spectrum is observed for the less stable species. However, at room temperature, all of the spectra show a peak of intensity at around 300 GHz; this will be helpful for future laboratory studies.
For the 10 most stable isomers of the C3H3NO family, our computational methodology allowed for the accurate determination of anharmonic vibrational frequencies in addition to the rotational parameters, which can be useful in guiding the assignment of high-resolution infrared (IR) spectra. Table 3 reports our best anharmonic values for oxazole and t-vinylisocyanate, obtained with a hybrid approach that combines fc-CCSD(T)/cc-pVTZ harmonic frequencies with anharmonic revDSD/junTZ contributions. The computed anharmonic vibrational frequencies for oxazole reproduce experimental data with errors below 5 cm−1 in the 600–1400 cm−1 range, as evident from the comparison with available experimental data. For t-vinylisocyanate, vibrational frequencies have not been measured experimentally, and thus our values represent the most accurate data available in the literature.
For the other eight species lying below 30 kJ mol−1, the anharmonic revDSD/junTZ vibrational frequencies are reported in the SM. In this case, the comparison with experimental data can be performed only for pyruvonitrile. For this isomer, it is noted that the assignment of the band observed at 535 cm−1 to the ν 16 fundamental [64] casts some doubts. In fact, for the same mode, our calculations predict a value of 584.86 cm−1 and, thus, 50 cm−1 above the observed one. However, such a deviation is five times larger than the average deviation observed for the other fundamental bands.

3. Materials and Methods

In this work, a three-step protocol based on the MEP was introduced and exploited to identify potential candidates for astronomical detection among the isomers of the C3H3NO family. The procedure is schematized in Figure 5 and explained in detail in the subsequent sections. All DFT calculations were performed using the Gaussian16 software suite [65], while the CFour program package [66,67] and MRCC [68] (interfaced with CFour) were used for methodologies based on the coupled cluster (CC) approach [69].

3.1. Step 1: Preliminary Investigation

The first step of the procedure is the preliminary investigation of a large number of molecules having the same molecular formula, i.e., C3H3NO in this case. This step involves a large number of compounds and needs to be carried out using an affordable computational method that can provide reliable data. For this reason, as mentioned in Section 2, all species of the C3H3NO family were analyzed at the revDSD/junTZ level. All DFT computations were corrected for Grimme’s D3 empirical dispersion, including the BJ dumping function [70,71]. The jun-cc-pVTZ set is a partially augmented triple-zeta basis set that balances computational cost and the inclusion of diffuse functions. For each isomer, the revDSD/junTZ level of theory was used to obtain the equilibrium molecular structure, the corresponding electronic energy, and the analytic harmonic force field. The force field was used to ensure the nature of the stationary point and to provide the harmonic zero-point energy correction (hZPE) to the electronic energy.
The selection of initial structures to be considered is left to the user, but in this case, Sci-Finder was employed to obtain a large number of initial structures. An automatic procedure based on genetic algorithms guiding semi-empirical methods [72,73,74] is currently under development.

3.2. Step 2: Energy Refinement

For the isomers within a specific energy range (for example, 0–100 kJ mol−1, based on the preliminary investigation at step 1), the electronic energy is further improved by exploiting a composite scheme rooted in the CC theory on top of the revDSD/junTZ geometries. By definition, composite schemes include different contributions, each of them computed at the best possible level of theory according to the size of the system, in order to minimize errors associated with quantum-chemical calculations. In view of the number of isomers considered and their sizes, the present protocol uses the CBS + CV approach [33,56,75,76,77,78] for the electronic energy ( E ( CBS + CV ) ):
E ( CBS + CV ) = E HF - SCF + E CCSD ( T ) + Δ E ( CV )
The first term on the right-hand side of Equation (1) is the extrapolation to the CBS limit of the Hartree–Fock self-consistent field (HF-SCF) energy, which is evaluated with the exponential formula by Feller [79], requiring three energy computations. These are carried out with the cc-pVnZ family of basis sets, using n = T, Q, and 5 [49]. The second contribution in Equation (1) accounts for the CBS limit of the CCSD(T) correlation energy, which is estimated via the two-point n 3 formula [80], using the cc-pVTZ and cc-pVQZ basis sets. Since these computations are carried out within the frozen-core (fc) approximation, the last term on the right-hand side of Equation (1) incorporates the effects of correlating the inner-shell electrons, which is the CV correlation term. The Δ E ( CV ) is computed as the energy difference between all-electron (ae) and fc calculations, both carried out with the same basis set, cc-pCVTZ [81]. The CBS + CV electronic energies augmented for hZPE were then employed for a new classification of the C3H3NO isomers, and only those in the 0–30 kJ mol−1 energy range were retained in the third step.

3.3. Step 3: Structure Improvement & Final Energy Refinement

The third step of our computational protocol is used to obtain improved structural determinations with the aim of further improving the electronic energies, as well as laying the basis for an accurate spectroscopic characterization. For this purpose, the CBS + CV energy of Equation (1) is used to build an energy gradient ( d E ( CBS + CV ) / d x ), which is minimized to obtain the corresponding equilibrium geometry [82,83]:
d E ( CBS + CV ) d x = d E ( HF - SCF ) d x + d Δ E ( CCSD ( T ) ) d x + d Δ E ( CV ) d x .
In addition to the CBS + CV equilibrium structures, the above scheme provides CBS + CV electronic energies on top of the geometries optimized at the same level of theory. These energies, augmented for aZPE (CBS + CV + aZPE), are used to determine the final stability order of the C3H3NO isomers considered in step 3.
Specifically, this step was applied to the 10 isomers lying within 30 kJ mol−1 above the most stable species. The accurate CBS + CV equilibrium structures of these isomers straightforwardly provide their equilibrium rotational constants ( B e ). Using the vibrational perturbation theory to the second order (VPT2) [84], the rotational constants of the vibrational ground state consist of two terms:
B 0 γ = B e γ + Δ B 0 γ ,
where B 0 γ = A 0 , B 0 , and C 0 for γ = a , b , and c, respectively. As mentioned above, the B e γ term depends on the equilibrium molecular structure, but also on its atomic composition. The vibrational correction ( Δ B 0 ) is obtained as the half-sum over the normal modes (r) of the vibration–rotation interaction constants ( α r γ ). [84]:
Δ B 0 γ = 1 2 r α r γ .
The B e γ value accounts for about 99% of B 0 γ and, thus, strongly affects its accuracy [62,85]. Therefore, the equilibrium geometries of the species of interest have been evaluated at the CBS + CV level to ensure the accuracy of the B 0 value. To calculate the Δ B 0 term, an anharmonic force field is required. Although its contribution is small, the vibrational correction is evaluated at the revDSD/junTZ level due to its affordability. Additionally, these anharmonic computations provide other spectroscopic parameters for rotational spectroscopy such as the centrifugal distortion constants up to sextic terms, electric dipole moment, and nuclear quadrupole coupling constants, along with anharmonic vibrational frequencies and the corresponding ZPE contribution (aZPE) [86].
The final energy refinement is part of the third step of the protocol, but it is not mandatory. Instead, it should be employed whenever two or more relevant isomers are particularly close in energy so that a sub-kJ mol−1 accuracy is required to assure the relative energy. In this case, the electronic energy at the CBS + CV level ( E ( CBS + CV ) ) is augmented by higher-order terms to exploit the HEAT-like protocol (E(HEAT-like)) [58,59]:
E ( HEAT - like ) = E ( CBS + CV ) + Δ E ( fT ) + Δ E ( pQ ) + Δ E ( DBOC ) + Δ E ( rel ) .
The second term on the right-hand side of Equation (5) represents the contribution from the full treatment of triple excitations, which is computed as the energy difference between CCSDT [87,88] and CCSD(T) energies using the cc-pVTZ basis set. Here, CCSDT denotes the CC singles, doubles, and triples method. The term Δ E ( pQ ) takes into account the effect of quadruple excitations and is computed using the CCSDT(Q) method [89,90,91] (which is CCSDT augmented by a perturbative treatment of quadruples), employing the cc-pVDZ basis set. The remaining two terms account for the diagonal Born–Oppenheimer correction (DBOC), computed at the HF-SCF level using the aug-cc-pVTZ basis set [92], and the scalar relativistic term, computed at the CCSD(T) level in conjunction with the aug-cc-pCVTZ basis set. This last term includes only one-electron mass-velocity and Darwin corrections and is obtained using perturbative techniques [93,94].
The E ( HEAT - like ) energy needs to be corrected for the ZPE term, which is purposely improved by using the harmonic fc-CCSD(T)/cc-pVTZ data and the anharmonic contribution at the revDSD/junTZ level. This hybrid anharmonic force field has already been used in Reference [33] and allows retrieving accurate vibrational frequencies of the species considered in the final refinement.

4. Conclusions

The computational protocol developed in this study provides a systematic and simple methodology to investigate families of isomers of astrochemical interest using the MEP as a guiding reference.
The first step of the procedure involves a preliminary investigation of all possible isomers and conformers, which were 67 in the case of the C3H3NO family. This step is carried out using a double-hybrid DFT functional in conjunction with a partially augmented triple-zeta basis set, which is able to provide a reliable energy "ladder" when the hZPE correction is included. The results from this step are used to select a smaller group of species, for example, those within 100 kJ mol−1 above the most stable isomer, which are retained for the next step. In the second step, a more accurate computational methodology is used to improve the electronic energy of the selected isomers. To accomplish this, a composite scheme rooted in the CC theory, the CBS + CV approach, is used to compute the electronic energy while retaining the double-hybrid reference geometries and the corresponding hZPE. Moving from the first to the second step, the relative energy can change up to 10 kJ mol−1, as noted for oxazole, and several changes in the relative energy scale might occur. The third step of the protocol involves a limited number of species, i.e., those lying below 30 kJ mol−1, according to the results of the previous step. At this stage, the aim is to further improve the energetics while providing an accurate set of rotational and vibrational spectroscopic data. This is accomplished by optimizing the molecular structures at the CBS + CV level and extending the determination of the force field by incorporating the anharmonic terms at the revDSD/junTZ level. Overall, these computations provide a comprehensive set of rotational parameters, ranging from rotational to centrifugal distortion and nuclear quadrupole coupling constants, as well as anharmonic vibrational frequencies, which are useful for interpreting IR spectra. The third step of our protocol results in only small energy changes, with oxazole becoming the second most stable isomer, situated between two vinylisocyanate conformers. To establish the most stable species with sub-kJ mol−1 accuracy, a non-mandatory step of the procedure, referred to as final energy refinement, is carried out using the HEAT-like composite scheme. At this level, oxazole is predicted to be the most stable species by about 3 kJ mol−1, but the smaller zero-point energy contribution of t-vinylisocyanate, in comparison to that of the aromatic ring, establishes t-CH2CHNCO as the most stable isomer within the C3H3NO family.
In conclusion, based on the MEP, the first isomer of the C3H3NO family that should be detected in the ISM is the trans form of vinylisocyanate. To date, this species has been searched for in the G + 0.693-0.027 molecular cloud [95] and the Sgr B2(N) region [45]. In both cases, only upper limits on its abundance could be derived, which are 2.5 × 10 12 cm−2 and 2.4 × 10 16 cm−2 for G + 0.693 and Sgr B2(N), respectively. To the best of our knowledge, no searches have been reported towards the cold core of the Taurus Molecular Cloud (TMC-1), one of the richest and most studied dark clouds. However, the high sensitivity of the QUIJOTE [96] and GOTHAM [97] line surveys (with a noise level below 1 mK) should allow for the identification of trans vinylisocyanate even if its column density is as low as ∼ 10 11 cm−2. Furthermore, our protocol indicates that cyanovinylalcohol and cyanoacetaldehyde are potential candidates for detection in the ISM. However, laboratory spectroscopic studies are still missing for cyanovinylalcohol, thus our study provides a comprehensive set of rotational and vibrational spectroscopic parameters to assist future experimental work.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28073226/s1. The information includes: the list of revDSD/junTZ and revDSD/junTZ + hZPE relative energies, together with the corresponding isomeric labels and structures (Table S1). CBS + CV and CBS + CV + aZPE relative energies for the ten most stable isomers obtained on top of the CBS + CV geometries (Table S2). Optimized CBS + CV geometries for the ten most stable isomers in Cartesian coordinates (Tables S3–S7) together with the anharmonic frequencies computed at the revDSD/junTZ level for the same species. Harmonic and anharmonic ZPE corrections computed at the revDSD/junTZ level (Tables S8–S11).

Author Contributions

S.A. worked on conceptualization of the work as well as investigation. S.A., L.B. and M.M. carried out the formal analysis and wrote a first draft of the manuscript. V.M.R. and C.P. reviewed the manuscript and provided the funding. Final editing was carried out by C.P., V.M.R., L.B., M.M. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by MUR (PRIN grant number 202082CE3T) and by the University of Bologna (RFO funds). V.M.R. received support from the project RYC2020-029387-I funded by MCIN/AEI/10.13039/501100011033. The COST Action CA21101 “COSY-Confined molecular systems: from a new generation of materials to the stars’’ is also acknowledged.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article Supplementary Materials (for details see above), further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ISMinterstellar medium
MEPminimum energy principle
DFTdensity functional theory
ZPEzero-point vibrational energy
HF-SCFHartree–Fock self-consistent field
CCcoupled cluster
CCSD(T) coupled cluster with single, double, and perturbative triple excitations
CBScomplete basis set
CVcore-valence
CCSDTcoupled-cluster with single, double, and triple excitations
CCSDT(Q)coupled-cluster single, double, triple, and perturbative quadrupole excitations
DBOCdiagonal Born–Oppenheimer correction
MP2second-order Møller-Plesset perturbation theory
VPT2second-order vibrational perturbation theory (VPT2)
HFFharmonic force field

References

  1. McGuire, B.A. 2021 Census of Interstellar, Circumstellar, Extragalactic, Protoplanetary Disk, and Exoplanetary Molecules. Astrophys. J. Suppl. Ser. 2022, 259, 30. [Google Scholar] [CrossRef]
  2. CDMS, Cologne Database for Molecular Spectroscopy. Available online: cdms.astro.uni-koeln.de/classic/molecules (accessed on 1 March 2023).
  3. The Astrochymist: An Internet Resource for Astrochemists and Interested Bystanders. Available online: astrochymist.org (accessed on 1 March 2023).
  4. Garrod, R.; Belloche, A.; Müller, H.; Menten, K. Exploring molecular complexity with ALMA (EMoCA): Simulations of branched carbon-chain chemistry in Sgr B2 (N). Astron. Astrophys. 2017, 601, A48. [Google Scholar] [CrossRef] [Green Version]
  5. Calcutt, H.; Jørgensen, J.; Müller, H.; Kristensen, L.; Coutens, A.; Bourke, T.; Garrod, R.; Persson, M.; Van der Wiel, M.; Van Dishoeck, E.; et al. The ALMA-PILS survey: Complex nitriles towards IRAS 16293–2422. Astron. Astrophys. 2018, 616, A90. [Google Scholar] [CrossRef] [Green Version]
  6. Rivilla, V.M.; Jimenez-Serra, I.; Martín-Pintado, J.; Colzi, L.; Tercero, B.; De Vicente, P.; Zeng, S.; Martín, S.; García De La Concepción, J.; Bizzocchi, L.; et al. Molecular precursors of the RNA-world in space: New nitriles in the G + 0.693 − 0.027 molecular cloud. Front. Astron. Space Sci. 2022, 9, 160. [Google Scholar] [CrossRef]
  7. Godfrey, P.; Brown, R.; Robinson, B.; Sinclair, M. Discovery of interstellar methanimine (formaldimine). Astrophys. Lett. 1973, 13, 119. [Google Scholar]
  8. Zaleski, D.P.; Seifert, N.A.; Steber, A.L.; Muckle, M.T.; Loomis, R.A.; Corby, J.F.; Martinez, O.; Crabtree, K.N.; Jewell, P.R.; Hollis, J.M.; et al. Detection of E-cyanomethanimine toward Sagittarius B2 (N) in the Green Bank Telescope PRIMOS survey. Astrophys. J. Lett. 2013, 765, L10. [Google Scholar] [CrossRef]
  9. Loomis, R.A.; Zaleski, D.P.; Steber, A.L.; Neill, J.L.; Muckle, M.T.; Harris, B.J.; Hollis, J.M.; Jewell, P.R.; Lattanzi, V.; Lovas, F.J.; et al. The detection of interstellar ethanimine (CH3CHNH) from observations taken during the GBT PRIMOS survey. Astrophys. J. Lett. 2013, 765, L9. [Google Scholar] [CrossRef] [Green Version]
  10. Ligterink, N.F.W.; Calcutt, H.; Coutens, A.; Kristensen, L.; Bourke, T.; Drozdovskaya, M.N.; Müller, H.; Wampfler, S.; van Der Wiel, M.; Van Dishoeck, E.; et al. The ALMA-PILS survey: Stringent limits on small amines and nitrogen-oxides towards IRAS 16293–2422B. Astron. Astrophys. 2018, 619, A28. [Google Scholar] [CrossRef] [Green Version]
  11. Rivilla, V.; Martín-Pintado, J.; Jiménez-Serra, I.; Zeng, S.; Martín, S.; Armijos-Abenda no, J.; Requena-Torres, M.; Aladro, R.; Riquelme, D. Abundant Z-cyanomethanimine in the interstellar medium: Paving the way to the synthesis of adenine. Mon. Not. R. Astron. Soc. Lett. 2019, 483, L114–L119. [Google Scholar] [CrossRef] [Green Version]
  12. Bizzocchi, L.; Prudenzano, D.; Rivilla, V.M.; Pietropolli-Charmet, A.; Giuliano, B.M.; Caselli, P.; Martín-Pintado, J.; Jiménez-Serra, I.; Martín, S.; Requena-Torres, M.A.; et al. Propargylimine in the laboratory and in space: Millimetre-wave spectroscopy and its first detection in the ISM. Astron. Astrophys. 2020, 640, A98. [Google Scholar] [CrossRef]
  13. Alberton, D.; Bizzocchi, L.; Jiang, N.; Melosso, M.; Rivilla, V.; Charmet, A.P.; Giuliano, B.; Caselli, P.; Puzzarini, C.; Alessandrini, S.; et al. Laboratory spectroscopy of allylimine and tentative detection towards the G + 0.693 − 0.027 molecular cloud. Astron. Astrophys. 2023, 669, A93. [Google Scholar] [CrossRef]
  14. Kaifu, N.; Morimoto, M.; Nagane, K.; Akabane, K.; Iguchi, T.; Takagi, K. Detection of interstellar methylamine. Astrophys. J. 1974, 191, L135–L137. [Google Scholar] [CrossRef]
  15. Zeng, S.; Jiménez-Serra, I.; Rivilla, V.M.; Martín-Pintado, J.; Rodríguez-Almeida, L.F.; Tercero, B.; de Vicente, P.; Rico-Villas, F.; Colzi, L.; Martín, S.; et al. Probing the chemical complexity of amines in the ISM: Detection of vinylamine (C2H3NH2) and tentative detection of ethylamine (C2H5NH2). Astrophys. J. Lett. 2021, 920, L27. [Google Scholar] [CrossRef]
  16. Müller, H.S.; Belloche, A.; Xu, L.H.; Lees, R.M.; Garrod, R.T.; Walters, A.; van Wijngaarden, J.; Lewen, F.; Schlemmer, S.; Menten, K.M. Exploring molecular complexity with ALMA (EMoCA): Alkanethiols and alkanols in Sagittarius B2 (N2). Astron. Astrophys. 2016, 587, A92. [Google Scholar] [CrossRef] [Green Version]
  17. Rivilla, V.M.; Colzi, L.; Jiménez-Serra, I.; Martín-Pintado, J.; Megías, A.; Melosso, M.; Bizzocchi, L.; López-Gallifa, Á.; Martínez-Henares, A.; Massalkhi, S.; et al. Precursors of the RNA World in Space: Detection of (Z)-1,2-ethenediol in the Interstellar Medium, a Key Intermediate in Sugar Formation. Astrophys. J. Lett. 2022, 929, L11. [Google Scholar] [CrossRef]
  18. Jiménez-Serra, I.; Rodríguez-Almeida, L.F.; Martín-Pintado, J.; Rivilla, V.M.; Melosso, M.; Zeng, S.; Colzi, L.; Kawashima, Y.; Hirota, E.; Puzzarini, C.; et al. Precursors of fatty alcohols in the ISM: Discovery of n-propanol. Astron. Astrophys. 2022, 663, A181. [Google Scholar] [CrossRef]
  19. Colzi, L.; Rivilla, V.; Beltrán, M.; Jiménez-Serra, I.; Mininni, C.; Melosso, M.; Cesaroni, R.; Fontani, F.; Lorenzani, A.; Sánchez-Monge, A.; et al. The GUAPOS project-II. A comprehensive study of peptide-like bond molecules. Astron. Astrophys. 2021, 653, A129. [Google Scholar] [CrossRef]
  20. Belloche, A.; Meshcheryakov, A.; Garrod, R.; Ilyushin, V.; Alekseev, E.; Motiyenko, R.; Margulès, L.; Müller, H.; Menten, K. Rotational spectroscopy, tentative interstellar detection, and chemical modeling of N-methylformamide. Astron. Astrophys. 2017, 601, A49. [Google Scholar] [CrossRef]
  21. Belloche, A.; Garrod, R.; Müller, H.; Menten, K.; Medvedev, I.; Thomas, J.; Kisiel, Z. Re-exploring Molecular Complexity with ALMA (ReMoCA): Interstellar detection of urea. Astron. Astrophys. 2019, 628, A10. [Google Scholar] [CrossRef] [Green Version]
  22. Martín-Doménech, R.; Rivilla, V.; Jiménez-Serra, I.; Quénard, D.; Testi, L.; Martín-Pintado, J. Detection of methyl isocyanate (CH3NCO) in a solar-type protostar. Mon. Not. R. Astron. Soc. 2017, 469, 2230–2234. [Google Scholar] [CrossRef] [Green Version]
  23. Ligterink, N.F.; El-Abd, S.J.; Brogan, C.L.; Hunter, T.R.; Remijan, A.J.; Garrod, R.T.; McGuire, B.M. The family of amide molecules toward NGC 6334I. Astrophys. J. 2020, 901, 37. [Google Scholar] [CrossRef]
  24. Irvine, W.M.; Brown, R.D.; Cragg, D.M.; Friberg, P.; Godfrey, P.D.; Kaifu, N.; Matthews, H.E.; Ohishi, M.; Suzuki, H.; Takeo, H. A New Interstellar Polyatomic Molecule: Detection of Propynal in the Cold Cloud TMC-1. Astrophys. J. Lett. 1988, 335, L89. [Google Scholar] [CrossRef]
  25. Hollis, J.M.; Lovas, F.J.; Jewell, P.R. Interstellar Glycolaldehyde: The First Sugar. Astrophys. J. 2000, 540, L107. [Google Scholar] [CrossRef]
  26. Belloche, A.; Müller, H.S.P.; Menten, K.M.; Schilke, P.; Comito, C. Complex organic molecules in the interstellar medium: IRAM 30 m line survey of Sagittarius B2(N) and (M). Astron. Astrophys. 2013, 559, A47. [Google Scholar] [CrossRef] [Green Version]
  27. Manigand, S.; Coutens, A.; Loison, J.C.; Wakelam, V.; Calcutt, H.; Müller, H.S.P.; Jørgensen, J.K.; Taquet, V.; Wampfler, S.F.; Bourke, T.L.; et al. The ALMA-PILS survey: First detection of the unsaturated 3-carbon molecules Propenal (C2H3CHO) and Propylene (C3H6) towards IRAS 16293-2422 B. Astron. Astrophys. 2021, 645, A53. [Google Scholar] [CrossRef]
  28. Lattelais, M.; Pauzat, F.; Ellinger, Y.; Ceccarelli, C. Interstellar complex organic molecules and the minimum energy principle. Astrophys. J. 2009, 696, L133. [Google Scholar] [CrossRef]
  29. Loomis, R.A.; McGuire, B.A.; Shingledecker, C.; Johnson, C.H.; Blair, S.; Robertson, A.; Remijan, A.J. Investigating the minimum energy principle in searches for new molecular species—The case of H2C3O isomers. Astrophys. J. 2015, 799, 34. [Google Scholar] [CrossRef] [Green Version]
  30. Lattelais, M.; Pauzat, F.; Ellinger, Y.; Ceccarelli, C. A new weapon for the interstellar complex organic molecule hunt: The minimum energy principle. Astron. Astrophys. 2010, 519, A30. [Google Scholar] [CrossRef]
  31. Goerigk, L.; Grimme, S. A thorough benchmark of density functional methods for general main group thermochemistry, kinetics, and noncovalent interactions. Phys. Chem. Chem. Phys. 2011, 13, 6670–6688. [Google Scholar] [CrossRef]
  32. Melli, A.; Melosso, M.; Bizzocchi, L.; Alessandrini, S.; Jiang, N.; Tonolo, F.; Boi, S.; Castellan, G.; Sapienza, C.; Guillemin, J.C.; et al. Rotational Spectra of Unsaturated Carbon Chains Produced by Pyrolysis: The Case of Propadienone, Cyanovinylacetylene, and Allenylacetylene. J. Phys. Chem. A 2022, 126, 6210–6220. [Google Scholar] [CrossRef]
  33. Alessandrini, S.; Melosso, M.; Jiang, N.; Bizzocchi, L.; Dore, L.; Puzzarini, C. Conformational stability of cyclopropanecarboxaldehyde is ruled by vibrational effects. Mol. Phys. 2021, 119, e1955988. [Google Scholar] [CrossRef]
  34. Møllendal, H.; Margulés, L.; Motiyenko, R.A.; Larsen, N.W.; Guillemin, J.C. Rotational spectrum and conformational composition of cyanoacetaldehyde, a compound of potential prebiotic and astrochemical interest. J. Phys. Chem. A 2012, 116, 4047–4056. [Google Scholar] [CrossRef]
  35. Behnke, M.; Medvedev, I.; Winnewisser, M.; De Lucia, F.C.; Herbst, E. The Millimeter-and Submillimeter-Wave Spectrum of Oxiranecarbonitrile. Astrophys. J. Suppl. S. 2004, 152, 97. [Google Scholar] [CrossRef] [Green Version]
  36. Mackrodt, W.; Wardley, A.; Curnuck, P.; Owen, N.; Sheridan, J. Microwave spectra of oxazole and isoxazole. Chem. Comm. 1966, 692–693. [Google Scholar] [CrossRef]
  37. Sauer, E.G.; Oja, T. Nitrogen-14 nuclear quadrupole resonance in isoxazole. J. Chem. Phys. 1973, 58, 2650–2651. [Google Scholar] [CrossRef]
  38. Stiefvater, O.; Nösberger, P.; Sheridan, J. Microwave spectrum and structure of isoxazole. Chem. Phys. 1975, 9, 435–444. [Google Scholar] [CrossRef]
  39. Kumar, A.; Sheridan, J.; Stiefvater, O.L. The microwave spectrum of oxazole I. The complete structure by DRM microwave spectroscopy. Z. Naturforsch. A 1978, 33, 145–152. [Google Scholar] [CrossRef]
  40. Kumar, A.; Sheridan, J.; Stiefvater, O.L. The Microwave Spectrum of Oxazole: II. Dipole Moment and Quadrupole Coupling Constants. Z. Naturforsch. A 1978, 33, 549–558. [Google Scholar] [CrossRef] [Green Version]
  41. Kraśnicki, A.; Pszczółkowski, L.; Kisiel, Z. Analysis of the rotational spectrum of pyruvonitrile up to 324 GHz. J. Mol. Spectrosc. 2010, 260, 57–65. [Google Scholar] [CrossRef]
  42. Bouchy, A.; Roussy, G. Geometrical and electronic structures of two conformers of vinyl isocyanate by microwave spectroscopy. J. Mol. Spectrosc. 1977, 68, 156–165. [Google Scholar] [CrossRef]
  43. Kirby, C.; Kroto, H. Microwave and photoelectron study of cis-and trans-isocyanato ethene, CH2=CHNCO (vinylisocyanate). J. Mol. Spectrosc. 1978, 70, 216–228. [Google Scholar] [CrossRef]
  44. Bouchy, A.; Roussy, G. The determination of the dipole moment in microwave spectroscopy: The case of cis vinyl isocyanate. J. Mol. Spectrosc. 1979, 77, 145–146. [Google Scholar] [CrossRef]
  45. Vávra, K.; Kolesniková, L.; Belloche, A.; Garrod, R.; Kouckỳ, J.; Uhlíková, T.; Luková, K.; Guillemin, J.C.; Kania, P.; Müller, H.; et al. Millimeter wave spectrum and search for vinyl isocyanate toward Sgr B2 (N)with ALMA. Astron. Astrophys. 2022, 666, A50. [Google Scholar] [CrossRef]
  46. Vijayakumar, S.; Kolandaivel, P. Isomerization of C3H3NO isomers: Ab initio study. Mol. Phys. 2006, 104, 1401–1411. [Google Scholar] [CrossRef]
  47. Chong, D.P.; Hu, C.H. Accurate density functional calculation of core-electron binding energies with a scaled polarized triple-zeta basis set. IV. Application to isomers of C3H6O, C3H3NO, and C6H6. J. Chem. Phys. 1998, 108, 8950–8956. [Google Scholar] [CrossRef]
  48. Santra, G.; Sylvetsky, N.; Martin, J.M.L. Minimally Empirical Double-Hybrid Functionals Trained against the GMTKN55 Database: RevDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4. J. Phys. Chem. A 2019, 123, 5129–5143. [Google Scholar] [CrossRef]
  49. Dunning, T.H., Jr. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  50. Papajak, E.; Truhlar, D.G. Convergent Partially Augmented Basis Sets for Post-Hartree-Fock Calculations of Molecular Properties and Reaction Barrier Heights. J. Chem. Theory Comput. 2011, 7, 10–18. [Google Scholar] [CrossRef]
  51. 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. [Google Scholar] [CrossRef]
  52. Alessandrini, S.; Melosso, M. Fate of the Gas-Phase Reaction Between Oxirane and the CN Radical in Interstellar Conditions. Front. Astron. Space Sci. 2021, 8, 754977. [Google Scholar] [CrossRef]
  53. Turner, B.E.; Apponi, A.J. Microwave detection of interstellar vinyl alcohol, CH2CHOH. Astrophys. J. 2001, 561, L207. [Google Scholar] [CrossRef]
  54. Agúndez, M.; Marcelino, N.; Tercero, B.; Cabezas, C.; de Vicente, P.; Cernicharo, J. O-bearing complex organic molecules at the cyanopolyyne peak of TMC-1: Detection of C2H3CHO, C2H3OH, HCOOCH3, and CH3OCH3. Astron. Astrophys. 2021, 649, L4. [Google Scholar] [CrossRef] [PubMed]
  55. Tonolo, F.; Lupi, J.; Puzzarini, C.; Barone, V. The quest for a plausible formation route of formyl cyanide in the interstellar medium: A state-of-the-art quantum-chemical and kinetic approach. Astrophys. J. 2020, 900, 85. [Google Scholar] [CrossRef]
  56. Barone, V.; Biczysko, M.; Bloino, J.; Puzzarini, C. Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: The challenge of the conformational equilibrium in glycine. Phys. Chem. Chem. Phys. 2013, 15, 10094–10111. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Puzzarini, C.; Bloino, J.; Tasinato, N.; Barone, V. Accuracy and Interpretability: The Devil and the Holy Grail. New Routes across Old Boundaries in Computational Spectroscopy. Chem. Rev. 2019, 119, 8131–8191. [Google Scholar] [CrossRef]
  58. Puzzarini, C.; Barone, V. The challenging playground of astrochemistry: An integrated rotational spectroscopy–quantum chemistry strategy. Phys. Chem. Chem. Phys. 2020, 22, 6507–6523. [Google Scholar] [CrossRef]
  59. Lupi, J.; Puzzarini, C.; Cavallotti, C.; Barone, V. State-of-the-Art Quantum Chemistry Meets Variable Reaction Coordinate Transition State Theory to Solve the Puzzling Case of the H2S + Cl System. J. Chem. Theory Comput. 2020, 16, 5090–5104. [Google Scholar] [CrossRef]
  60. Little, G.B.; Gerry, M. The microwave spectrum, structure, and dipole moment of propiolamide, HCCC(O)NH2. J. Mol. Spectrosc. 1978, 71, 321–342. [Google Scholar] [CrossRef]
  61. Hegelund, F.; Larsen, R.W.; Palmer, M. High-resolution infrared and theoretical study of gaseous oxazole in the 600–1400 cm−1 region. J. Mol. Spectr. 2007, 241, 26–44. [Google Scholar] [CrossRef]
  62. Puzzarini, C.; Stanton, J.F.; Gauss, J. Quantum-chemical calculation of spectroscopic parameters for rotational spectroscopy. Int. Rev. Phys. Chem. 2010, 29, 273–367. [Google Scholar] [CrossRef]
  63. Barone, V.; Alessandrini, S.; Biczysko, M.; Cheeseman, 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. [Google Scholar] [CrossRef]
  64. Heise, H.; Scappini, F.; Dreizler, H. Acetyl Cyanide III. Vibrational Spectrum and Vibrational Analysis. Z. Naturforsch. A 1976, 31, 1408–1412. [Google Scholar] [CrossRef]
  65. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian16 Revision C.01, 2016; Gaussian Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  66. 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. [Google Scholar] [CrossRef] [PubMed]
  67. Stanton, J.F.; Gauss, J.; Cheng, L.; Harding, M.E.; Matthews, D.A.; Szalay, P.G. CFOUR, Coupled-Cluster Techniques for Computational Chemistry, a Quantum-Chemical Program Package. Available online: http://www.cfour.de (accessed on 1 March 2023).
  68. Kállay, M.; Nagy, P.R.; Mester, D.; Rolik, Z.; Samu, G.; Csontos, J.; Csóka, J.; Szabó, P.B.; Gyevi-Nagy, L.; Hégely, B.; et al. The MRCC program system: Accurate quantum chemistry from water to proteins. J. Chem. Phys. 2020, 152, 074107. [Google Scholar] [CrossRef] [PubMed]
  69. Shavitt, I.; Bartlett, R.J. Many-Body Methods in Chemistry and Physics: MBPT and Coupled-Cluster Theory; Cambridge University Press: Cambridge, UK, 2009. [Google Scholar]
  70. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456–1465. [Google Scholar] [CrossRef]
  72. Mancini, G.; Fusè, M.; Lazzari, F.; Chandramouli, B.; Barone, V. Unsupervised search of low-lying conformers with spectroscopic accuracy: A two-step algorithm rooted into the island model evolutionary algorithm. J. Chem. Phys. 2020, 153, 124110. [Google Scholar] [CrossRef]
  73. León, I.; Fusè, M.; Alonso, E.R.; Mata, S.; Mancini, G.; Puzzarini, C.; Alonso, J.l.; Barone, V. Unbiased disentanglement of conformational baths with the help of microwave spectroscopy, quantum chemistry, and artificial intelligence: The puzzling case of homocysteine. J. Chem. Phys. 2022, 157, 074107. [Google Scholar] [CrossRef]
  74. Barone, V.; Fusè, M.; Aguado, R.; Potenti, S.; Leon, I.; Alonso, E.R.; Mata, S.; Lazzari, F.; Mancini, G.; Spada, L.; et al. Bringing Machine-Learning Enhanced Quantum Chemistry and Microwave Spectroscopy to Conformational Landscape Exploration: The Paradigmatic Case of 4-Fluoro-Threonine. Chem. Eur. J. 2023; in press. [Google Scholar] [CrossRef]
  75. Puzzarini, C.; Gambi, A. The energetics and structural properties of diazomethyl (HCNN) and cyanomidyl (HNCN) radicals and their related cations and anions from ab initio calculations. J. Chem. Phys. 2005, 122, 064316. [Google Scholar] [CrossRef] [PubMed]
  76. Feller, D.; Peterson, K.A.; Dixon, D.A. A survey of factors contributing to accurate theoretical predictions of atomization energies and molecular structures. J. Chem. Phys. 2008, 129, 204105. [Google Scholar] [CrossRef]
  77. Puzzarini, C. Accurate thermochemistry and spectroscopy of the oxygen-protonated sulfur dioxide isomers. Phys. Chem. Chem. Phys. 2011, 13, 21319–21327. [Google Scholar] [CrossRef]
  78. Alessandrini, S.; Puzzarini, C. Structural and energetic characterization of prebiotic molecules: The case study of formamide and its dimer. J. Phys. Chem. A 2016, 120, 5257–5263. [Google Scholar] [CrossRef] [PubMed]
  79. Feller, D. The Use of Systematic Sequences of Wave Functions for Estimating the Complete Basis Set, Full Configuration Interaction Limit in Water. J. Phys. Chem. 1993, 98, 7059–7071. [Google Scholar] [CrossRef]
  80. Helgaker, T.; Klopper, W.; Koch, H.; Noga, J. Basis-Set Convergence of Correlated Calculations on Water. J. Chem. Phys. 1997, 106, 9639. [Google Scholar] [CrossRef]
  81. Woon, D.E.; Dunning, T.H., Jr. Gaussian Basis Sets for Use in Correlated Molecular Calculations. V. Core-Valence Basis Sets for Boron through Neon. J. Chem. Phys. 1995, 103, 4572. [Google Scholar] [CrossRef] [Green Version]
  82. Heckert, M.; Kállay, M.; Gauss, J. Molecular Equilibrium Geometries Based on Coupled-Cluster Calculations Including Quadruple Excitations. Mol. Phys. 2005, 103, 2109–2115. [Google Scholar] [CrossRef]
  83. 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. [Google Scholar] [CrossRef] [PubMed]
  84. Mills, I.M. Vibration–Rotation Structure in Asymmetric- and Symmetric-Top Molecules. In Molecular Spectroscopy: Modern Research; Rao, K.N., Matthews, C.W., Eds.; Academic Press: Cambridge, MA, USA, 1972; Volume 1, pp. 115–140. [Google Scholar]
  85. Puzzarini, C.; Stanton, J.F. Connections between the accuracy of rotational constants and equilibrium molecular structures. Phys. Chem. Chem. Phys. 2023, 25, 1421–1429. [Google Scholar] [CrossRef] [PubMed]
  86. Piccardo, M.; Bloino, J.; Barone, V. Generalized vibrational perturbation theory for rotovibrational energies of linear, symmetric and asymmetric tops: Theory, approximations, and automated approaches to deal with medium-to-large molecular systems. Int. J. Quantum Chem. 2015, 115, 948–982. [Google Scholar] [CrossRef] [Green Version]
  87. Noga, J.; Bartlett, R.J. The full CCSDT model for molecular electronic structure. J. Chem. Phys. 1987, 86, 7041–7050. [Google Scholar] [CrossRef]
  88. Scuseria, G.E.; Schaefer, H.F., III. A new implementation of the full CCSDT model for molecular electronic structure. Chem. Phys. Lett. 1988, 152, 382–386. [Google Scholar] [CrossRef]
  89. Bomble, Y.J.; Stanton, J.F.; Kállay, M.; Gauss, J. Coupled-cluster methods including noniterative corrections for quadruple excitations. J. Chem. Phys. 2005, 123, 054101. [Google Scholar] [CrossRef]
  90. Kállay, M.; Gauss, J. Approximate treatment of higher excitations in coupled-cluster theory. J. Chem. Phys. 2005, 123, 214105. [Google Scholar] [CrossRef]
  91. Kállay, M.; Gauss, J. Approximate treatment of higher excitations in coupled-cluster theory. II. Extension to general single-determinant reference functions and improved approaches for the canonical Hartree–Fock case. J. Chem. Phys. 2008, 129, 144101. [Google Scholar] [CrossRef]
  92. Handy, N.C.; Yamaguchi, Y.; Schaefer, H.F., III. The diagonal correction to the Born–Oppenheimer approximation: Its effect on the singlet–triplet splitting of CH2 and other molecular effects. J. Chem. Phys. 1986, 84, 4481–4484. [Google Scholar] [CrossRef]
  93. Klopper, W. Simple recipe for implementing computation of first-order relativistic corrections to electron correlation energies in framework of direct perturbation theory. J. Comp. Chem. 1997, 18, 20–27. [Google Scholar] [CrossRef]
  94. Cowan, R.D.; Griffin, D.C. Approximate relativistic corrections to atomic radial wave functions. J. Opt. Soc. Am. 1976, 66, 1010–1014. [Google Scholar] [CrossRef]
  95. Rodríguez-Almeida, L.; Rivilla, V.; Jiménez-Serra, I.; Melosso, M.; Colzi, L.; Zeng, S.; Tercero, B.; de Vicente, P.; Martín, S.; Requena-Torres, M.; et al. First Detection of C2H5NCO in the ISM and Search of Other Isocyanates towards the G + 0.693-0.027 Molecular Cloud. Astron. Astrophys. 2021, 654, L1. [Google Scholar] [CrossRef]
  96. Cernicharo, J.; Agúndez, M.; Cabezas, C.; Marcelino, N.; Tercero, B.; Pardo, J.R.; Fuentetaja, R.; de Vicente, P. The QUIJOTE line survey of TMC-1. In Proceedings of the EPJ Web of Conferences, EDP Sciences, Sacramento, CA, USA, 24–29 July 2022; Volume 265. [Google Scholar]
  97. McGuire, B.A.; Burkhardt, A.M.; Loomis, R.A.; Shingledecker, C.N.; Lee, K.L.K.; Charnley, S.B.; Cordiner, M.A.; Herbst, E.; Kalenskii, S.; Momjian, E.; et al. Early science from GOTHAM: Project overview, methods, and the detection of interstellar propargyl cyanide (HCCCH2CN) in TMC-1. Astrophys. J. Lett. 2020, 900, L10. [Google Scholar] [CrossRef]
Figure 1. Relative energies of the 67 isomers of the C3H3NO family as obtained from the revDSD/junTZ level of theory, including the hZPE contribution. The reader may refer to the SM for the correspondence between the labeling and the molecular structure.
Figure 1. Relative energies of the 67 isomers of the C3H3NO family as obtained from the revDSD/junTZ level of theory, including the hZPE contribution. The reader may refer to the SM for the correspondence between the labeling and the molecular structure.
Molecules 28 03226 g001
Figure 2. Schematic representation of the 20 isomers considered in the second step of our procedure (energy refinement). The relative energies are at the CBS + CV + hZPE level and computed on top of the revDSD/junTZ geometries. The 10 most stable isomers are reported in blue, while light blue is employed for the remaining species.
Figure 2. Schematic representation of the 20 isomers considered in the second step of our procedure (energy refinement). The relative energies are at the CBS + CV + hZPE level and computed on top of the revDSD/junTZ geometries. The 10 most stable isomers are reported in blue, while light blue is employed for the remaining species.
Molecules 28 03226 g002
Figure 3. Total electric dipole moment (revDSD/junTZ) plotted against relative energy (CBS + CV + hZPE) for the 20 most stable isomers of the C3H3NO family. For labels, see Figure 2 and the text.
Figure 3. Total electric dipole moment (revDSD/junTZ) plotted against relative energy (CBS + CV + hZPE) for the 20 most stable isomers of the C3H3NO family. For labels, see Figure 2 and the text.
Molecules 28 03226 g003
Figure 4. Simulated rotational spectra of the four isomers of cyanovinylalcohol based on the computed constants reported in Table 2. The calculations were performed for J ranging in the 0–100 interval and T = 300 K. The y-axes are normalized with respect to the most intense Z t spectrum. The relative energy at the CBS + CV + aZPE level of theory is provided with respect to the most stable isomer ( Z g ).
Figure 4. Simulated rotational spectra of the four isomers of cyanovinylalcohol based on the computed constants reported in Table 2. The calculations were performed for J ranging in the 0–100 interval and T = 300 K. The y-axes are normalized with respect to the most intense Z t spectrum. The relative energy at the CBS + CV + aZPE level of theory is provided with respect to the most stable isomer ( Z g ).
Molecules 28 03226 g004
Figure 5. Schematic representation of the computational protocol developed in this work for the individuation of potential candidates for astronomical observations.
Figure 5. Schematic representation of the computational protocol developed in this work for the individuation of potential candidates for astronomical observations.
Molecules 28 03226 g005
Table 1. Experimental a (Exp.) and computed b (Theory) rotational constants of the vibrational ground state (MHz) for vinylisocyanate (t and g), oxazole, cyanoacetaldehyde (t and g) and pyruvonitrile.
Table 1. Experimental a (Exp.) and computed b (Theory) rotational constants of the vibrational ground state (MHz) for vinylisocyanate (t and g), oxazole, cyanoacetaldehyde (t and g) and pyruvonitrile.
CH2CHNCOc-N(CH)2OCH
transgauche
TheoryExp.TheoryExp.TheoryExp.
A 0 62,759.3 (0.28%)62,586.320,306.5 (0.81%)20,144.110,058.0 (0.07%)10,051.0
B 0 2437.2 (0.02%)2437.83097.7 (0.31%)3107.59659.1 (0.14%)9645.8
C 0 2346.2 (0.01%)2346.52684.7 (0.17%)2689.44924.1 (0.10%)4919.4
NCCH2CHOCH3COCN
transgauche
TheoryExp.TheoryExp.TheoryExp.
A 0 27,147.3 (1.46%)267,56.812,827.5 (0.12%)12,812.110,194.2 (0.08%)10,185.9
B 0 2569.5 (0.25%)2576.03657.5 (0.03%)3658.74157.6 (0.004%)4157.7
C 0 2416.9 (0.28%)2423.62894.5 (0.002%)2894.63003.5 (0.03%)3002.8
a See text for experimental references. b CBS + CV equilibrium rotational constants corrected for the vibrational contribution at the revDSD/junTZ level. Relative deviations from the experiment are given in parentheses.
Table 2. Computed rotational parameters for the (E)- and (Z)-cyanovinylalcohol forms and their relative energies.
Table 2. Computed rotational parameters for the (E)- and (Z)-cyanovinylalcohol forms and their relative energies.
Parameter aUnits(Z)-Cyanovinylalcohol(E)-Cyanovinylalcohol
transgauchetransgauche
A 0 MHz13,893.012,869.245,505.944,451.1
B 0 MHz3554.53845.52387.02384.3
C 0 MHz2827.22956.32267.22262.0
D J kHz3.14.60.330.33
D J K kHz−27.3−32.6−35.7−32.54
D K kHz117.1100.32330.01981.70
d 1 kHz−0.95−1.5−0.041−0.040
d 2 Hz−60.6−0.1−0.92−0.96
H J Hz0.0110.0250.000330.00032
H J K Hz−0.072−0.18−0.072−0.067
H K J Hz−0.52−0.282.22.2
H K Hz3.32.469.832.2
h 1 Hz0.00540.0120.00010.0001
h 2 mHz0.881.80.00310.0031
h 3 mHz0.160.350.00130.0013
χ a a MHz−2.1−1.6−1.3−3.6
χ b b χ c c MHz−1.8−1.670.10.4
χ a b MHz3.03.30.71.85
| μ a | D5.452.604.784.03
| μ b | D2.161.830.931.72
Δ E kJ mol−114.10.011.49.8
a Rotational constants of the vibrational ground state obtained from the CBS + CV equilibrium terms augmented for the vibrational correction at the revDSD/junTZ level. Centrifugal distortion and nuclear quadrupole coupling constants as well as vibrationally averaged dipole moment components from the revDSD/junTZ anharmonic force field calculations. Watson’s Hamiltonian, S-reduction, I I I r representation. Relative energies at the CBS + CV + aZPE level of theory.
Table 3. Computed vibrational frequencies of oxazole and t-vinylisocyanate compared to the available experimental values. All values are in cm−1.
Table 3. Computed vibrational frequencies of oxazole and t-vinylisocyanate compared to the available experimental values. All values are in cm−1.
oxazole
CC/TZ a Δ ν anh bbest EstimateExp. c
ν 1 3309.09−126.553182.54
ν 2 3285.40−126.813158.59
ν 3 3275.62−125.383150.24
ν 4 1575.76−35.971539.79
ν 5 1540.35−43.161497.19
ν 6 1358.40−32.971325.431329.75165(3)
ν 7 1280.23−27.821252.41
ν 8 1177.36−34.361143.001142.50528(3)
ν 9 1112.64−26.011086.631091.12069(5)
ν 10 1108.71−28.001080.701081.29060(4)
ν 11 1074.50−23.581050.921051.75844(3)
ν 12 919.16−13.30905.86909.28465(6)
ν 13 908.75−13.73895.01899.33009(6)
ν 14 872.69−18.31854.38859.19(1)
ν 15 851.29−18.58832.71832.01870(3)
ν 16 763.10−17.15745.95749.31060(3)
ν 17 655.72−11.31644.41646.35537(3)
ν 18 618.04−10.21607.83
t-vinylisocyanate
CC/TZ a Δ ν anh bbest estimateLit. d
ν 1 3270.82−137.943132.873149.9
ν 2 3186.24−132.333053.923079.0
ν 3 3168.87−106.843062.033051.5
ν 4 2324.44−60.722263.712301.7
ν 5 1689.63−53.051636.581635.6
ν 6 1484.92−28.651456.271490.3
ν 7 1405.70−33.451372.251408.3
ν 8 1335.24−26.771308.481328.3
ν 9 1115.46−22.701092.771109.3
ν 10 972.93−14.57958.36972.1
ν 11 899.79−3.96895.84915.6
ν 12 854.820.76855.58853.5
ν 13 694.681.63696.31698.4
ν 14 639.59−24.99614.60656.0
ν 15 573.32−17.43555.89601.1
ν 16 446.77−11.57435.21453.1
ν 17 143.04−2.45140.58136.5
ν 18 83.223.4586.6778.3
a Harmonic frequencies at the fc-CCSD(T)/cc-pVTZ level of theory. b Anharmonic corrections at the revDSD/junTZ level. c Experimental data taken from Reference [61]. d The literature values computed at the CCSD/cc-pVTZ level (see Reference [45]).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alessandrini, S.; Melosso, M.; Rivilla, V.M.; Bizzocchi, L.; Puzzarini, C. Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers. Molecules 2023, 28, 3226. https://doi.org/10.3390/molecules28073226

AMA Style

Alessandrini S, Melosso M, Rivilla VM, Bizzocchi L, Puzzarini C. Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers. Molecules. 2023; 28(7):3226. https://doi.org/10.3390/molecules28073226

Chicago/Turabian Style

Alessandrini, Silvia, Mattia Melosso, Víctor M. Rivilla, Luca Bizzocchi, and Cristina Puzzarini. 2023. "Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers" Molecules 28, no. 7: 3226. https://doi.org/10.3390/molecules28073226

APA Style

Alessandrini, S., Melosso, M., Rivilla, V. M., Bizzocchi, L., & Puzzarini, C. (2023). Computational Protocol for the Identification of Candidates for Radioastronomical Detection and Its Application to the C3H3NO Family of Isomers. Molecules, 28(7), 3226. https://doi.org/10.3390/molecules28073226

Article Metrics

Back to TopTop