Carbamic acid and its dimer: A computational study.

A recent work by Marks et al. on the formation of carbamic acid in NH 3 $$ {}_3 $$ -CO 2 $$ {}_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.

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.

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%.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Intersystem crossing in the entrance channel of the reaction of O(3P) with pyridine.

Two quantum effects can enable reactions to take place at energies below the barrier separating reactants from products: tunnelling and intersystem crossing between coupled potential energy surfaces. Here we show that intersystem crossing in the region between the pre-reactive complex and the reaction barrier can control the rate of bimolecular reactions for weakly coupled potential energy surfaces, even in the absence of heavy atoms. For O(3P) plus pyridine, a reaction relevant to combustion, astrochemistry and biochemistry, crossed-beam experiments indicate that the dominant products are pyrrole and CO, obtained through a spin-forbidden ring-contraction mechanism. The experimental findings are interpreted-by high-level quantum-chemical calculations and statistical non-adiabatic computations of branching fractions-in terms of an efficient intersystem crossing occurring before the high entrance barrier for O-atom addition to the N-atom lone pair. At low to moderate temperatures, the computed reaction rates prove to be dominated by intersystem crossing.