VMS-ROT: A New Module of the Virtual Multifrequency Spectrometer for Simulation, Interpretation, and Fitting of Rotational Spectra

Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy

The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.

Quantifying Conformational Isomerism in Chain Molecules by Linear Raman Spectroscopy: The Case of Methyl Esters

The conformational preferences of the ester group have the potential to facilitate the large amplitude folding of long alkyl chains in the gas phase. They are monitored by Raman spectroscopy in supersonic jet expansions for the model system methyl butanoate, after establishing a quantitative relationship with quantum-chemical predictions for methyl methanoate. This requires a careful analysis of experimental details, and a simulation of the rovibrational contours for near-symmetric top molecules. The technique is shown to be complementary to microwave spectroscopy in quantifying coexisting conformations. It confirms that a C-O-C(=O)-C-C chain segment can be collapsed into a single all-trans conformation by collisional cooling, whereas alkyl chain isomerism beyond this five-membered chain largely survives the jet expansion. This sets the stage for the investigation of linear alkyl alkanoates in terms of dispersion-induced stretched-chain to hairpin transitions by Raman spectroscopy.

Looking for the Elusive Imine Tautomer of Creatinine: Different States of Aggregation Studied by Quantum Chemistry and Molecular Spectroscopy

New spectroscopic experiments and state-of-the-art quantum-chemical computations of creatinine in different aggregation states unequivocally unveiled a significant tuning of tautomeric equilibrium by the environment: from the exclusive presence of the amine tautomer in the solid state and aqueous solution to a mixture of amine and imine tautomers in the gas phase. Quantum-chemical calculations predict the amine species as the most stable tautomer by about 30 kJ mol^-1 in condensed phases. On the contrary, moving to the isolated forms, both Z and E imine isomers become more stable by about 7 kJ mol^-1 . Since the imine isomers and one amine tautomer are separated by significant energy barriers, all of them should be present in the gas phase. This prediction has indeed been confirmed by high-resolution rotational spectroscopy, which provides the first experimental characterization of the elusive imine tautomer. The interpretation of the complicated hyperfine structure of the rotational spectrum, originated by three ¹⁴ N nuclei, makes it possible to use the spectral signatures as a sort of fingerprint for each individual tautomer in the complex sample.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals

Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.

Structure, Stability, and Spectroscopic Properties of Small Acetonitrile Cation Clusters

Ionization and fragmentation pathways induced by ionizing agents are key to understanding the formation of complex molecules in astrophysical environments. Acetonitrile (CH₃CN), the simplest organic nitrile, is an important molecule present in the interstellar medium. In this work, DFT and MP2 calculations were performed in order to obtain the low energy structures of the most relevant cations formed from electron-stimulated ion desorption of CH₃CN ices. Selected reaction pathways and spectroscopic properties were also calculated. Our results indicate that the most stable acetonitrile cation structure is CH₂CNH⁺ and that hydrogenation can occur successively without isomerization steps until its complete saturation. Moreover, the stability of distinct cluster families formed from the interaction of acetonitrile with small fragments, such as CH_n⁺, C₂H_n⁺, and CH_nCNH⁺, is discussed in terms of their respective binding energies. Some of these molecular clusters are stabilized by hydrogen bonds, leading to species whose infrared features are characterized by a strong redshift of the N-H stretching mode. Finally, the rotational spectra of CH₃CN and protonated acetonitrile, CH₃CNH⁺, were simulated using distinct computational protocols based on DFT, MP2, and CCSD(T) considering centrifugal distortion, vibrational-rotational coupling, and vibrational anharmonicity corrections. By adopting an empirical scaling procedure for calculating spectroscopic parameters, we were able to estimate the rotational frequencies of CH₃CNH⁺ with an expected average error below 1 MHz for J values up to 10.

The Role of State-of-the-Art Quantum-Chemical Calculations in Astrochemistry: Formation Route and Spectroscopy of Ethanimine as a Paradigmatic Case

The gas-phase formation and spectroscopic characteristics of ethanimine have been re-investigated as a paradigmatic case illustrating the accuracy of state-of-the-art quantum-chemical (QC) methodologies in the field of astrochemistry. According to our computations, the reaction between the amidogen, NH, and ethyl, C2H5, radicals is very fast, close to the gas-kinetics limit. Although the main reaction channel under conditions typical of the interstellar medium leads to methanimine and the methyl radical, the predicted amount of the two E,Z stereoisomers of ethanimine is around 10%. State-of-the-art QC and kinetic models lead to a [E-CH3CHNH]/[Z-CH3CHNH] ratio of ca. 1.4, slightly higher than the previous computations, but still far from the value determined from astronomical observations (ca. 3). An accurate computational characterization of the molecular structure, energetics, and spectroscopic properties of the E and Z isomers of ethanimine combined with millimeter-wave measurements up to 300 GHz, allows for predicting the rotational spectrum of both isomers up to 500 GHz, thus opening the way toward new astronomical observations.

The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy

While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.

Spectroscopic evidence of a particular intermolecular interaction in iodomethane–ethanol mixtures: the cooperative effect of halogen bonding, hydrogen bonding, and the solvent effect

The particular intermolecular interaction of an iodomethane–ethanol mixture is revealed by NIR, Raman, DFT calculation, and 2D correlation analysis.

State-of-the-art computation of the rotational and IR spectra of the methyl-cyclopropyl cation: hints on its detection in space

Recent measurements by the Cassini Ion Neutral Mass Spectrometer demonstrated the presence of numerous carbocations in Titan's upper atmosphere. In [Ali et al., Planet. Space Sci., 2013, 87, 96], an analysis of these measurements revealed the formation of the three-membered cyclopropenyl cation and its methyl derivatives. As a starting point of a future coordinated effort of laboratory experiments, quantum-chemical calculations, and astronomical observations, in the present work the molecular structure and spectroscopic properties of the methyl-cyclopropenyl cation have been investigated by means of state-of-the-art computational approaches in order to simulate its rotational and infrared spectra. Rotational parameters have been predicted with an expected accuracy better than 0.1% for rotational constants and on the order of 1-2% for centrifugal-distortion terms. As for the infrared spectrum, despite the challenge of a large amplitude motion, fundamental transitions have been computed to a good accuracy, i.e., the uncertainties are expected to be smaller than 5-10 wavenumbers.

Accuracy of Rotational Parameters Predicted by High-Level Quantum-Chemical Calculations: Case Study of Sulfur-Containing Molecules of Astrochemical Interest

The accuracy of rotational parameters obtained from high-level quantum-chemical calculations is discussed for molecules containing second-row atoms. The main focus is on computed rotational constants for which two statistical analyses have been carried out. A first benchmark study concerns sulfur-bearing species and involves 15 molecules (for a total of 74 isotopologues). By comparing 15 different computational approaches, all of them based on the coupled-cluster singles and doubles approach (CCSD) augmented by a perturbative treatment of triple excitations, CCSD(T), we have analyzed the effects on computed rotational constants due to ( i) extrapolation to the complete basis-set limit, ( ii) correlation of core electrons, and ( iii) vibrational corrections to rotational constants. To extend the analysis to other molecules containing second-row elements, as well as to understand the effect of higher excitations, a second benchmark study involving 11 molecules (for a total of 54 isotopologues) has been performed. Finally, the rotational parameters of seven sulfur-containing molecules of astrochemical interest (CCS, C₃S, C₄S, C₅S, HCCS⁺, HC₄S⁺, and HOCS⁺/HSCO⁺) have been computed and compared to experimental data, when available, also addressing the direct comparison of simulated and experimental rotational spectra.

Towards the SMART workflow system for computational spectroscopy

Is it possible to convert highly specialized research in the field of computational spectroscopy into robust and user-friendly aids to experiments and industrial applications?