Minimally Empirical Double-Hybrid Functionals Trained against the GMTKN55 Database: revDSD-PBEP86-D4, revDOD-PBE-D4, and DOD-SCAN-D4

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.

Neighboring-Group Participation by C-2 Acyloxy Groups: Influence of the Nucleophile and Acyl Group on the Stereochemical Outcome of Acetal Substitution Reactions

A single acyloxy group at C-2 can control the outcome of nucleophilic substitution reactions of pyran-derived acetals, but the extent of the neighboring-group participation depends on a number of factors. We show here that neighboring-group participation does not necessarily control the stereochemical outcome of acetal substitution reactions with weak nucleophiles. The 1,2-trans selectivity increased with increasing reactivity of the incoming nucleophile. This trend suggests the intermediacy of both cis-fused dioxolenium ions and oxocarbenium ions in the stereochemistry-determining step. In addition, as the electron-donating ability of the neighboring group decreased, the preference for the 1,2-trans products increased. Computational studies show how the barriers for the ring-opening reaction on the dioxolenium ions and the transition states to provide the oxocarbenium ions change with the electron-donating capacity of the C-2-acyloxy group and the reactivity of the nucleophile.

High-Level Quantum Chemical Prediction of C-F Bond Dissociation Energies of Perfluoroalkyl Substances

In this study, 550 C-F bond dissociation energies (BDEs) of a variety of per- and polyfluoroalkyl substances (PFASs) obtained from high-level DLPNO-CCSD(T)/CBS calculations were used to assess the accuracy of contemporary density functional theory (DFT) and semiempirical methods. DLPNO-CCSD(T)/CBS gas phase C-F BDEs fall between 404.9-550.7 kJ mol-1 and M06-2X and ωB97M-V in conjunction with the aug-cc-pVTZ basis set predicted BDEs closest to the benchmark level with a mean absolute deviation (MAD) of 7.3 and 8.3 kJ mol-1, respectively. It was observed that DFT prediction errors increase with the degree of fluorination and system size. As such, previous model chemistry recommendations based on smaller nonfluorinated systems may not be carried over to modeling the energetics of PFASs and related systems.

Regularized and Opposite Spin-Scaled Functionals from Møller-Plesset Adiabatic Connection─Higher Accuracy at Lower Cost

Noncovalent interactions (NCIs) play a crucial role in biology, chemistry, material science, and everything in between. To improve pure quantum-chemical simulations of NCIs, we propose a methodology for constructing approximate correlation energies by combining an interpolation along the Møller-Plesset adiabatic connection (MP AC) with a regularization and spin-scaling strategy applied to MP2 correlation energies. This combination yields cosκos-SPL2, which exhibits superior accuracy for NCIs compared to any of the individual strategies. With the N4 formal scaling, cosκos-SPL2 is competitive or often outperforms more expensive dispersion-corrected double hybrids for NCIs. The accuracy of cosκos-SPL2 particularly shines for anionic halogen bonded complexes, where it surpasses standard dispersion-corrected DFT by a factor of 3 to 5.

Hunting for interstellar molecules: rotational spectra of reactive species

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

Stereodivergent Zweifel Olefination and its Mechanistic Dichotomy

Stereoselective Zweifel olefination using boronate complexes carrying two different reactive π-systems was achieved to synthesize vinyl heteroarenes and conjugated 1,3-dienes in good yield and up to 100 % stereoselectivity, which remains unexplored until now. Most importantly, we report the unprecedented formation of E vs. Z-vinyl heteroarenes for different heteroarenes under identical conditions. Density functional theory (DFT) investigations unveil the mechanistic dichotomy between olefin and heteroarene activation followed by 1,2-migration, leading to E or Z-vinyl heteroarenes respectively. We also report a previously unknown reversal of stereoselectivity by using 2,3-Dichloro-5,6-dicyano-1,4-benzoquinone (DDQ) as an electrophile. The Zweifel olefination using a boronate complex that carries two different olefins was previously unexplored due to significant challenges associated with the site-selective activation of olefins. We have solved this problem and reported the site-selective activation of olefins for the stereoselective synthesis of 1,3-dienes.

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.

Origins of Enhanced Enantioselectivity in the Pd-Catalyzed Decarboxylative Allylic Alkylation of N-Benzoyl Lactams

We explore the origins of the marked improvement in enantioselectivity in the inner-sphere (PHOX)Pd-catalyzed allylic alkylation of N-benzoyl lactam nucleophiles over their carbocyclic counterparts. We employ density functional theory calculations to aid in the interpretation of experimental results. Ultimately, we propose that the enhancement in enantioselectivity arises primarily from noncovalent interactions between the substrate and ligand rather than secondary substrate chelation, as previously hypothesized.

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.

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.

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH)4H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (∆E) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ∆E estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree-Fock included (0%-54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ∆E is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn-Sham orbitals gives results within 4-11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ∼5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

Molecular Geometries and Vibrational Contributions to Reaction Thermochemistry Are Surprisingly Insensitive to the Choice of Basis Sets

Calculation of molecular geometries and harmonic vibrational frequencies are pre-requisites for thermochemistry calculations. Contrary to conventional wisdom, this paper demonstrates that quantum chemical predictions of the thermochemistry of many gas and solution phase chemical reactions appear to be very insensitive to the choice of basis sets. For a large test set of 80 diverse organic and transition-metal-containing reactions, variations in reaction free energy based on geometries and frequencies calculated using a variety of double and triple-zeta basis sets from the Pople, Jensen, Ahlrichs, and Dunning families are typically less than 4 kJ mol^-1, especially when the quasiharmonic oscillator correction is applied to mitigate the effects of low-frequency modes. Our analysis indicates that for many organic molecules and their transition states, high-level revDSD-PBEP86-D4 and DLPNO-CCSD(T)/(aug-)cc-pVTZ single-point energies usually vary by less than 2 kJ mol^-1 on density functional theory geometries optimized using basis sets ranging from 6-31+G(d) to aug-pcseg-2 and aug-cc-pVTZ. In cases where these single-point energies vary significantly, indicating sensitivity of molecular geometries to the choice of basis set, there is often substantial cancellation of errors when the reaction energy or barrier is calculated. The study concludes that the choice of basis set for molecular geometry and frequencies, particularly those considered in this study, is not critical for the accuracy of thermochemistry calculations in the gas or solution phase.

High throughput computations of the effective removal of liquified gases by novel perchlorate hybrid material

The utilization of hybrid materials in separation technology, sorbents, direct air capture (DAC) technology, sensors, adsorbents, and chiral material recognition has increased in the past decade due to the recognized impact of atmospheric pollutants and hazardous industrial gases on climate change. A novel hybrid material, perchlorate hybrid (PClH), has been proposed in this study for the effective sensory detection and trapping of atmospheric pollutants and industrial hazardous gases. The study evaluated the structural properties, adsorption mechanism, electronic sensitivity, and topological analysis of PClH using highly accurate computational methods (M062X-D3BJ/def2-ccpVTZ and DSDPBEP86/def2-ccpVTZ). The computational analysis demonstrated that PClH has considerable adsorption energies and favorable interaction with CO2, NO2, SO2, COCl2, and H2S. PClH is more suitable for detecting liquefiable gases such as COCl2, CO2, and SO2, and can be easily recovered under ambient conditions. Developing such materials can contribute to reducing hazardous gases and pollutants in the atmosphere, leading to a cleaner and safer environment.

The entropic penalty for associative reactions and their physical treatment during routine computations

A systematic study of the entropic penalty for associative reactions is presented. It is shown that computed solution-phase Gibbs free energies typically overestimate entropic contributions. This entropic penalty for associative reactions in solution, i.e., if the number of particles decreases along the reaction coordinate (sum of stoichiometric numbers ), originates from the insufficient treatment of entropic effects by implicit solvent models. We propose an additive correction scheme to Gibbs free energies that is suitable for routine applications by non-expert users. This correction is based on Garza's formalism for the solution-phase entropy [A. J. Garza, J. Chem. Theory Comput., 2019, 15, 3204.] that is physically sound and embedded into an efficient black-box type algorithm. To critically evaluate the entropic penalty and its proposed treatment, we compiled an experimental benchmark set of 31 Δ_rG and 22 in 15 different solvents. Using a representative best-practice computational protocol (at wave function theory (WFT) based DLPNO-CCSD(T) and density functional theory (DFT) based revDSD-PBEP86-D4 level with an implicit solvent model), we determined a sizeable entropic penalty ranging from 2-11 kcal mol^-1. Using the correction scheme presented herein, the entropic penalty is corrected to the chemical accuracy of ≤1 kcal mol^-1 (WFT and DFT). The same applies to at the WFT level. Barriers at the DFT level are overestimated by 2 kcal mol^-1 (classic) and underestimated by 2 kcal mol^-1 (corrected). This effect is attributed to the finding that barriers computed at the DFT level are systematically 2-3 kcal mol^-1 lower than barriers obtained with WFT.

Mechanistic Study and Conceptual Chemical Reactivity Analysis of Hydroboration of Carbon Dioxide Catalyzed by a Manganese(I)-PNP-Pincer Complex

Designing efficient and selective catalysts for carbon dioxide reduction is an intensive research area in the recent literature on homogeneous catalysis. In this work, we study the catalytic activity of a newly reported Mn(I)-PNP-pincer catalyst with an embedded aromatic ring. First, we systematically examine its capability to yield different products and highlight the importance of ligand aromaticity and steric effects on metal-ligand cooperativity. We then further conceptually probe its reactivity with descriptors from both conceptual density functional theory and an information-theoretic approach, thereby proposing a novel partitioning of the reaction coordinate into three relevant regions. Our results show that the reactivity in these different regions is governed by different properties such as steric effects, electrophilicity/nucleophilicity, or aromaticity. We anticipate that this methodology, with the analytical tools employed in this study, can be generalized and extended to other catalytic systems and find applications in designing better catalysts.

Aromatic diglycosides from Sophora tonkinensis and a multi-step conformer filtering procedure for TDDFT calculation of flexible glycoside

Three previously undescribed aromatic diglycosides (1, 5, and 8) and six known analogs (2-4, 6, 7, and 9) were isolated from the roots and rhizomes of Sophora tonkinensis Gagnep. Their structures were elucidated by detailed spectroscopic analysis. The absolute configuration of compound 8 was determined by comparing the experimental and TDDFT calculated ECD spectra of 8 and aglycone 8a. Furthermore, a multistep conformer filtering procedure for TDDFT calculation of flexible glycoside was proposed, which afforded high accuracy with acceptable computing cost in determining the absolute configuration of glycosides using quantum calculated ECD.

Structural confirmation of synthetic cannabinoids in seized electronic cigarette oil: A combined mass spectrometric and computational study

RATIONALE

Synthetic cannabinoids are some of the most used and abused new psychoactive substances, because they can produce a stronger intense pleasure than natural cannabis. Most of the new synthetic cannabinoids are structurally similar to existing synthetic cannabinoids and can be obtained by modifying partial structures of the latter without changing their effects. Therefore, the derivatization rules and common fragmentation patterns of synthetic cannabinoids could be used for rapid screening and structural identification of them.

METHODS

The derivatization rules of synthetic cannabinoids are summarized, and the common fragmentation pattern of synthetic cannabinoids including three typical cleavage pathways was explored and extended in this work based on combined mass spectrometry (MS) and density functional theory studies. Five synthetic cannabinoids in electronic cigarette oil from a drug case were separated and characterized using gas chromatography with MS and liquid chromatography coupled to high-resolution quadrupole Orbitrap MS.

RESULTS

The structures of five synthetic cannabinoids in seized electronic cigarette oil were deduced from electron impact ion source (EI) MS and high-resolution electrospray ionization (ESI) MSⁿ data, along with the derivatization rules and common fragmentation pattern of synthetic cannabinoids. The proposed structures of these synthetic cannabinoids were further verified via reference substances. Computational study showed that selective cleavage of these compounds was mainly controlled by spin population in EI-MS, but a tunneling effect arose from proton transfer in ESI-MSⁿ detection, which has been rarely reported in previous works.

CONCLUSIONS

Our results showed that EI-MS was suitable for identifying synthetic cannabinoids with aromatic ketone structure, which could also be extended to adamantane linked group. Nevertheless, synthetic cannabinoids with carbamoyl linked group were better characterized by high-resolution ESI-MSⁿ compared to EI-MS. This study demonstrated a method with promising potential for rapid and reliable screening of synthetic cannabinoids in mixtures with enhanced detection throughput and operation simplicity.

A double-hybrid density functional based on good local physics with outstanding performance on the GMTKN55 database

In two recent papers [A. D. Becke, J. Chem. Phys. 156, 214101 (2022) and A. D. Becke, J. Chem. Phys. 157, 234102 (2022)], we compared two Kohn-Sham density functionals based on physical modeling and theory with the best density-functional power-series fits in the literature. The best error statistics reported to date for a hybrid functional on the general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) chemical database of Goerigk et al. [Phys. Chem. Chem. Phys. 19, 32184 (2017)] were obtained. In the present work, additional second-order perturbation-theory terms are considered. The result is a 12-parameter double-hybrid density functional with the lowest GMTKN55 WTMAD2 "weighted total mean absolute deviation" error (1.76 kcal/mol) yet seen for any hybrid or double-hybrid density-functional approximation. We call it "DH23."

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.

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

The C3H3NO family of isomers is relevant in astrochemistry, even though its members are still elusive in the interstellar medium. To identify the best candidate for astronomical detection within this family, we developed a new computational protocol based on the minimum-energy principle. This approach aims to identify the most stable isomer of the family and consists of three steps. The first step is an extensive investigation that characterizes the vast number of compounds having the C3H3NO chemical formula, employing density functional theory for this purpose. The second step is an energy refinement, which is used to select isomers and relies on coupled cluster theory. The last step is a structural improvement with a final energy refinement that provides improved energies and a large set of accurate spectroscopic parameters for all isomers lying within 30 kJ mol−1 above the most stable one. According to this protocol, vinylisocyanate is the most stable isomer, followed by oxazole, which is about 5 kJ mol−1 higher in energy. The other stable species are pyruvonitrile, cyanoacetaldehyde, and cyanovinylalcohol. For all of these species, new computed rotational and vibrational spectroscopic data are reported, which complement those already available in the literature or fill current gaps.

β-Boron Effect Enables Regioselective and Stereospecific Electrophilic Addition to Alkenes

Electrophilic addition to alkenes is a textbook-taught reaction, yet it is not always possible to control the regioselectivity of addition to unsymmetrical 1,2-disubstituted substrates. We report the observation and applications of the β-boron effect that accounts for high regioselectivity in electrophilic addition reactions to allylic MIDA (N-methyliminodiacetic acid) boronates. While the well-established β-silicon effect bears partial resemblance to the observed reactivity, the silyl group is typically lost during functionalization. In contrast, the boryl moiety is retained in the product when B(MIDA) is used as the nucleophilic stabilizer. Mechanistic studies elucidate the origin of this effect and demonstrate how σ(C-B) hyperconjugation helps stabilize the incipient carbocation. This transformation represents a rare example of the stereospecific hydrohalogenation of secondary allyl MIDA-boronates that proceeds in a syn-fashion.

Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy

Pair atomic density fitting (PADF) has been identified as a promising strategy to reduce the scaling with system size of quantum chemical methods for the calculation of the correlation energy like the direct random-phase approximation (RPA) or second-order Møller-Plesset perturbation theory (MP2). PADF can however introduce large errors in correlation energies as the two-electron interaction energy is not guaranteed to be bounded from below. This issue can be partially alleviated by using very large fit sets, but this comes at the price of reduced efficiency and having to deal with near-linear dependencies in the fit set. One posibility is to use global density fitting (DF), but in this work, we introduce an alternative methodology to overcome this problem that preserves the intrinsically favorable scaling of PADF. We first regularize the Fock matrix by projecting out parts of the basis set which gives rise to orbital products that are hard to describe by PADF. After having thus obtained a reliable self-consistent field solution, we then also apply this projector to the orbital coefficient matrix to improve the precision of PADF-MP2 and PADF-RPA. We systematically assess the accuracy of this new approach in a numerical atomic orbital framework using Slater type orbitals (STO) and correlation consistent Gaussian type basis sets up to quintuple-ζ quality for systems with more than 200 atoms. For the small and medium systems in the S66 database we show the maximum deviation of PADF-MP2 and PADF-RPA relative correlation energies to DF-MP2 and DF-RPA reference results to be 0.07 and 0.14 kcal/mol, respectively. When the new projector method is used, the errors only slightly increase for large molecules and also when moderately sized fit sets are used the resulting errors are well under control. Finally, we demonstrate the computational efficiency of our algorithm by calculating the interaction energies of large, non-covalently bound complexes with more than 1000 atoms and 20000 atomic orbitals at the RPA@PBE/CC-pVTZ level of theory.

Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities

Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.

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.

Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions

A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet ¹L_a and ¹L_b states of linear oligoacenes as well as the triplet ³L_a and ³L_b states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S₀-T₁ singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.

Matrix-isolation and cryosolution-VCD spectra of α-pinene as benchmark for anharmonic vibrational spectra calculations

The inclusion of anharmonicity in vibrational spectral analyis remains associated to small molecular systems with up to a dozen of atoms, with half a dozen of non-hydrogen atoms, typically thesize of propylene oxide. One may see two reasons for this: first of all, larger systems are often thought to be computationally too demanding (high computational costs) for a full anharmonic vibrational analysis. Second, the identification of resonances and their correction is often considered something only expert theoreticians could address because of the lack of unequivocal criteria. In this contribution, we illustrate that resonances can indeed become a complex problem, which can be handled almost transparently thanks to recent advances in vibrational perturbation theory (VPT2). The study also emphasizes the importance and the central role played by experiment in benchmarking novel theoretical approaches. In fact, we herein provide the currently highest resolution VCD spectra available for α- and β-pinene obtained under matrix-isolation conditions and in liquid Xenon as solvent. They are interpreted by VPT2 with novel tests for the identification of resonances. Hence, the study demonstrates the mutual stimulation of advances in both experimental techniques and computational models.

Probing the Size Limit of Dispersion Energy Donors with a Bifluorenylidene Balance: Magic Cyclohexyl

We report the synthesis of 14 2,2'-disubstituted 9,9'-bifluorenylidenes as molecular balances for the quantification of London dispersion interactions between various dispersion energy donors. For all balances, we measured ΔG_Z/E at 333 K using ¹H NMR in seven organic solvents. For various alkyl and aryl substituents, we generally observe a preference for the "folded" Z-isomer due to attractive London dispersion interactions. The cyclohexyl-substituted system shows the largest Z-preference in this study with ΔG_Z/E = -0.6 ± 0.05 kcal mol^-1 in all solvents, owing to the rotational freedom of cyclohexyl groups paired with their large polarizability that maximizes London dispersion interactions. On the other hand, rigid and sterically more demanding substituents like tert-butyl unexpectedly favor the unfolded E-isomer. This is a result of the close relative position in which the functional groups are positioned in this molecular balance. This close proximity is the reason for the increase of Pauli repulsion in the Z-isomers with large rigid substituents (tert-butyl, adamantyl, and diamantyl) which leads to an equilibrium shift toward the unfolded E-form. While we were able to reproduce most of our experimental trends qualitatively using contemporary computational chemistry methods, quantitative accuracy of the employed methods still needs further improvement.

A comprehensive benchmark investigation of quantum chemical methods for carbocations

The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol^-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol^-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol^-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.

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.

Synthesis of Sterically Encumbered Thiourea S-Oxides through Direct Thiourea Oxidation

Thiourea S-oxides can be viewed as formal analogs of the currently unknown diamino-substituted Criegee intermediates (urea O-oxides). However, the preparation of such S-oxides is rather challenging, and the direct oxidation of thioureas typically only leads to formation of desulfurized products. Employing the accurate revDSD-PBEP86-D4 double hybrid density functional, it was found that the peracid mediated oxidation of thiourea S-oxides exhibits a lower reaction barrier than the oxidation of the corresponding thiourea itself in contrast to most other ordinary thioketones. The undesired overoxidation reactivity, which is associated with strong π-donation from the thiourea's nitrogen atoms, can be partially suppressed by introduction of bulky substituents and the utilization of protic solvents. In this regard, we managed to prepare two sterically encumbered thiourea S-oxides in isolated yields of 35-40 %. The S-oxides are stable in the solid state and in alcoholic solutions at room temperature for extended periods of time, but swiftly decompose in aprotic solvents by disproportionation. A dimesityl-substituted thiourea S-oxide complexed with residual mCBA could be characterized by means of X-ray crystallography, confirming the importance of hydrogen bonding in the stabilization of the amino-substituted C=S⁺ -O^- moiety.

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.

Building on the strengths of a double-hybrid density functional for excitation energies and inverted singlet-triplet energy gaps

It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet-triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over "safe" excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet-triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature-the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S₁ and T₁ excited states.

Performance of Localized-Orbital Coupled-Cluster Approaches for the Conformational Energies of Longer n-Alkane Chains

We report an update and enhancement of the ACONFL (conformer energies of large alkanes [J. Phys. Chem. A2022,126, 3521-3535]) dataset. For the ACONF12 (n-dodecane) subset, we report basis set limit canonical coupled-cluster with singles, doubles, and perturbative triples [i.e., CCSD(T)] reference data obtained from the MP2-F12/cc-pV{T,Q}Z-F12 extrapolation, [CCSD(F12*)-MP2-F12]/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/aug-cc-pV{D,T}Z calculations. Then, we explored the performance of a variety of single and composite localized-orbital CCSD(T) approximations, ultimately finding an affordable localized natural orbital CCSD(T) [LNO-CCSD(T)]-based post-MP2 correction that agrees to 0.006 kcal/mol mean absolute deviation with the revised canonical reference data. In tandem with canonical MP2-F12 complete basis set extrapolation, this was then used to re-evaluate the ACONF16 and ACONF20 subsets for n-hexadecane and n-icosane, respectively. Combining those with the revised canonical reference data for the dodecane conformers (i.e., ACONF12 subset), a revised ACONFL set was obtained. It was then used to assess the performance of different localized-orbital coupled-cluster approaches, such as pair natural orbital localized CCSD(T) [PNO-LCCSD(T)] as implemented in MOLPRO, DLPNO-CCSD(T₀) and DLPNO-CCSD(T₁) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, at their respective "Normal", "Tight", "vTight", and "vvTight" accuracy settings. For a given accuracy threshold and basis set, DLPNO-CCSD(T₁) and DLPNO-CCSD(T₀) perform comparably. With "VeryTightPNO" cutoffs, explicitly correlated DLPNO-CCSD(T₁)-F12/VDZ-F12 is the best pick among all the DLPNO-based methods tested. To isolate basis set incompleteness from localized-orbital-related truncation errors (domain, LNOs), we have also compared the localized coupled-cluster approaches with canonical DF-CCSD(T)/aug-cc-pVTZ for the ACONF12 set. We found that gradually tightening the cutoffs improves the performance of LNO-CCSD(T), and using a composite scheme such as vTight + 0.50[vTight - Tight] improves things further. For DLPNO-CCSD(T₁), "TightPNO" and "VeryTightPNO" offer a statistically similar accuracy, which gets slightly better when T_CutPNO is extrapolated to the complete PNO space limit. Similar to Brauer et al.'s [Phys. Chem. Chem. Phys.2016,18 (31), 20905-20925] previous report for the S66x8 noncovalent interactions, the dispersion-corrected direct random phase approximation (dRPA)-based double hybrids perform remarkably well for the ACONFL set. While the revised reference data do not affect any conclusions on the less accurate methods, they may upend orderings for more accurate methods with error statistics on the same order as the difference between reference datasets.

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.

Reliability of Computing van der Waals Bond Lengths of Some Rare Gas Diatomics

When the bond lengths of 11 molecules containing van der Waals bonds are optimized by 572 methods and 20 basis sets, it is found that the best mean absolute deviations (MADs) of density-functional theory (DFT) methods are 0.005 Å (shown by APFD/6-311++G**), 0.007 Å (B2PLYPD3(Full)/aug-cc-pVQZ), and 0.010 Å (revDSDPBEP86/aug-cc-pVQZ), while the best MADs of ab initio methods are 0.008 Å (BD(T)/aug-cc-pVTZ) and 0.016 Å (MP4/aug-cc-pVQZ). Moreover, the best MADs calculated by 54 selected methods in combination with 60 other basis sets (such as 6-311++G, 6-31++G(3d'f,3p'd), and UGBS1V++) are not better. Therefore, these bond lengths can be calculated with extremely high accuracy by some special methods and basis sets, and CCSD(T) is also not as good as expected because its best MAD is only 0.023 Å (CCSD(T)/aug-cc-pVQZ).

Halogen Atom Participation in Guiding the Stereochemical Outcomes of Acetal Substitution Reactions

Acetal substitution reactions of α-halogenated five- and six-membered rings can be highly stereoselective. Erosion of stereoselectivity occurs as nucleophilicity increases, which is consistent with additions to a halogen-stabilized oxocarbenium ion, not a three-membered-ring halonium ion. Computational investigations confirmed that the open-form oxocarbenium ions are the reactive intermediates involved. Kinetic studies suggest that hyperconjugative effects and through-space electrostatic interactions can both contribute to the stabilization of halogen-substituted oxocarbenium ions.

Best-Practice DFT Protocols for Basic Molecular Computational Chemistry

Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.

Reduced-Scaling Double Hybrid Density Functional Theory with Rapid Basis Set Convergence through Localized Pair Natural Orbital F12

Following earlier work [Mehta, N.; Martin, J. M. L. J. Chem. Theory Comput.2022, 10.1021/acs.jctc.2c00426] that showed how the slow basis set convergence of the double hybrid density functional theory can be obviated by the use of F12 explicit correlation in the GLPT2 step (second order Görling-Levy perturbation theory), we demonstrate here for the very large and chemically diverse GMTKN55 benchmark suite that the CPU time scaling of this step can be reduced (asymptotically linearized) using the localized pair natural orbital (PNO-L) approximation at negligible cost in accuracy.

Silyl Groups Are Strong Dispersion Energy Donors

We present an experimental and computational study to investigate noncovalent interactions between silyl groups that are often employed as "innocent" protecting groups. We chose an extended cyclooctatetraene (COT)-based molecular balance comprising unfolded (1,4-disubstituted) and folded (1,6-disubstituted) valance bond isomers that typically display remote and close silyl group contacts, respectively. The thermodynamic equilibria were determined using nuclear magnetic resonance measurements. Additionally, we utilized Boltzmann weighted symmetry-adapted perturbation theory (SAPT) at the sSAPT0/aug-cc-pVDZ level of theory to dissect and quantify noncovalent interactions. Apart from the extremely bulky tris(trimethylsilyl)silyl "supersilyl" group, there is a preference for the folded 1,6-COT valence isomer, with London dispersion interactions being the main stabilizing factor. This makes silyl groups excellent dispersion energy donors, a finding that needs to be taken into account in synthesis planning.

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.

Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies

With increasing interelectronic distance, the screening of the electron–electron interaction by the presence of other electrons becomes the dominant source of electron correlation. This effect is described by the random phase approximation (RPA) which is therefore a promising method for the calculation of weak interactions. The success of the RPA relies on the cancellation of errors, which can be traced back to the violation of the crossing symmetry of the 4-point vertex, leading to strongly overestimated total correlation energies. By the addition of second-order screened exchange (SOSEX) to the correlation energy, this issue is substantially reduced. In the adiabatic connection (AC) SOSEX formalism, one of the two electron–electron interaction lines in the second-order exchange term is dynamically screened (SOSEX(W, v_c)). A related SOSEX expression in which both electron–electron interaction lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, v_c), the evaluation of this correlation energy expression does not require an expensive numerical frequency integration and is therefore advantageous from a computational perspective. We compare the accuracy of the statically screened variant to RPA and RPA+SOSEX(W, v_c) for a wide range of chemical reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, v_c) for charged excitations and noncovalent interactions where they lead to major improvements over RPA.

Anharmonic Aspects in Vibrational Circular Dichroism Spectra from 900 to 9000 cm-1 for Methyloxirane and Methylthiirane

Vibrational circular dichroism (VCD) spectra and the corresponding IR spectra of the chiral isomers of methyloxirane and of methylthiirane have been reinvestigated, both experimentally and theoretically, with particular attention to accounting for anharmonic corrections, as calculated by the GVPT2 approach. De novo recorded VCD spectra in the near IR (NIR) range regarding CH-stretching overtone transitions, together with the corresponding NIR absorption spectra, were also considered and accounted for, both with the GVPT2 and with the local mode approaches. Comparison of the two methods has permitted us to better describe the nature of active “anharmonic” modes in the two molecules and the role of mechanical and electrical anharmonicity in determining the intensities of VCD and IR/NIR data. Finally, two nonstandard IR/NIR regions have been investigated: the first one about ≈2000 cm^–1, involving mostly two-quanta bending mode transitions, the second one between 7000 and 7500 cm^–1 involving three-quanta transitions containing CH-stretching overtones and HCC/HCH bending modes.

Explicitly Correlated Double-Hybrid DFT: A Comprehensive Analysis of the Basis Set Convergence on the GMTKN55 Database

Double-hybrid density functional theory (DHDFT) offers a pathway to accuracy approaching composite wavefunction approaches such as G4 theory. However, the Görling–Levy second-order perturbation theory (GLPT2) term causes them to partially inherit the slow ∝L^–3 (with L the maximum angular momentum) basis set convergence of correlated wavefunction methods. This could potentially be remedied by introducing F12 explicit correlation: we investigate the basis set convergence of both DHDFT and DHDFT-F12 (where GLPT2 is replaced by GLPT2-F12) for the large and chemically diverse general main-group thermochemistry, kinetics, and noncovalent interactions (GMTKN55) benchmark suite. The B2GP-PLYP-D3(BJ) and revDSD-PBEP86-D4 DHDFs are investigated as test cases, together with orbital basis sets as large as aug-cc-pV5Z and F12 basis sets as large as cc-pVQZ-F12. We show that F12 greatly accelerates basis set convergence of DHDFs, to the point that even the modest cc-pVDZ-F12 basis set is closer to the basis set limit than cc-pV(Q+d)Z or def2-QZVPPD in orbital-based approaches, and in fact comparable in quality to cc-pV(5+d)Z. Somewhat surprisingly, aug-cc-pVDZ-F12 is not required even for the anionic subsets. In conclusion, DHDF-F12/VDZ-F12 eliminates concerns about basis set convergence in both the development and applications of double-hybrid functionals. Mass storage and I/O bottlenecks for larger systems can be circumvented by localized pair natural orbital approximations, which also exhibit much gentler system size scaling.