Diving for Accurate Structures in the Ocean of Molecular Systems with the Help of Spectroscopy and Quantum Chemistry

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

Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions

The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.

Accurate Structures and Rotational Constants of Steroid Hormones at DFT Cost: Androsterone, Testosterone, Estrone, β-Estradiol, and Estriol

A comprehensive analysis of the structural, conformational, and spectroscopic properties in the gas phase has been performed for five prototypical steroid hormones, namely, androsterone, testosterone, estrone, β-estradiol, and estriol. The revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion and a suitable triple-ζ basis set provides accurate conformational energies and equilibrium molecular structures, with the latter being further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental ground state rotational constants are reached when adding to those equilibrium values vibrational corrections obtained at the cost of standard harmonic frequencies thanks to the use of a new computational tool. Together with the intrinsic interest of the studied hormones, the accuracy of the results obtained at DFT cost for molecules containing about 50 atoms paves the way toward the accurate investigations of other flexible bricks of life.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

Accurate Geometries of Large Molecules by Integration of the Pisa Composite Scheme and the Templating Synthon Approach

An effective yet reliable computational workflow is proposed, which permits the computation of accurate geometrical structures for large flexible molecules at an affordable cost thanks to the integration of machine learning tools and DFT models together with reduced scaling computations of vibrational averaging effects. After validation of the different components of the overall strategy, a panel of molecules of biological interest have been analyzed. The results confirm that very accurate geometrical parameters can be obtained at reasonable cost for molecules including up to about 50 atoms, which are the largest ones for which comparison with high-resolution rotational spectra is possible. Since the whole computational workflow can be followed employing standard electronic structure codes, accurate results for large-sized molecules can be obtained at DFT cost also by nonspecialists.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

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.

Precise Equilibrium Structure of Benzene

Recent advances in gas-phase structure determination afford outstanding agreement between the CCSD(T)/cc-pCVTZ-corrected semi-experimental (reSE) equilibrium structures and their corresponding best theoretical estimates (BTEs) of the equilibrium structures (re) based upon corrections to the CCSD(T)/cc-pCV5Z geometries for the aromatic heterocycles pyrimidine and pyridazine. Herein, that same analysis is extended to the fundamental aromatic molecule benzene, using published experimental spectroscopic data for a total of 11 available isotopologues. The incorporation of rotational constants from all of these isotopologues and CCSD(T) corrections to address the impacts of both the vibration-rotation interaction and electron-mass distribution results in a highly precise and accurate reSE structure. The CCSD(T)/cc-pCV5Z optimized geometry has been further corrected to address a finite basis set, untreated electron correlation, relativistic effects, and a breakdown of the Born-Oppenheimer approximation. This analysis achieves outstanding agreement between the re (BTE) and reSE structural parameters of benzene to a highly satisfying level (0.0001 Å), an agreement that surpasses our recently published structures of the aforementioned nitrogen-substituted benzene analogues. The D6h geometry of benzene is now known to an unprecedented precision: RC-C = 1.3913 (1) Å and RC-H = 1.0809 (1) Å. The mutual agreement between theory and experiment presented in this work validates both, substantially resolving all discrepancies between the reSE and theoretical re structures available in the literature.

Potential Energy Surfaces Sampled in Cremer-Pople Coordinates and Represented by Common Force Field Functionals for Small Cyclic Molecules

The complex conformations of the cyclic moieties impact the physical and chemical properties of molecules. In this work, we chose 22 molecules of four-, five-, and six-membered rings and performed a thorough conformational sampling using Cremer-Pople coordinates. With consideration of symmetries, we obtained a total of 1504 conformational structures for four-membered, 5576 for five-membered, and 13509 for six-membered rings. All well-known and many less well-known conformers for each molecule were identified. We represented the potential energy surfaces (PESs) by fitting the data to common analytical force field (FF) functional forms. We found that the general features of PESs can be described by the essential FF functional forms; however, the accuracy of representation can be improved remarkably by including the torsion-bond and torsion-angle coupling terms. The best fit yields R-squared (R²) values close to 1.0 and mean absolute errors in energy less than 0.3 kcal/mol.

Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy

The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.

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.

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.

Integration of Quantum Chemistry, Statistical Mechanics, and Artificial Intelligence for Computational Spectroscopy: The UV-Vis Spectrum of TEMPO Radical in Different Solvents

The ongoing integration of quantum chemistry, statistical mechanics, and artificial intelligence is paving the route toward more effective and accurate strategies for the investigation of the spectroscopic properties of medium-to-large size chromophores in condensed phases. In this context we are developing a novel workflow aimed at improving the generality, reliability, and ease of use of the available computational tools. In this paper we report our latest developments with specific reference to unsupervised atomistic simulations employing non periodic boundary conditions (NPBC) followed by clustering of the trajectories employing optimized feature spaces. Next accurate variational computations are performed for a representative point of each cluster, whereas intracluster fluctuations are taken into account by a cheap yet reliable perturbative approach. A number of methodological improvements have been introduced including, e.g., more realistic reaction field effects at the outer boundary of the simulation sphere, automatic definition of the feature space by continuous perception of solute–solvent interactions, full account of polarization and charge transfer in the first solvation shell, and inclusion of vibronic contributions. After its validation, this new approach has been applied to the challenging case of solvatochromic effects on the UV–vis spectra of a prototypical nitroxide radical (TEMPO) in different solvents. The reliability, effectiveness, and robustness of the new platform is demonstrated by the remarkable agreement with experiment of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional quantum mechanics and molecular mechanics models without any accuracy reduction.

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 PF₃, a species for which any isotopic substitution is not possible.

Tackling an accurate description of molecular reactivity with double-hybrid density functionals

In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the “chemical accuracy.” In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties.

Quantification of Geometric Errors Made Simple: Application to Main-Group Molecular Structures

Nearly all electronic structure simulations begin with obtaining approximate geometries, making a systematic quantification of errors in approximate molecular structures of key importance. Recently, the geometric energy offset (GEO) framework based on a single and natural measure for quantifying and analyzing these errors has been proposed ( J. Phys. Chem. Lett. 2020, 11, 99579964). An accurate and far less costly approximation to GEO is utilized here to readily quantify errors in main-group structures and analyze them in a chemically intuitive way. The use of semiexperimental geometries as a reference further simplifies the analysis. The analysis reveals new insights into the geometric performance of methods, their rankings, as well as patterns across different classes of methods and basis sets that arise from the analysis.

Ab Initio Study of Fine and Hyperfine Interactions in Triplet POH

Phosphorous-containing molecules have a great relevance in prebiotic chemistry in view of the fact that phosphorous is a fundamental constituent of biomolecules, such as RNA, DNA, and ATP. Its biogenic importance has led astrochemists to investigate the possibility that P-bearing species could have formed in the interstellar medium (ISM) and subsequently been delivered to early Earth by rocky bodies. However, only two P-bearing molecules have been detected so far in the ISM, with the chemistry of interstellar phosphorous remaining poorly understood. Here, in order to shed further light on P-carriers in space, we report a theoretical spectroscopic characterisation of the rotational spectrum of POH in its 3A″ ground electronic state. State-of-the-art coupled-cluster schemes have been employed to derive rotational constants, centrifugal distortion terms, and most of the fine and hyperfine interaction parameters, while the electron spin-spin dipolar coupling has been investigated using the multi-configuration self-consistent-field method. The computed spectroscopic parameters have been used to simulate the appearance of triplet POH rotational and ro-vibrational spectra in different conditions, from cold to warm environments, either in gas-phase experiments or in molecular clouds. Finally, we point out that the predicted hyperfine structures represent a key pattern for the recognition of POH in laboratory and interstellar spectra.

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

Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations

We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.

F12-TZ-cCR: A Methodology for Faster and Still Highly Accurate Quartic Force Fields

The F12-TZ-cCR quartic force field (QFF) methodology, defined here as CCSD(T)-F12b/cc-pCVTZ-F12 with further corrections for relativity, is introduced as a cheaper and even more accurate alternative to more costly composite QFF methods like those containing complete basis set extrapolations within canonical coupled cluster theory. F12-TZ-cCR QFFs produce B₀ and C₀ vibrationally averaged principal rotational constants within 7.5 MHz of gas-phase experimental values for tetraatomic and larger molecules, offering higher accuracy in these constants than the previous composite methods. In addition, F12-TZ-cCR offers an order of magnitude decrease in the computational cost of highly accurate QFF methodologies accompanying this increase in accuracy. An additional order of magnitude in cost reduction is achieved in the F12-DZ-cCR method, while also matching the accuracy of the traditional composite method's B₀ and C₀ constants. Finally, F12-DZ and F12-TZ are benchmarked on the same test set, revealing that both methods can provide anharmonic vibrational frequencies that are comparable in accuracy to all three of the more expensive methodologies, although their rotational constants lag behind. Hence, the present work demonstrates that highly accurate theoretical rovibrational spectral data can be obtained for a fraction of the cost of conventional QFF methodologies, extending the applicability of QFFs to larger molecules.

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)^MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol^-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH₂, Ac-Ala-NH₂, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)^MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.

Extending the Applicability of the Semi-experimental Approach by Means of "Template Molecule" and "Linear Regression" Models on Top of DFT Computations

The accurate determination of equilibrium structures for isolated molecules plays a central role in the evaluation and interpretation of stereoelectronic, thermodynamic, and spectroscopic properties. For small semi-rigid systems, state-of-the-art quantum-chemical computations can rival the most sophisticated experimental results. For larger molecules, cheaper yet accurate approaches need to be defined. The double-hybrid rev-DSD-PBEP86 functional already delivers remarkable results that can be further improved by means of a "Lego brick" model. This is based on the idea that a molecular system can be seen as formed by different fragments (the "Lego bricks"), whose accurate semi-experimental (SE) equilibrium geometries are available. The template molecule (TM) approach can be used to account for the modifications occurring when going from the isolated fragment to the molecular system under investigation, with the linear regression (LR) model employed to correct the linkage between the different fragments. The resulting TM-SE_LR approach has been tested with respect to available SE equilibrium structures and rotational constants. Indeed, the latter parameters straightforwardly depend on the equilibrium geometry of the system under consideration. The main outcome of our study is the reliability, robustness, and accuracy of this novel approach. The molecular systems considered for benchmarking the TM-SE_LR scheme are those formally issued from addition/elimination reactions of nucleophilic unsaturated radicals (e.g., CN, C2H, and phenyl) to alkenes, imines, and aldehydes, whose rotational spectra have been investigated, but accurate structural determinations are not yet available.

junChS and junChS-F12 Models: Parameter-free Efficient yet Accurate Composite Schemes for Energies and Structures of Noncovalent Complexes

A recently developed model chemistry (denoted as junChS [Alessandrini, S.; et al. J. Chem. Theory Comput. 2020, 16, 988-1006]) has been extended to the employment of explicitly correlated (F12) methods. This led us to propose a family of effective, reliable, and parameter-free schemes for the computation of accurate interaction energies of molecular complexes ruled by noncovalent interactions. A thorough benchmark based on a wide range of interactions showed that the so-called junChS-F12 model, which employs cost-effective revDSD-PBEP86-D3(BJ) reference geometries, has an improved performance with respect to its conventional counterpart and outperforms well-known model chemistries. Without employing any empirical parameter and at an affordable computational cost, junChS-F12 reaches subchemical accuracy. Accurate characterizations of molecular complexes are usually limited to energetics. To take a step forward, the conventional and F12 composite schemes developed for interaction energies have been extended to structural determinations. A benchmark study demonstrated that the most effective option is to add MP2-F12 core-valence correlation corrections to fc-CCSD(T)-F12/jun-cc-pVTZ geometries without the need of recovering the basis set superposition error and the extrapolation to the complete basis set.

Accurate Biomolecular Structures by the Nano-LEGO Approach: Pick the Bricks and Build Your Geometry

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical–chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of the role played by dynamical and environmental effects in tuning their overall behavior. For small semi-rigid systems in the gas phase, state-of-the-art quantum chemical computations rival the most sophisticated experimental (from, for example, high-resolution spectroscopy) results. For larger molecules, more effective computational approaches must be devised. To this end, we have further enlarged the compilation of available semi-experimental (SE) equilibrium structures, now covering the most important fragments containing H, B, C, N, O, F, P, S, and Cl atoms collected in the new SE100 database. Next, comparison with geometries optimized by methods rooted in the density functional theory showed that the already remarkable results delivered by PW6B95 and, especially, rev-DSDPBEP86 functionals can be further improved by a linear regression (LR) approach. Use of template fragments (taken from the SE100 library) together with LR estimates for the missing interfragment parameters paves the route toward accurate structures of large molecules, as witnessed by the very small deviations between computed and experimental rotational constants. The whole approach has been implemented in a user-friendly tool, termed nano-LEGO, and applied to a number of demanding case studies.

Theoretical insights into the effect of halogenated substituent on the electronic structure and spectroscopic properties of the favipiravir tautomeric forms and its implications for the treatment of COVID-19

In this study, we systematically investigated the electronic structure, spectroscopic (nuclear magnetic resonance, infrared, Raman, electron ionization mass spectrometry, UV-Vis, circular dichroism, and emission) properties, and tautomerism of halogenated favipiravir compounds (fluorine, chlorine, and bromine) from a computational perspective. Additionally, the effects of hydration on the proton transfer mechanism of the tautomeric forms of the halogenated favipiravir compounds are discussed. Our results suggest that spectroscopic properties allow for the elucidation of such tautomeric forms. As is well-known, the favipiravir compound has excellent antiviral properties and hence was recently tested for the treatment of new coronavirus (SARS-CoV-2). Through in silico modeling, in the current study, we evaluate the role of such tautomeric forms in order to consider the effect of drug-metabolism in the inhibition process of the main protease (M^pro) and RNA-dependent RNA polymerase (RdRp) of SARS-CoV-2 virus. According to the molecular docking, all halogenated compounds presented a better interaction energy than the co-crystallized active ligand (−3.5 kcal mol⁻¹) in the viral RdRp, in both wild-type (−6.3 to −6.5 kcal mol⁻¹) and variant (−5.4 to −5.6 kcal mol⁻¹) models. The variant analyzed for RdRp (Y176C) decreases the affinity of the keto form of the compounds in the active site, and prevented the ligands from interacting with RNA. These findings clearly indicated that all these compounds are promising as drug candidates for this molecular target.

In this study, we systematically investigated the electronic structure, spectroscopic properties, and tautomerism of halogenated favipiravir compounds (fluorine, chlorine, and bromine) from a computational perspective.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

Semi-Experimental Equilibrium (reSE) and Theoretical Structures of Pyridazine (o-C4H4N2)

A semi-experimental equilibrium structure (r_e^SE) of pyridazine (o-C₄H₄N₂) has been determined using the rotational spectra of 18 isotopologues. Spectroscopic constants of four isotopologues are reported for the first time (measured from 235 to 360 GHz), while spectroscopic constants for previously reported isotopologues are improved by extending the frequency coverage (measured from 130 to 375 GHz). The experimental values of the ground-state rotational constants (A₀, B₀, and C₀) from each isotopologue were converted to determinable constants (A₀^″, B₀^″, and C₀^″), which were then corrected for the effects of vibration-rotation interactions and electron-mass distributions using CCSD(T)/cc-pCVTZ calculations. The resultant r_e^SE for pyridazine determines bond distances to within 0.001 Å and bond angles within 0.04°, a reduction in the statistical uncertainties by at least a factor of two relative to the previously reported r_e^SE. The improvement in precision appears to be largely due to the use of higher-level theoretical calculations of the vibration-rotation and electron-mass effects, though the incorporation of the newly measured isotopologues ([4-²H, 4-¹³C]-, [4-²H, 5-¹³C]-, [4-²H, 6-¹³C]-, and [4,5-²H, 4-¹³C]-pyridazine) is partially responsible for the improved determination of the hydrogen-containing bond angles. The computed equilibrium structure (r_e) (CCSD(T)/cc-pCV5Z) and a "best theoretical estimate" of the equilibrium structure (r_e) both agree with the updated r_e^SE structure within the statistical experimental uncertainty (2σ) of each structural parameter.

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.

Unveiling Bifunctional Hydrogen Bonding with the Help of Quantum Chemistry: The Imidazole-Water Adduct as Test Case

The ubiquitous role of water and its amphiprotic nature call for a deeper insight into the physical-chemical properties of hydrogen-bonded complexes formed with building blocks of biomolecules. In this work, the semiexperimental (SE) approach combined with the template model (TM) protocol allowed the accurate determination of the equilibrium structure of two isomeric forms of the imidazole-water complex. In this procedure, the integration of experiment (thanks to a recent rotational spectroscopy investigation) and theory is exploited, also providing the means of assessing the reliability and accuracy of different quantum-chemical approaches. Overall, this study demonstrated the robustness of the combined SE-TM approach, which can provide accurate results using affordable quantum-chemical methods. Finally, the structural and energetic characteristics of these complexes have been examined in detail and compared with those of analogous heterocycle-water adducts, also exploiting energy decomposition analyses.

Highly-accurate quartic force fields for the prediction of anharmonic rotational constants and fundamental vibrational frequencies

The CcCR quartic force field (QFF) methodology is capable of computing B₀ and C₀ rotational constants to within 35 MHz (0.14%) of experiment for triatomic and larger molecules with at least two heavy atoms. Additionally, the same constants for molecules with four or more atoms agree to within 20 MHz (0.12%) of experiment for the current test set. This work also supports previous claims that the same QFF methodology can produce fundamental vibrational frequencies with a deviation less than 5.7 cm^-1 from experiment. Consequently, this approach of augmenting complete basis set extrapolated energies with treatments of core electron correlation and scalar relativity produces some of the most accurate rovibrational spectroscopic data available.

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.

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.

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.

Surprisingly broad applicability of the cc-pVnZ-F12 basis set for ground and excited states

Excellent convergence properties for the (aug-)cc-pVnZ-F12 basis set family, purpose-made for explicitly correlated calculations, are demonstrated with conventional wave function methods and Kohn-Sham density functional theory for various ground and excited-state calculations. Among the ground-state properties studied are dipole moments, covalent bond lengths, and interaction and reaction energies. For excited states, we looked at vertical excitation energies, UV absorption, and excited-state absorption spectra. Convergence is compared against the basis sets cc-pVnZ, def2-nVD, aug-pcseg-n, and nZaPa-NR. It is established that the cc-pVnZ-F12 family consistently yields results of n + 1 quality and better. Especially, the cc-pVDZ-F12 basis set is found to be a basis set of good cost vs performance trade-off.

The challenge of non-covalent interactions: theory meets experiment for reconciling accuracy and interpretation

In the past decade, many gas-phase spectroscopic investigations have focused on the understanding of the nature of weak interactions in model systems. Despite the fact that non-covalent interactions play a key role in several biological and technological processes, their characterization and interpretation are still far from being satisfactory. In this connection, integrated experimental and computational investigations can play an invaluable role. Indeed, a number of different issues relevant to unraveling the properties of bulk or solvated systems can be addressed from experimental investigations on molecular complexes. Focusing on the interaction of biological model systems with solvent molecules (e.g., water), since the hydration of the biomolecules controls their structure and mechanism of action, the study of the molecular properties of hydrated systems containing a limited number of water molecules (microsolvation) is the basis for understanding the solvation process and how structure and reactivity vary from gas phase to solution. Although hydrogen bonding is probably the most widespread interaction in nature, other emerging classes, such as halogen, chalcogen and pnicogen interactions, have attracted much attention because of the role they play in different fields. Their understanding requires, first of all, the characterization of the directionality, strength, and nature of such interactions as well as a comprehensive analysis of their competition with other non-covalent bonds. In this review, it is shown how state-of-the-art quantum-chemical computations combined with rotational spectroscopy allow for fully characterizing intermolecular interactions taking place in molecular complexes from both structural and energetic points of view. The transition from bi-molecular complex to microsolvation and then to condensed phase is shortly addressed.

Theory meets experiment for elucidating the structure and stability of non-covalent complexes: water–amine interaction as a proof of concept

A joint experimental-theoretical spectroscopic investigation has focused on a better understanding of the nature of weak, non-covalent interactions in amine-water model systems.

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.

Conformational Equilibria and Molecular Structures of Model Sulfur-Sulfur Bridge Systems: Diisopropyl Disulfide

The conformations and molecular structures of diisopropyl disulfide have been studied by high-resolution microwave spectroscopy and quantum chemical calculations. Three conformers, G'GG', G'GT, and GGG', have been observed in the jet expansion. The global minimum, G'GG', adopts a configuration with the G' orientation of H-C-S-S and S-S-C-H and the G orientation of C-S-S-C showing the C₂ symmetry. The rotational spectra of monosubstituted ¹³C and ³⁴S isotopologues have also been recorded for G'GG', leading to an accurate structural determination of this conformer. Two additional ³⁴S isotopologues have also been measured for G'GT. The relative energies of three observed conformers calculated at the MP2/6-311++(d,p) level of theory are within 2 kJ mol^-1, while the relative intensity measurements suggested their population ratio to be N_G'GG'/N_G'GT/N_GGG' ≈ 5:3:2.

Astrochemistry and Astrobiology: Materials Sciencein Wonderland?

Astrochemistry and astrobiology, the fascinating disciplines that strive to unravel the origin of life, have opened unprecedented and unpredicted vistas into exotic compounds as well as extreme or complex reaction conditions of potential relevance for a broad variety of applications. Representative, and so far little explored sources of inspiration include complex organic systems, such as polycyclic aromatic hydrocarbons (PAHs) and their derivatives; hydrogen cyanide (HCN) and formamide (HCONH2) oligomers and polymers, like aminomalononitrile (AMN)-derived species; and exotic processes, such as solid-state photoreactions on mineral surfaces, phosphorylation by minerals, cold ice irradiation and proton bombardment, and thermal transformations in fumaroles. In addition, meteorites and minerals like forsterite, which dominate dust chemistry in the interstellar medium, may open new avenues for the discovery of innovative catalytic processes and unconventional methodologies. The aim of this review was to offer concise and inspiring, rather than comprehensive, examples of astrochemistry-related materials and systems that may be of relevance in areas such as surface functionalization, nanostructures, and hybrid material design, and for innovative technological solutions. The potential of computational methods to predict new properties from spectroscopic data and to assess plausible reaction pathways on both kinetic and thermodynamic grounds has also been highlighted.

Potential-Energy Surfaces for Ring-Puckering Motions of Flexible Cyclic Molecules through Cremer-Pople Coordinates: Computation, Analysis, and Fitting

Ring-puckering motion in 12 flexible cyclic molecules is investigated by calculation and analysis of two-dimensional potential-energy surfaces (PESs) using the so-called ring-puckering coordinates proposed by Cremer and Pople. The PESs are calculated by means of density-functional theory using a B2PLYP-D3BJ exchange-correlation functional with a maug-cc-pVTZ basis set, and results are compared to the available experimental and theoretical data. Special care is devoted to the aspect of symmetry in such two-dimensional PESs, which are here reported for the first time also for molecules whose planar form has symmetry lower than D_{5 h} or C_{2 v}. The issue of PES fitting and that of solving the nuclear dynamics using ring-puckering coordinates are also addressed. Analytical formulations of the computed PESs using suitable functional forms with a limited set of parameters are provided.