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.

A new chapter in the never ending story of cycloadditions: The puzzling case of SO2 and acetylene

A comprehensive study of the different classes of cycloaddition reactions ([3+2], [2+2], and [2+1]) of SO2 to acetylene and ethylene has been performed using density functional theory (DFT) and composite wavefunction methods. The [3+2] cycloaddition reaction, that was previously explored in the context of the cycloaddition of thioformaldehyde S-methylide (TSM) to ethylene and acetylene, proceeds in a concerted way to the formation of stable heterocycles. In this paper, we extend our study to the [2+2] and [2+1] cycloadditions of SO2 to acetylene, which would produce 1,1-oxathiete-2-oxide and thiirene-1,1-dioxide, respectively. One of the main conclusions is that cyclic 1,1-oxathiete-2-oxide can open through a relatively easy breaking of the SO single bond and rearrange toward sulfinyl acetaldehyde (SA). The SA molecule can easily undergo several internal rearrangements, which eventually lead to sulfenic acid and sulfoxide derivatives of ethenone, 1,2,3-dioxathiole, and CO plus sulfinylmethane. The most probable path, however, produces 2-thioxoacetic acid, whose derivatives (or those of the corresponding acetate) are usually obtained by Willgerodt-Kindler-type sulfuration of acetates. This product can in turn decompose, leading to the final products CO2 and H2CS. Comparison of this decomposition path with that of 2-amino-2-thioxoacetic acid shows that the process occurs through different H-transfer processes.

Revisiting Stereoselective Propene Polymerization Mechanisms: Insights through the Activation Strain Model

The stereoelectronic factors responsible for stereoselectivity in propene polymerization with several metallocene and post-metallocene transition metal catalysts have been revisited using a combined approach of DFT calculations, the Activation Strain Model, Natural Energy Decomposition Analysis and a molecular descriptor (%VBur). There are in most cases two different paths leading to the formation of stereoerrors (SE), and the classical model does not suffice to fully understand stereoregulation. Improving stereoselectivity requires raising the energies of both SE insertion transition states. Our analyses show that the degrees of deformation of the active site (catalyst+chain) and the prochiral monomer differ for these two paths, and between different catalyst classes. Based on such analyses we discuss: a) the subtle differences in SE formation between stereoselective catalysts with different ligand frameworks; b) the reason for exceptional stereoselectivity reported for a special ansa-metallocene catalyst; c) the (double) stereocontrol origin for isoselective catalysts; d) the electronic contribution for isoselective catalysts generating SE by a modification of the ligand wrapping mode during the polymerization. Although this study will not immediately suggest new catalyst structures, we believe that understanding stereoregulation in great detail will increase our chances of success.

Predicting the propene stereoselectivity on transition metal catalysts: A daunting task for density functional theory

Thanks to recent developments in hardware and software, quantum chemical methods are increasingly used for interpreting the complex mechanisms underlying polymerization reaction by homogeneous catalysis. Unfortunately, the dimensions of even the smallest realistic models are too large to permit the use of state-of-the-art composite wave function methods. Under these circumstances, density functional theory still offers the best compromise between cost and accuracy. However, comprehensive benchmarks of different functionals are not yet available for this important research field. The main aim of the present paper is to fill this gap by performing an unbiased comparison of several density functionals and continuum solvent models for the stereo-control in the propylene polymerization on prototypical catalysts inducing different reaction mechanisms. While it was not possible to define a unique computational protocol providing the best results in all the situations, the B3PW91 functional in conjunction with D3 empirical dispersions and the solvent model density solvent model performs remarkably well for three out of the four investigated catalysts. Under such circumstances, it is recommended to compare the results delivered by different models when approaching additional classes of catalysts.

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.

Hunting for Complex Organic Molecules in the Interstellar Medium: The Role of Accurate Low-Cost Theoretical Geometries and Rotational Constants

A new approach to computation at affordable cost of accurate geometrical structures and rotational constants for medium-sized molecules in the gas phase is further improved and applied to a large panel of interstellar complex organic molecules. The most distinctive feature of the new model is the effective inclusion of core-valence correlation and vibrational averaging effects in the framework of density functional theory (DFT). In particular, a double-hybrid functional in conjunction with a quadruple-ζ valence/triple-ζ polarization basis set is employed for geometry optimizations, whereas a cheaper hybrid functional in conjunction with a split-valence basis set is used for the evaluation of vibrational corrections. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme approaches the accuracy of state-of-the-art wave function methods with the computational cost of the standard methods (DFT or MP2) routinely employed in the interpretation of microwave spectra. Since the whole computational workflow involves the postprocessing of the output of standard electronic structure codes by a new freely available web utility, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.

LEADER local action groups and inner areas. An Italian case study

In this study we investigate marginality in European rural areas to understand how EU policies for local development can stimulate novel approaches to social innovation and economic growth. The method applied is a Social Network Analysis (SNA) on the call tenders of an Italian Local Action Group managing the LEADER programme in the Piacentino and Parmense Apennine (Emilia-Romagna). The implementation of the Inner Area Strategy in a part of the LEADER territory (11 out of the 59 municipalities of the LEADER area) allows to verify complementarities, interrelations and overlays between the two development policies.

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.

Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps

An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.

Accurate Thermochemical and Kinetic Parameters at Affordable Cost by Means of the Pisa Composite Scheme (PCS)

A new strategy for the computation at an affordable cost of geometrical structures, thermochemical parameters, and rate constants for medium-sized molecules in the gas phase is proposed. The most distinctive features of the new model are the systematic use of cc-pVnZ-F12 basis sets, the addition of MP2 core-valence correlation in geometry optimizations by a double-hybrid functional, the separate extrapolation of MP2 and post-MP2 contributions, and the inclusion of anharmonic contributions in zero-point energies and thermodynamic functions. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme outperforms the most well-known model chemistries without the need for any empirical parameter. Additional tests show that the computed zero-point energies and thermal contributions can be confidently used for obtaining accurate thermochemical and kinetic parameters. Since the whole computational workflow is translated in a black-box procedure, which can be followed with standard electronic structure codes, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.

Accurate thermochemistry at affordable cost by means of an improved version of the JunChS-F12 model chemistry

The junChS-F12 composite method has been improved by means of the latest implementation of the CCSD(F12*)(T+) ansatz and validated for the thermochemistry of molecules containing atoms of the first three rows of the periodic table. A thorough benchmark showed that this model, in conjunction with cost-effective revDSD-PBEP86-D3(BJ) reference geometries, offers an optimal compromise between accuracy and computational cost. If improved geometries are sought, the most effective option is to add MP2-F12 core-valence correlation corrections to CCSD(T)-F12b/jun-cc-pVTZ geometries without the need of performing any extrapolation to the complete basis set limit. In the same vein, CCSD(T)-F12b/jun-cc-pVTZ harmonic frequencies are remarkably accurate without any additional contribution. Pilot applications to noncovalent intermolecular interactions, conformational landscapes, and tautomeric equilibria confirm the effectiveness and reliability of the model.

DFT Meets Wave-Function Methods for Accurate Structures and Rotational Constants of Histidine, Tryptophan, and Proline

A new computational strategy has been applied to the conformational and spectroscopic properties in the gas phase of amino acids with very distinctive features, ranging from different tautomeric forms (histidine) to ring puckering (proline), and heteroaromatic structures with non-equivalent rings (tryptophan). The integration of modern double-hybrid functionals and wave-function composite methods has allowed us to obtain accurate results for a large panel of conformers with reasonable computer times. The remarkable agreement between computations and microwave experiments allows an unbiased interpretation of the latter in terms of stereoelectronic effects.

Accurate structures and spectroscopic parameters of α,α-dialkylated α-amino acids in the gas-phase: a joint venture of DFT and wave-function composite methods

Accurate computations of structural, conformational and spectroscopic properties in the gas phase have been performed for two α,α-dialkylated α-amino acids, namely aminoisobutyric acid and cyclopropylglycine. Thanks to the integration of modern double hybrid functionals and wave-function methods, several low-energy structures of the title molecules could be analyzed employing standard computer resources. The computed features of all the most stable conformers of the target amino acids closely match the corresponding spectroscopic parameters issued from microwave spectroscopic studies in the gas-phase. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the way for accurate investigations of other flexible bricks of life.

New prebiotic molecules in the interstellar medium from the reaction between vinyl alcohol and CN radicals: unsupervised reaction mechanism discovery, accurate electronic structure calculations and kinetic simulations

Vinyl alcohol (VyA) and cyanide (CN) radicals are relatively abundant in the interstellar medium (ISM). VyA is the enolic tautomer of acetaldehyde and has two low-lying conformers, characterized by the syn or anti placement of hydroxyl hydrogen with respect to the double bond. In this paper, we present a gas-phase model of the barrierless reactions of both VyA's conformers with CN employing accurate quantum chemical computations in the framework of a master equation approach based on the transition state theory. Our results indicate that both VyA conformers feature a similar reactivity with CN, starting with a barrierless addition to the double bond and followed by different isomerization, dissociation, and/or hydrogen elimination steps. The rate constants computed for temperatures up to 600 K show that several reaction channels are open even under the harsh conditions of the ISM, with the favoured one providing the first feasible formation route of a prebiotic molecule not yet detected in the ISM, namely cyanoacetaldehyde. This finding suggests looking for cyanoacetaldehyde in regions where both VyA and CN have already been detected, like, e.g., Sagittarius B2N or G+0.693-0.027.

PCS/Bonds and PCS0: Pick your molecule and get its accurate structure and ground state rotational constants at DFT cost

An unsupervised computational protocol is proposed with the aim of obtaining accurate structures of large molecules in the gas phase at the cost of standard density functional theory (DFT) computations. The whole workflow is fully automated and provides optimized equilibrium geometries and ground state rotational constants to be directly compared with experiments. The results for a panel of molecules of biological or medicinal interest show that the accuracy of the results delivered by the new tool at the cost of a single DFT geometry optimization is close to that delivered by state-of-the-art composite wavefunction methods for small semi-rigid molecules.

Accurate Structures and Spectroscopic Parameters of Guanine Tautomers in the Gas Phase by the Pisa Conventional and Explicitly Correlated Composite Schemes (PCS and PCS-F12)

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase is proposed and validated for tautomeric equilibria. The main features of the new model are the inclusion of core-valence correlation in geometry optimizations by a double hybrid functional and the systematic use of wave-function composite methods in conjunction with cc-pVnZ-F12 basis sets with separate extrapolation of MP2 and post-MP2 contributions. The resulting Pisa composite scheme employing conventional (PCS) or explicitly correlated (PCS-F12) approaches is applied to the challenging problem of guanine tautomers in the gas phase. The results are in remarkable agreement with the experimental structures, relative stabilities, and spectroscopic signatures of different tautomers. The accuracy of the results obtained at reasonable cost by means of black-box parameter-free approaches paves the way toward systematic investigations of other molecular bricks of life also by non-specialists.

DFT Meets Wave-Function Composite Methods for Characterizing Cytosine Tautomers in the Gas Phase

A general strategy for the accurate computation of structural and spectroscopic properties of biomolecule building blocks in the gas phase has been further improved and validated with a special reference to tautomeric equilibria. The main improvements concern the use of the cc-pVTZ-F12 basis set in both DFT and CCSD(T)-F12 computations, the inclusion of core-valence correlation in geometry optimizations by double hybrid functionals, and the use of the cc-pVQZ-F12 basis set for complete basis set extrapolation at the MP2-F12 level. The resulting model chemistry is applied to the challenging problem of cytosine tautomers in the gas phase. The results are in remarkable agreement with experiment concerning both rotational and vibrational spectroscopic parameters and permit their unbiased interpretation in terms of structural and thermochemical features. Together with the intrinsic interest of the studied molecule, the accuracy of the results obtained at reasonable cost without any empirical parameter suggests that the proposed composite method can be profitably employed for accurate investigations of other molecular bricks of life.

Accuracy Meets Feasibility for the Structures and Rotational Constants of the Molecular Bricks of Life: A Joint Venture of DFT and Wave-Function Methods

A fully unsupervised computational protocol is proposed with the aim of obtaining reliable structural properties for molecular bricks of life in the gas phase. The results of the new composite scheme approach spectroscopic accuracy at a moderate cost without any empirical parameter in addition to those of the underlying electronic structure method. The whole workflow is fully automated and provides optimized geometries and equilibrium rotational constants. Direct comparison with experimental ground state rotational constants can be performed thanks to the effective computation of vibrational corrections in the framework of second-order vibrational perturbation theory. The results for all the nucleic acid bases and several flexible molecules of biological or medicinal interest show that the accuracy of the new tool is close to that delivered by state-of-the-art composite wave function methods for small semirigid molecules.

Toward Spectroscopic Accuracy for the Structures of Large Molecules at DFT Cost: Refinement and Extension of the Nano-LEGO Approach

The SE100 database collecting accurate equilibrium geometries of medium size molecules obtained by the semiexperimental (SE) approach has been extended to species containing Br and I atoms. This has allowed the determination of accurate linear regressions between DFT and SE values for all the main bonds and angles involving H, B, C, N, O, F, P, S, Cl, Br, and I atoms. An improved Nano-LEGO tool has been developed, which is based on suitable hybrid and double hybrid functionals and combines in a fully coherent way the templating molecule and linear regression approaches. A number of case studies show that the new Nano LEGO tool provides geometrical parameters on par with state-of-the-art composite wave function methods, but can be routinely applied to medium- to large-size molecules. The accuracy reached for structural parameters is mirrored on rotational constants that can be predicted with an average error within 0.2%.

Reliable Gas Phase Reaction Rates at Affordable Cost by Means of the Parameter-Free JunChS-F12 Model Chemistry

A recently developed strategy for the computation at affordable cost of reliable barrier heights ruling reactions in the gas phase (junChS, [Barone, V.; J. Chem. Theory Comput. 2021, 17, 4913-4928]) has been extended to the employment of explicitly correlated (F12) methods. A thorough benchmark based on a wide range of prototypical reactions shows that the new model (referred to as junChS-F12), which employs cost-effective revDSD-PBEP86-D3(BJ) reference geometries, has an improved performance with respect to its conventional counterpart and outperforms the most well-known model chemistries without the need of any empirical parameter and at an affordable computational cost. Several benchmarks show that revDSD-PBEP86-D3(BJ) structures and force fields provide zero point energies and thermal contributions, which can be confidently used, together with junChS-F12 electronic energies, for obtaining accurate reaction rates in the framework of the master equation approach based on the ab initio transition-state theory.

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 spectroscopy of some representative molecular bricks of life.

Expected final online publication date for the Annual Review of Physical Chemistry, Volume 74 is April 2023. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Accurate Structures and Spectroscopic Parameters of Phenylalanine and Tyrosine in the Gas Phase: A Joint Venture of DFT and Composite Wave-Function Methods

A general strategy for the accurate computation of conformational and spectroscopic properties of flexible molecules in the gas phase is applied to two representative proteinogenic amino acids with aromatic side chains, namely, phenylalanine and tyrosine. The main features of all the most stable conformers predicted by this computational strategy closely match those of the species detected in microwave and infrared experiments. Together with their intrinsic interest, the accuracy of the results obtained with reasonable computer times paves the route for accurate investigations of other flexible bricks of life.

Bringing Machine-Learning Enhanced Quantum Chemistry and Microwave Spectroscopy to Conformational Landscape Exploration: the Paradigmatic Case of 4-Fluoro-Threonine

A combined experimental and theoretical study has been carried out on 4-fluoro-threonine, the only naturally occurring fluorinated amino acid. Fluorination of the methyl group significantly increases the conformational complexity with respect to the parent amino acid threonine. The conformational landscape has been characterized in great detail, with special attention given to the inter-conversion pathways between different conformers. This led to the identification of 13 stable low-energy minima. The equilibrium population of so many conformers produces a very complicated and congested rotational spectrum that could be assigned through a strategy that combines several levels of quantum chemical calculations with the principles of machine learning. Twelve conformers out of 13 could be experimentally characterized. The results obtained from the analysis of the intra-molecular interactions can be exploited to accurately model fluorine-substitution effects in biomolecules.

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.

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.

Perturb-Then-Diagonalize Vibrational Engine Exploiting Curvilinear Internal Coordinates

The present paper is devoted to the implementation and validation of a second-order perturbative approach to anharmonic vibrations, followed by variational treatment of strong couplings (GVPT2) based on curvilinear internal coordinates. The main difference with respect to the customary Cartesian-based formulation is that the kinetic energy operator is no longer diagonal, and has to be expanded as well, leading to additional terms which have to be taken into proper account. It is, however, possible to recast all the equations as well-defined generalizations of the corresponding Cartesian-based counterparts, thus achieving a remarkable simplification of the new implementation. Particular attention is paid to the treatment of Fermi resonances with significant number of test cases analyzed fully, validating the new implementation. The results obtained in this work confirm that curvilinear coordinates strongly reduce the strength of inter-mode couplings compared to their Cartesian counterparts. This increases the reliability of low-order perturbative treatments for semi-rigid molecules and paves the way toward the reliable representation of more flexible molecules where small- and large-amplitude motions can be safely decoupled and treated at different levels of theory.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

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.

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.

Unbiased disentanglement of conformational baths with the help of microwave spectroscopy, quantum chemistry and artificial intelligence: the puzzling case of homocysteine

An integrated experimental-computational strategy for the accurate characterization of the conformational landscape of flexible biomolecule building blocks is proposed. This is based on the combination of rotational spectroscopy with quantum-chemical computations guided by artificial intelligence tools. The first step of the strategy is the conformer search and relative stability evaluation performed by means of an evolutionary algorithm. In this step, last generation semiempirical methods are exploited together with hybrid and double-hybrid density functionals. Next, the barriers ruling the interconversion between the low-lying conformers are evaluated in order to unravel possible fast relaxation paths. The relative stabilities and spectroscopic parameters of the ``surviving' conformers are then refined using state-of-the-art composite schemes. The reliability of the computational procedure is further improved by the inclusion of vibrational and thermal effects. The final step of the strategy is the comparison between experiment and theory without any ad hoc adjustment, which allows an unbiased assignment of the spectroscopic features in terms of different conformers and their spectroscopic parameters. The proposed approach has been tested and validated for homocysteine, a highly flexible non-proteinogenic alpha-amino acid. The synergism of the integrated strategy allowed the characterization of five conformers stabilized by bifurcated N-H-O=C hydrogen bonds, together with an additional conformer involving a more conventional HNH-O hydrogen bond. The stability order estimated from the experimental intensities as well as the number and type of conformers observed in the gas phase are in full agreement with the theoretical predictions. Analogously, a good match has been found for the spectroscopic parameters.

Dipolar 1,3-cycloaddition of thioformaldehyde S-methylide (CH2 SCH2 ) to ethylene and acetylene. A comparison with (valence) isoelectronic O3 , SO2 , CH2 OO and CH2 SO

Methods rooted in the density functional theory and in the coupled cluster ansatz were employed to investigate the cycloaddition reactions to ethylene and acetylene of 1,3-dipolar species including ozone and the derivatives issued from replacement of the central oxygen atom by the valence-isoelectronic sulfur atom, and/or of one or both terminal oxygen atoms by the isoelectronic CH₂ group. This gives rise to five different 1,3-dipolar compounds, namely ozone itself (O₃ ), sulfur dioxide (SO₂ ), the simplest Criegee intermediate (CH₂ OO), sulfine (CH₂ SO), and thioformaldehyde S-methylide (CH₂ SCH₂ , TSM). The experimental and accurate theoretical data available for some of those molecules were employed to assess the accuracy of two last-generation composite methods employing conventional or explicitly correlated post-Hartree-Fock contributions (jun-Cheap and SVECV-f12, respectively), which were then applied to investigate the reactivity of TSM. The energy barriers provided by both composite methods are very close (the average values for the two composite methods are 7.1 and 8.3 kcal mol^-1 for the addition to ethylene and acetylene, respectively) and comparable to those ruling the corresponding additions of ozone (4.0 and 7.7 kcal mol^-1 , respectively). These and other evidences strongly suggest that, at least in the case of cycloadditions, the reactivity of TSM is similar to that of O₃ and very different from that of SO₂ .

Spectroscopic Characterization of 3-Aminoisoxazole, a Prebiotic Precursor of Ribonucleotides

The processes and reactions that led to the formation of the first biomolecules on Earth play a key role in the highly debated theme of the origin of life. Whether the first chemical building blocks were generated on Earth (endogenous synthesis) or brought from space (exogenous delivery) is still unanswered. The detection of complex organic molecules in the interstellar medium provides valuable support to the latter hypothesis. To gather more insight, here we provide the astronomers with accurate rotational frequencies to guide the interstellar search of 3-aminoisoxazole, which has been recently envisaged as a key reactive species in the scenario of the so-called RNA-world hypothesis. Relying on an accurate computational characterization, we were able to register and analyze the rotational spectrum of 3-aminoisoxazole in the 6–24 GHz and 80–320 GHz frequency ranges for the first time, exploiting a Fourier-transform microwave spectrometer and a frequency-modulated millimeter/sub-millimeter spectrometer, respectively. Due to the inversion motion of the −NH₂ group, two states arise, and both of them were characterized, with more than 1300 lines being assigned. Although the fit statistics were affected by an evident Coriolis interaction, we were able to produce accurate line catalogs for astronomical observations of 3-aminoisoxazole.

Gliding on Ice in Search of Accurate and Cost-Effective Computational Methods for Astrochemistry on Grains: The Puzzling Case of the HCN Isomerization

The isomerization of hydrogen cyanide to hydrogen isocyanide on icy grain surfaces is investigated by an accurate composite method (jun-Cheap) rooted in the coupled cluster ansatz and by density functional approaches. After benchmarking density functional predictions of both geometries and reaction energies against jun-Cheap results for the relatively small model system HCN···(H2O)2, the best performing DFT methods are selected. A large cluster containing 20 water molecules is then employed within a QM/QM' approach to include a realistic environment mimicking the surface of icy grains. Our results indicate that four water molecules are directly involved in a proton relay mechanism, which strongly reduces the activation energy with respect to the direct hydrogen transfer occurring in the isolated molecule. Further extension of the size of the cluster up to 192 water molecules in the framework of a three-layer QM/QM'/MM model has a negligible effect on the energy barrier ruling the isomerization. Computation of reaction rates by the transition state theory indicates that on icy surfaces, the isomerization of HNC to HCN could occur quite easily even at low temperatures thanks to the reduced activation energy that can be effectively overcome by tunneling.

FC073: An Epigenetically Driven Mechanism Triggered by Viral and Bacterial RNA Regulates the IL-6 Levels in IGA Nephropathy

BACKGROUND AND AIMS

Several models have been proposed to describe the pathogenesis of Immunoglobulin A nephropathy (IgAN) and, among them, the multihit model where the gut-microbiota may play an important role. These models explain the pathogenesis of IgAN caused by the production of aberrant IgA, but it is believed that further predisposing factors are present, including immunological, genetic, environmental or nutritional factors.

Recently, the role of IL-6 in IgAN pathogenesis is becoming increasingly important. It is essential for the deposition of glomerular immunoglobulin A and the development of renal disease in Cd37-deficient mice, although the pathogenetic mechanisms that determine its increase are not well known.

A possible hypothesis emerges from our recent work on genome-wide DNA methylation screening in patients with IgAN, which identified, among other findings, a hypermethylated region comprising Vault 2–1 RNA (VTRNA2-1), a non-RNA coding also known as a precursor of miR-886 (pre-mi-RNA). Consistently, VTRNA2-1 expression was found downregulated in IgAN patients.

Here we studied the involvement of the VTRNA2-1/PKR/CREB/IL-6 pathway in IgAN.

METHOD

Total RNA were isolated from PBMCs of IgAN patients, transplanted IgAN patients (TP-IgAN), non-IgAN transplanted patients (TP) and healthy subjects (HS). VTRNA2-1, CREB and PKR transcripts were evaluated by RT-PCR. Total and phosphorylated PKR, CREB and Il-6 proteins were evaluated by ELISA. Poly (I: C), a synthetic analogue of dsRNA and Pfizer-BioNTech COVID-19 COMIRNATY vaccine were used to transfect patient PBMCs. PKR inhibitor imoxin (C16) 1 µM was used to stimulate patient PBMCs.

RESULTS

Here we confirm that VTRNA2-1 transcript was down-regulated in native and transplanted IgAN subjects compared to HS and non IgAN transplanted patients, with a decrease of 30- and 100-folds, respectively (P < 0.05, and P < 0.0001). IgAN patients with downregulated VTRNA2-1 showed a PKR overactivation (fold increase of phosphorilation of 2.6- in IgAN and 2-folds in TP-IgAN patients; P < 0.05), coherently with the role played by VTRNA2-1 that binds to PKR and inhibits its phosphorylation. Then, we found that PKR causes the activation of CREB, a classical cAMP-inducible CRE-binding factor (fold increase of phosphorilation of 3- in IgAN and 2.67-folds in TP-IgAN patients; P < 0.01). CREB, interacting with a region of the IL-6 promoter, led to IL-6 production. Indeed, in IgAN patients we showed a IL-6 mean increase to 120 pg/mL compared to the respective controls (P < 0.05). Moreover, the IL-6 levels correlated with CREB and PKR phosphorylation (r = 0.97; P = 0.0006 and r = 0.89; P = 0.0064, respectively, for IgAN and TP-IgAN patients).

Since PKR is normally activated by bacterial and viral RNA, we hypothesized that these microorganisms can further activate the PKR/CREB/IL-6 pathway leading to an excess of IL-6 production. This may explain both the high levels of IL-6, and infection involvement in the disease, and cases of IgAN associated with COVID-19 infection or with COVID-19 RNA-vaccination, and recent data showing microbiota involvement in IgAN. Effectively, we found that IgAN PMBCs stimulated with RNA poly(I: C) or the COVID-19 RNA-vaccine showed a significant increase in IL-6 levels compared to not-stimulated PBMCs (P < 0.05), supporting the pathogentic role played by viral RNA in IgAN pathogenesis and explaining the cases of IgAN patients developing episodes of macrohematuria after a COVID-19 infection or vaccination.

Finally, we showed that the IL-6 secretion can be reduced by the PKR inhibitor imoxin (fold decrease of 5-folds in IgAN and TP-IgAN patients; P < 0.05).

CONCLUSION

In conclusion, the discovery of the upregulated VTRNA2-1/PKR/CREB/IL-6 pathway in IgAN patients may provide a new pathogenic mechanism in IgAN and may be useful for the development of novel therapeutic approaches, likely by modulating the VTRNA2-1 methylation level in IgAN patients.

Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection

The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH₂OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree–Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.

Development, Validation, and Pilot Application of a Generalized Fluctuating Charge Model for Computational Spectroscopy in Solution

A general approach enforcing nonperiodic boundary conditions for the computation of spectroscopic properties in solution has been improved including an effective description of charge-transfer contributions and coordination number adjustment for explicit solvent molecules. Both contributions are obtained from a continuous description of intermolecular hydrogen bonds, which has been employed also for an effective clustering of molecular dynamics trajectories. Fine tuning of the model has been performed for several water clusters, and then its efficiency and reliability have been demonstrated by computing the absorption spectra of different creatinine tautomers in aqueous solution.

Molecular Dynamics Simulations Enforcing Nonperiodic Boundary Conditions: New Developments and Application to the Solvent Shifts of Nitroxide Magnetic Parameters

Multiscale methods combining quantum mechanics and molecular mechanics (QM/MM) have become the most suitable and effective 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 a general protocol based on atomistic simulations, carried out under nonperiodic boundary conditions (NPBC). In particular, we add to our in house MD engine a new efficient treatment of mean field electrostatic contributions to energy and forces, together with the capability of performing the simulations either in the canonical (NVT) or in the isothermal–isobaric (NPT) ensemble. Next, we provide convincing evidence that the NBPC approach enhanced by specific tweaks for rigid body propagation, allows for the simulation of solute–solvent systems with a minimum number of degrees of freedom and large integration time step. After its validation, this new approach is applied to the challenging case of solvatochromic effects on the electron paramagnetic resonance (EPR) spectrum of a prototypical nitroxide radical. To this end, we propose and validate also an automated protocol to extract and weight simulation snapshots, making use of a continuous description of the strength of solute–solvent hydrogen bridges. While further developments are being worked on, the effectiveness of our approach, even in its present form, is demonstrated by the accuracy of the results obtained through an unsupervised approach characterized by a strongly reduced computational cost as compared to that of conventional QM/MM models, without any appreciable deterioration of accuracy.

Gas-phase identification of (Z)-1,2-ethenediol, a key prebiotic intermediate in the formose reaction

Prebiotic sugars are thought to be formed on primitive Earth by the formose reaction. However, their formation is not fully understood and it is plausible that key intermediates could have formed in extraterrestrial environments and subsequently delivered on early Earth by cometary bodies. 1,2-Ethenediol, the enol form of glycolaldehyde, represents a highly reactive intermediate of the formose reaction and is likely detectable in the interstellar medium. Here, we report the identification and first characterization of (Z)-1,2-ethenediol by means of rotational spectroscopy. The title compound has been produced in the gas phase by flash vacuum pyrolysis of bis-exo-5-norbornene-2,3-diol at 750 °C, through a retro-Diels-Alder reaction. The spectral analysis was guided by high-level quantum-chemical calculations, which predicted spectroscopic parameters in very good agreement with the experiment. Our study provides accurate spectral data to be used for searches of (Z)-1,2-ethenediol in the interstellar space.

Probing Liquid-Ordered and Disordered Phases in Lipid Model Membranes: A Combined Theoretical and Spectroscopic Study of a Fluorescent Molecular Rotor

An integrated theoretical/experimental strategy has been applied to the study of environmental effects on the spectroscopic parameters of 4-(diphenylamino)phtalonitrile (DPAP), a fluorescent molecular rotor. The computational part starts from the development of an effective force field for the first excited electronic state of DPAP and proceeds through molecular dynamics simulations in solvents of different polarities toward the evaluation of Stokes shifts by quantum mechanics/molecular mechanics (QM/MM) approaches. The trends of the computed results closely parallel the available experimental results thus giving confidence to the interpretation of new experimental studies of the photophysics of DPAP in lipid bilayers. In this context, results show unambiguously that both flexible dihedral angles and global rotations are significantly retarded in a cholesterol/DPPC lipid matrix with respect to the DOPC matrix, thus confirming the sensitivity of DPAP to probe different environments and, therefore, its applicability as a probe for detecting different structures and levels of plasma membrane organization.

An improved study of HCO+ and He system: Interaction potential, collisional relaxation, and pressure broadening

In light of its ubiquitous presence in the interstellar gas, the chemistry and reactivity of the HCO⁺ ion requires special attention. The availability of up-to-date collisional data between this ion and the most abundant perturbing species in the interstellar medium is a critical resource in order to derive reliable values of its molecular abundance from astronomical observations. This work intends to provide improved scattering parameters for the HCO⁺ and He collisional system. We have tested the accuracy of explicitly correlated coupled-cluster methods for mapping the short- and long-range multi-dimensional potential energy surface of atom-ion systems. A validation of the methodology employed for the calculation of the potential well has been obtained from the comparison with experimentally derived bound-state spectroscopic parameters. Finally, by solving the close-coupling scattering equations, we have derived the pressure broadening and shift coefficients for the first six rotational transitions of HCO⁺ as well as inelastic state-to-state transition rates up to j = 5 in the 5-100 K temperature interval.

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β2-adrenergic receptor (β2AR). The 3.2 Å x-ray crystal structure of human β2AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β2AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β2AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β2AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β2AR-ligand complexes performed here.

Enhancing the Accuracy of Ab Initio Molecular Dynamics by Fine Tuning of Effective Two-Body Interactions: Acetonitrile as a Test Case

Grimme's dispersion-corrected density functional theory (DFT-D) methods have emerged among the most practical approaches to perform accurate quantum mechanical calculations on molecular systems ranging from small clusters to microscopic and mesoscopic samples, i.e., including hundreds or thousands of molecules. Moreover, DFT-D functionals can be easily integrated into popular ab initio molecular dynamics (MD) software packages to carry out first-principles condensed-phase simulations at an affordable computational cost. Here, starting from the well-established D3 version of the dispersion-correction term, we present a simple protocol to improve the accurate description of the intermolecular interactions of molecular clusters of growing size, considering acetonitrile as a test case. Optimization of the interaction energy was performed with reference to diffusion quantum Monte Carlo calculations, successfully reaching the same inherent accuracy of the latter (statistical error of ∼0.1 kcal/mol per molecule). The refined DFT-D3 model was then used to perform ab initio MD simulations of liquid acetonitrile, again showing significant improvements toward available experimental data with respect to the default correction.

A Computational Journey across Nitroxide Radicals: From Structure to Spectroscopic Properties and Beyond

Nitroxide radicals are characterized by a long-lived open-shell electronic ground state and are strongly sensitive to the chemical environment, thus representing ideal spin probes and spin labels for paramagnetic biomolecules and materials. However, the interpretation of spectroscopic parameters in structural and dynamic terms requires the aid of accurate quantum chemical computations. In this paper we validate a computational model rooted into double-hybrid functionals and second order vibrational perturbation theory. Then, we provide reference quantum chemical results for the structures, vibrational frequencies and other spectroscopic features of a large panel of nitroxides of current biological and/or technological interest.

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.