Carbamic acid and its dimer: A computational study.

A recent work by Marks et al. on the formation of carbamic acid in NH 3 -CO 2 interstellar ices pointed out its stability in the gas phase and the concomitant production of its dimer. Prompted by these results and the lack of information on these species, we have performed an accurate structural, energetic and spectroscopic investigation of carbamic acid and its dimer. For the former, the structural and spectroscopic characterization employed composite schemes based on coupled cluster (CC) calculations that account for the extrapolation to the complete basis set limit and core correlation effects. A first important outcome is the definitive confirmation of the nonplanarity of carbamic acid, then followed by an accurate estimate of its rotational and vibrational spectroscopy parameters. As far as the carbamic acid dimer is concerned, the investigation started from the identification of its most stable forms. For them, structure and vibrational properties have been evaluated using density functional theory, while a composite scheme rooted in CC theory has been employed for the energetic characterization. Our results allowed us to provide a better interpretation of the feature observed in the recent experiment mentioned above.

Isomerism of CH 2 SO : Accurate structural, energetic, and spectroscopic characterization.

A recent work [Ye et al. Mon. Not. R. Astron. Soc. 2023, 525, 1158] on the gas-phase formation of t-HC(O)SH, already detected in the interstellar medium, pointed out that the trans form of HC(S)OH is a potential candidate for astronomical observations. Prompted by these results, the CH 2 SO family of isomers has been investigated from an energetic point of view using a double-hybrid density functional in combination with a partially augmented triple-zeta basis set. This preliminary study showed that the most stable species of the family are the cis and trans forms of HC(O)SH and HC(S)OH. For their structural and spectroscopic characterization, a composite scheme based on coupled cluster (CC) calculations that incorporates up to the quadruple excitations and accounts for the extrapolation to the complete basis set limit and core correlation effects has been employed. This approach opens to the prediction of rotational constants with an accuracy of 0.1%. A hybrid scheme, based on harmonic frequencies computed using the CC singles, doubles and a perturbative treatment of triples method (CCSD(T)) in conjunction with a quadruple-zeta basis set, allowed us to obtain fundamental vibrational frequencies with a mean absolute error of about 1%.

Gas-Phase Computational Spectroscopy: The Challenge of the Molecular Bricks of Life.

Gas-phase molecular spectroscopy is a natural playground for accurate quantum-chemical computations. However, the molecular bricks of life (e.g., DNA bases or amino acids) are challenging systems because of the unfavorable scaling of quantum-chemical models with the molecular size (active electrons) and/or the presence of large-amplitude internal motions. From the theoretical point of view, both aspects prevent the brute-force use of very accurate but very expensive state-of-the-art quantum-chemical methodologies. From the experimental point of view, both features lead to congested gas-phase spectra, whose assignment and interpretation are not at all straightforward. Based on these premises, this review focuses on the current status and perspectives of the fully a priori prediction of the spectral signatures of medium-sized molecules (containing up to two dozen atoms) in the gas phase with special reference to rotational and vibrational spectroscopies of some representative molecular bricks of life.

Describing the Disulfide Bond: From the Density Functional Theory and Back through the "Lego Brick" Approach.

Selected molecular species containing the disulfide bond, RSSR, have been considered, these ranging from hydrogen disulfide, H2S2 (R = H), to diphenyl disulfide with R = C6H5. The aim of this work is two-fold: (i) to investigate different computational approaches in order to derive accurate equilibrium structures at an affordable cost, (ii) to employ the results from the first goal in order to benchmark cheaper methodologies rooted in the density functional theory. Among the strategies used for the accurate geometrical determinations, the semiexperimental approach has been exploited in combination with a reduced-dimensionality VPT2 model, without however obtaining satisfactory results. Instead, the so-called "Lego brick" approach turned out to be very effective despite the flexibility of the systems investigated. Concerning the second target of this work, the focus was mainly on the S-S bond and the structural parameters related to it. Among those tested, PBE0(-D3BJ), M06-2X(-D3) and DSD-PBEP86-D3BJ have been found to be the best-performing functionals.

The Semiexperimental Approach at Work: Equilibrium Structure of Radical Species.

The so-called semiexperimental (SE) approach is a powerful technique for obtaining highly accurate equilibrium structures for isolated systems. This Featured Article describes its extension to open-shell species, thus providing the first systematic investigation on radical equilibrium geometries to be used for benchmarking purposes. The small yet significant database obtained demonstrates that there is no reduction in accuracy when moving from closed-shell species to radicals. We also provide an extension of the applicability of the SE approach to medium-/large-sized radicals by exploiting the so-called "Lego-brick" approach, which is based on the assumption that a molecular system can be seen as formed by smaller fragments for which the SE equilibrium structure is available. In this Featured Article we show that this model can be successfully applied also to open-shell species.

Tracing Prebiotic Molecules: Rotational Spectroscopy of Deuterated Glycolaldehyde and (Z)-1,2-Ethenediol.

Glycolaldehyde, an important prebiotic molecule, along with its monodeuterated species and its higher energy tautomer, (Z)-1,2-ethenediol, has been detected in the interstellar medium. Although the elemental D/H ratio in the universe is only ∼1.6 × 10-5, the deuterium relative abundance in interstellar molecules might be by far larger than this. As such, it provides a remarkable and almost unique diagnostic tool. In particular, it might help elucidate the reaction mechanisms that lead to the formation of the so-called complex organic molecules. It is therefore crucial to extend the census of the interstellar deuterated molecules. To this aim, in this work, we present for the first time a spectroscopic investigation of the rotational spectra of the CHDOD-CHO bideuterated variant of glycolaldehyde and of mono- and bideuterated species of (Z)-1,2-ethenediol (CHOD═CHOD, CHOD═CHOH, and CHOH═CHOD rotamers). For each species, more than a hundred transitions have been assigned. Their analysis led to the accurate determination of all rotational constants as well as quartic and sextic centrifugal distortion terms, thus providing spectroscopic line catalogs suitable for supporting astronomical searches. In addition, the rotational constants of the bideuterated glycolaldehyde isotopologue studied in this work allowed us to improve the semiexperimental equilibrium structure determination for this molecule.

Hunting for interstellar molecules: rotational spectra of reactive species.

Interstellar molecules are often highly reactive species, which are unstable under terrestrial conditions, such as radicals, ions and unsaturated carbon chains. Their detection in space is usually based on the astronomical observation of their rotational fingerprints. However, laboratory investigations have to face the issue of efficiently producing these molecules and preserving them during rotational spectroscopy measurements. A general approach for producing and investigating unstable/reactive species is presented by means of selected case-study molecules. The overall strategy starts from quantum-chemical calculations that aim at obtaining accurate predictions of the missing spectroscopic information required to guide spectral analysis and assignment. Rotational spectra of these species are then recorded by exploiting the approach mentioned above, and their subsequent analysis leads to accurate spectroscopic parameters. These are then used for setting up accurate line catalogs for astronomical searches.

Connections between the accuracy of rotational constants and equilibrium molecular structures.

Rotational spectroscopy is the technique of choice for investigating molecular structures in the gas phase. Indeed, rotational constants are strongly connected to the geometry of the molecular system under consideration. Therefore, they are powerful tools for assessing the accuracy that quantum chemical approaches can reach in structural determinations. In this review article, it is shown how it is possible to measure the accuracy of a computed equilibrium geometry based on the comparison of rotational constants. But, it is also addressed what accuracy is required by computations for providing molecular structures and thus rotational constants that are useful to experiment. Quantum chemical methodologies for obtaining the "0.1% accuracy" for rotational constants are reviewed for systems ranging in size from small molecules to small polycyclic aromatic hydrocarbons. This accuracy for systems containing two dozen or so atoms opens the way towards future applications such as the accurate characterization of non-covalent interactions, which play a key role in several biological and technological processes.

Insights into the molecular structure and infrared spectrum of the prebiotic species aminoacetonitrile.

Aminoacetonitrile is an interstellar molecule with a prominent prebiotic role, already detected in the chemically-rich molecular cloud Sagittarius B2(N) and postulated to be present in the atmosphere of the largest Saturn's moon, Titan. To further support its observation in such remote environments and laboratory experiments aimed at improving our understanding of interstellar chemistry, we report a thorough spectroscopic and structural characterization of aminoacetonitrile. Equilibrium geometry, fundamental bands as well as spectroscopic and molecular parameters have been accurately computed by exploiting a composite scheme rooted in the coupled-cluster theory that accounts for the extrapolation to the complete basis set limit and core-correlation effects. In addition, a semi-experimental approach that combines ground-state rotational constants for different isotopic species and calculated vibrational corrections has been employed for the structure determination. From the experimental side, we report the analysis of the three strongest fundamental bands of aminoacetonitrile observed between 500 and 1000 cm-1 in high-resolution infrared spectra. More generally, all computed band positions are in excellent agreement with the present and previous experiments. The only exception is the ν15 band, for which we provide a revision of the experimental assignment, now in good agreement with theory.

A rotational investigation of the three isomeric forms of cyanoethynylbenzene (HCC-C6H4-CN): benchmarking experiments and calculations using the "Lego brick" approach.

We report the study of three structural isomers of phenylpropiolonitrile (3-phenyl-2-propynenitrile, C6H5-C3N) containing an alkyne function and a cyano group, namely ortho-, meta-, and para-cyanoethynylbenzene (HCC-C6H4-CN). The pure rotational spectra of these species have been recorded at room temperature in the millimeter-wave domain using a chirped-pulse spectrometer (75-110 GHz) and a source-frequency modulation spectrometer (140-220 GHz). Assignments of transitions in the vibrational ground state and several vibrationally excited states were supported by quantum chemical calculations using the so-called "Lego brick" approach [A. Melli, F. Tonolo, V. Barone and C. Puzzarini, J. Phys. Chem. A, 2021, 125, 9904-9916]. From these assignments, accurate spectroscopic (rotational and centrifugal distortion) constants have been derived: for all species and all observed vibrational states, predicted rotational constants show relative accuracy better than 0.1%, and often of the order of 0.01%, compared to the experimental values. The present work hence further validates the use of the "Lego brick" approach for predicting spectroscopic constants with high precision.

Rotational spectra and semi-experimental structures of furonitrile and its water cluster.

Rotational spectroscopy represents an invaluable tool for several applications: from the identification of new molecules in interstellar objects to the characterization of van der Waals complexes, but also for the determination of very accurate molecular structures and for conformational analyses. In this work, we used high-resolution rotational spectroscopic techniques in combination with high-level quantum-chemical calculations to address all these aspects for two isomers of cyanofuran, namely 2-furonitrile and 3-furonitrile. In particular, we have recorded and analyzed the rotational spectra of both of them from 6 to 320 GHz; rotational transitions belonging to several singly-substituted isotopologues have been identified as well. The rotational constants derived in this way have been used in conjunction with computed rotation-vibration interaction constants in order to derive a semi-experimental equilibrium structure for both isomers. Moreover, we observed the rotational spectra of four different intermolecular adducts formed by furonitrile and water, whose identification has been supported by a conformational analysis and a theoretical spectroscopic characterization. A semi-experimental determination of the intermolecular parameters has been achieved for all of them and the results have been compared with those obtained for the analogous system formed by benzonitrile and water.

Insights into the Molecular Structure and Spectroscopic Properties of HONCO: An Accurate Ab Initio Study.

In an effort to provide the first accurate structural and spectroscopic characterization of the quasi-linear chain HONCO in its electronic ground state, state-of-the-art computational approaches mainly based on coupled-cluster (CC) theory have been employed. Equilibrium geometries have been calculated by means of a composite scheme based on CC calculations that incorporates up to the quadruple excitations and accounts for the extrapolation to the complete basis set limit and core correlation effects. This approach is proven to provide molecular structures with an accuracy better than 0.001 Å and 0.05° for bond lengths and angles, respectively. Incorporation of vibrational effects permits this level of theory to predict rotational constants with an estimated accuracy of 0.1% or better. Vibrational fundamental bands have been evaluated by means of a hybrid scheme based on harmonic frequencies computed using the CC singles, doubles, and a perturbative treatment of the triples method (CCSD(T)) in conjunction with a quadruple-ζ basis set, with all electrons being correlated, and anharmonic corrections from CCSD(T) calculations using a triple-ζ basis set, within the frozen-core approximation. Such a hybrid approach allowed us to obtain fundamental frequencies with a mean absolute error of about 1%. To complete the spectroscopic characterization, vertical electronic excitation energies have been calculated for the lowest singlet and triplet states using the internally contracted multireference configuration interaction (MRCI) method. Computations show that HONCO dissociates into OH + NCO upon the absorption of UV-vis light. In conclusion, we are confident that the highly accurate spectroscopic data provided herein can be useful for guiding future experimental investigations and supporting the characterization of this molecule in atmospheric and astrophysical media, as well as in combustion.

Hyperfine-resolved spectra of HDS together with a global ro-vibrational analysis.

Despite their chemical simplicity, the spectroscopic investigation of light hydrides, such as hydrogen sulfide, is challenging due to strong hyperfine interactions and/or anomalous centrifugal-distortion effects. Several hydrides have already been detected in the interstellar medium, and the list includes H2S and some of its isotopologues. Astronomical observation of isotopic species and, in particular, those bearing deuterium is important to gain insights into the evolutionary stage of astronomical objects and to shed light on interstellar chemistry. These observations require a very accurate knowledge of the rotational spectrum, which is so far limited for mono-deuterated hydrogen sulfide, HDS. To fill this gap, high-level quantum-chemical calculations and sub-Doppler measurements have been combined for the investigation of the hyperfine structure of the rotational spectrum in the millimeter- and submillimeter-wave region. In addition to the determination of accurate hyperfine parameters, these new measurements together with the available literature data allowed us to extend the centrifugal analysis using a Watson-type Hamiltonian and a Hamiltonian-independent approach based on the Measured Active Ro-Vibrational Energy Levels (MARVEL) procedure. The present study thus permits to model the rotational spectrum of HDS from the microwave to far-infrared region with great accuracy, thereby accounting for the effect of the electric and magnetic interactions due to the deuterium and hydrogen nuclei.

Toward Accurate yet Effective Computations of Rotational Spectroscopy Parameters for Biomolecule Building Blocks.

The interplay of high-resolution rotational spectroscopy and quantum-chemical computations plays an invaluable role in the investigation of biomolecule building blocks in the gas phase. However, quantum-chemical methods suffer from unfavorable scaling with the dimension of the system under consideration. While a complete characterization of flexible systems requires an elaborate multi-step strategy, in this work, we demonstrate that the accuracy obtained by quantum-chemical composite approaches in the prediction of rotational spectroscopy parameters can be approached by a model based on density functional theory. Glycine and serine are employed to demonstrate that, despite its limited cost, such a model is able to predict rotational constants with an accuracy of 0.3% or better, thus paving the way toward the accurate characterization of larger flexible building blocks of biomolecules.

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

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.

Parity-pair-mixing effects in nonlinear spectroscopy of HDO.

A non-linear spectroscopic study of the HDO molecule is performed in the wavelength range of 1.36-1.42 µm using noise-immune cavity-enhanced optical-heterodyne molecular spectroscopy (NICE-OHMS). More than 100 rovibrational Lamb dips are recorded, with an experimental precision of 2-20 kHz, related to the first overtone of the O-H stretch fundamental of HD16O and HD18O. Significant perturbations, including distortions, shifts, and splittings, have been observed for a number of Lamb dips. These spectral perturbations are traced back to an AC-Stark effect, arising due to the strong laser field applied in all saturation-spectroscopy experiments. The AC-Stark effect mixes parity pairs, that is pairs of rovibrational states whose assignment differs solely in the Kc quantum number, where Kc is part of the standard J K a,K c asymmetric-top rotational label. Parity-pair mixing seems to be especially large for parity pairs with Ka ≥ 3, whereby their energy splittings become as small as a few MHz, resulting in multi-component asymmetric Lamb-dip profiles of gradually increasing complexity. These complex profiles often include crossover resonances. This effect is well known in saturation spectroscopy, but has not been reported in combination with parity-pair mixing. Parity-pair mixing is not seen in H2 16O and H2 18O, because their parity pairs correspond to ortho and para nuclear-spin isomers, whose interaction is prohibited. Despite the frequency shifts observed for HD16O and HD18O, the absolute accuracy of the detected transitions still exceeds that achievable by Doppler-limited techniques.

Exploiting the "Lego brick" approach to predict accurate molecular structures of PAHs and PANHs.

Polycyclic aromatic hydrocarbons (PAHs) and polycyclic aromatic nitrogen heterocycles (PANHs) are important and ubiquitous species in space. However, their accurate structural and spectroscopic characterization is often missing. To fill this gap, we exploit the so-called "Lego brick" approach [Melli et al., J. Phys. Chem. A, 2021, 125, 9904] to evaluate accurate rotational constants of some astrochemically relevant PAHs and PANHs. This model is based on the assumption that a molecular system can be seen as formed by smaller fragments for which a very accurate equilibrium structure is available. Within this model, the "template molecule" (TM) approach is employed to account for the modifications occurring when going from the isolated fragment to the molecular system under investigation, with the "linear regression" model being exploited to correct the linkage between different fragments. In the present work, semi-experimental equilibrium structures are used within the TM model. The performance of the "Lego brick" approach has been first tested for a set of small PA(N)Hs for which experimental data are available, thus leading to the conclusion that it is able to provide rotational constants with a relative accuracy well within 0.1%. Subsequently, it has been extended to the accurate prediction of the rotational constants for systems lacking any spectroscopic characterization.

Dipolar spin-spin coupling as an auxiliary tool for the structure determination of small isolated molecules.

The "gold standard" for obtaining accurate equilibrium structures is the so-called semi-experimental (SE) approach, which exploits the structural information contained in rotational constants. Within the SE approach, ground-state rotational constants-accurately obtained from high-resolution spectroscopic studies-are computationally corrected in order to remove vibrational effects. The resulting SE equilibrium rotational constants for a significant set of isotopic species allow for retrieving a unique set of equilibrium bond lengths and angles for the molecule under consideration. However, in some cases, the lack of isotopic substitution hampers or even prevents a rigorous and complete structure determination. In this perspective, we introduce the use of dipolar spin-spin coupling constants as an additional source of structural information in support of the standard SE approach. As a proof-of-concept, we tested this new strategy on some prototypical species, such as water, ammonia, phosphine, and their fluorinated counterparts. Our results indicate that-even when the molecular structure can be obtained from a large set of isotopic rotational constants-the use of dipolar spin-spin coupling constants guarantees a better accuracy and reduces the correlations among the geometrical parameters. Moreover, we point out that our approach offers the possibility to fully derive the molecular structure of PF3, a species for which any isotopic substitution is not possible.

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

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.

Rotational Spectra of Unsaturated Carbon Chains Produced by Pyrolysis: The Case of Propadienone, Cyanovinylacetylene, and Allenylacetylene.

Several interstellar molecules are highly reactive unsaturated carbon chains, which are unstable under terrestrial conditions. Laboratory studies in support of their detection in space thus face the issue of how to produce these species and how to correctly model their rotational energy levels. In this work, we introduce a general approach for producing and investigating unsaturated carbon chains by means of selected test cases. We report a comprehensive theoretical/experimental spectroscopic characterization of three species, namely, propadienone, cyanovinylacetylene, and allenylacetylene, all of them being produced by means of flash vacuum pyrolysis of a suitable precursor. For each species, quantum-chemical calculations have been carried out with the aim of obtaining accurate predictions of the missing spectroscopic information required to guide spectral analysis and assignment. Rotational spectra of the title molecules have been investigated up to 400 GHz by using a frequency-modulation millimeter-/submillimeter-wave spectrometer, thus significantly extending spectral predictions over a wide range of frequency and quantum numbers. A comparison between our results and those available in the literature points out the clear need of the reported laboratory measurements at higher frequencies for setting up accurate line catalogs for astronomical searches.