Seemingly Simple Stereoelectronic Effects in Alkane Isomers and the Implications for Kohn–Sham Density Functional Theory

Carbon-Fluorine Activation in the Gas Phase: The Reactions of Benzyl C-F Bonds and Silyl Cations

The activation of C-F bonds figures largely in both fundamental and applied chemical processes. Herein the activation of benzyl C-F bonds by silyl cations is examined both computationally and experimentally in the gas phase. The experimental rate constant values obtained herein have not heretofore been measured and provide insight into the intrinsic ability of silyl cations to activate C-F bonds. Trends in reactivity and correlations between theoretical and experimental data are discussed in the context of C-F bond cleavage.

How do spin-scaled double hybrids designed for excitation energies perform for noncovalent excited-state interactions? An investigation on aromatic excimer models

Time-dependent double hybrids with spin-component or spin-opposite scaling to their second-order perturbative correlation correction have demonstrated competitive robustness in the computation of electronic excitation energies. Some of the most robust are those recently published by our group (M. Casanova-Páez, L. Goerigk, J. Chem. Theory Comput. 2021, 20, 5165). So far, the implementation of these functionals has not allowed correctly calculating their ground-state total energies. Herein, we define their correct spin-scaled ground-state energy expressions which enables us to test our methods on the noncovalent excited-state interaction energies of four aromatic excimers. A range of 22 double hybrids with and without spin scaling are compared to the reasonably accurate wavefunction reference from our previous work (A. C. Hancock, L. Goerigk, RSC Adv. 2023, 13, 35964). The impact of spin scaling is highly dependent on the underlying functional expression, however, the smallest overall errors belong to spin-scaled functionals with range separation: SCS- and SOS- ω PBEPP86, and SCS-RSX-QIDH. We additionally determine parameters for DFT-D3(BJ)/D4 ground-state dispersion corrections of these functionals, which reduce errors in most cases. We highlight the necessity of dispersion corrections for even the most robust TD-DFT methods but also point out that ground-state based corrections are insufficient to completely capture dispersion effects for excited-state interaction energies.

Big data benchmarking: how do DFT methods across the rungs of Jacob's ladder perform for a dataset of 122k CCSD(T) total atomization energies?

Total atomization energies (TAEs) are a central quantity in density functional theory (DFT) benchmark studies. However, so far TAE databases obtained from experiment or high-level ab initio wavefunction theory included up to hundreds of TAEs. Here, we use the GDB-9 database of 133k CCSD(T) TAEs generated by Curtiss and co-workers [B. Narayanan, P. C. Redfern, R. S. Assary and L. A. Curtiss, Chem. Sci., 2019, 10, 7449] to evaluate the performance of 14 representative DFT methods across the rungs of Jacob's ladder (namely, PBE, BLYP, B97-D, M06-L, τ-HCTH, PBE0, B3LYP, B3PW91, ωB97X-D, τ-HCTHh, PW6B95, M06, M06-2X, and MN15). We first use the A25[PBE] diagnostic for nondynamical correlation to eliminate systems that potentially include significant multireference effects, for which the CCSD(T) TAEs might not be sufficiently reliable. The resulting database (denoted by GDB9-nonMR) includes 122k species. Of the considered functionals, B3LYP attains the best performance relative to the G4(MP2) reference TAEs, with a mean absolute deviation (MAD) of 4.09 kcal mol-1. This first-generation hybrid functional, in which the three mixing coefficients were fitted against a small set of TAEs, is one of the few functionals that are not systematically biased towards overestimating the G4(MP2) TAEs, as demonstrated by a mean-signed deviation (MSD) of 0.45 kcal mol-1. The relatively good performance of B3LYP is followed by the heavily parameterized M06-L meta-GGA functional, which attains a MAD of 6.24 kcal mol-1. The PW6B95, M06, M06-2X, and MN15 functionals tend to systematically overestimate the G4(MP2) TAEs and attain MADs ranging between 18.69 (M06) and 28.54 (MN15) kcal mol-1. However, PW6B95 and M06-2X exhibit particularly narrow error distributions. Thus, scaling their TAEs by an empirical scaling factor reduces their MADs to merely 3.38 (PW6B95) and 2.85 (M06-2X) kcal mol-1. Empirical dispersion corrections (e.g., D3 and D4) are attractive, and therefore, their inclusion worsens the performance of methods that systematically overestimate the TAEs.

Influence of Steric and Dispersion Interactions on the Thermochemistry of Crowded (Fluoro)alkyl Compounds

ConspectusAlkanes play a pivotal role in industrial, environmental, and biological processes. They are characterized by their carbon-carbon single-bond structure, remarkable stability, and conformational diversity. Fluorination of such compounds imparts unique physicochemical properties that often enhance pharmacokinetic profiles, metabolic stability, and receptor interactions while keeping beneficial properties. However, such per- and polyfluoroalkyl substances (PFAS) show a persistent presence in the environment and potential adverse health effects, which propelled them to the forefront of global environmental and health discussions. Alkyl compounds are also prototypical for stereoelectronic (SE) effects that are widely applied in chemistry. Substituents are typically described as electron-density-donating/withdrawing and/or responsible for sterically interacting with reagents or strategic groups in the molecule. That alkane branching can result in higher stability compared to less-branched isomers has been investigated in detail also by testing quantum chemical methods, in particular density functional theory (DFT). Alkane branching results in spatially compact structures with close intramolecular contacts so that at a specific size the detailed balance of attractive London dispersion and covalent versus repulsive Pauli exchange interactions shifts to new, chemically unfragile situations. This may lead to dissociation at room temperature and opens the central question: what is the smallest crowed alkane that cannot be made synthetically? In this Account, we try to shed light on the interplay among the various (free) energy components for crowded (fluoro)alkane dissociation. In this context, homolytic cleavage of the central C-C bond in a series of model alkanes of increasing size with tert-butyl (tBu), adamantyl (Ad), and [1.1.1]propellanyl (Prop) substituents is investigated. Reference energies are calculated at the PNO-LCCSD(T)-F12b level and used to benchmark the performance of contemporary DFT functionals. In line with previous conclusions, the application of dispersion corrections to density functionals is mandatory. For crowed structures, the accurate description of the midrange correlation effects, specifically repulsive van der Waals interactions, is crucial, and we observed that the density-dependent VV10 correction is superior to D4 in this context, although the asymptotic region is better described by the latter. The best available dispersion-inclusive functionals show systematic and reasonably small residual errors and can be safely applied to large systems (>100 atoms), for which coupled cluster methods with large basis sets are not computationally feasible anymore. For qualitatively correct predictions of synthetic accessibility under equilibrium conditions (free energy), the inclusion of thermostatistical (entropy) contributions is also essential. According to our results, tetra-tert-butylmethane (C17tBu) is the largest and most crowded system with a positive dissociation free energy and should be synthesizable. The difference between hydrogenated and perfluorinated systems originates from the increase in the steric repulsion of spatially close substituents, which is not compensated to the same extent by attractive orbital and dispersion interactions. A sometimes-assumed similar steric demand for fluorine and hydrogen atoms is not corroborated by our investigations on crowded systems. Perfluorination is found to substantially decrease thermal stability, rendering perfluorinated hexamethylethane (C8tBuF) the last potentially stable representative.

Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost

We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.

An Artificial Intelligence Approach for Tackling Conformational Energy Uncertainties in Chiroptical Spectroscopies

Determination of the absolute configuration of chiral molecules is a prerequisite for obtaining a fundamental understanding in any chirality-related field. The interaction with polarised light has proven to be a powerful means to determine this absolute configuration, but its application rests on the comparison between experimental and computed spectra for which the inherent uncertainty in conformational Boltzmann factors has proven to be extremely hard to tackle. Here we present a novel approach that overcomes this issue by combining a genetic algorithm that identifies the relevant conformers by accounting for the uncertainties in DFT relative energies, and a hierarchical clustering algorithm that analyses the trends in the spectra of the considered conformers and identifies on-the-fly when a given chiroptical technique is not able to make reliable predictions. The effectiveness of this approach is demonstrated by considering the challenging cases of papuamine and haliclonadiamine, two bis-indane natural products with eight chiral centres and considerable conformational heterogeneity that could not be assigned unambiguously with current approaches.

Simplified and enhanced VCD analysis of cyclic peptides guided by artificial intelligence

Cyclic peptides are privileged structures in medicinal chemistry; however, their solution-state structure characterization is difficult. Vibrational circular dichroism (VCD) spectroscopy is a powerful alternative to NMR, but requires challenging calculations. We present a VCD approach guided by a genetic algorithm, which is simple, more effective, and has a higher conformer resolution.

Interpretable Graph-Network-Based Machine Learning Models via Molecular Fragmentation

Chemists have long benefitted from the ability to understand and interpret the predictions of computational models. With the current shift to more complex deep learning models, in many situations that utility is lost. In this work, we expand on our previously work on computational thermochemistry and propose an interpretable graph network, FragGraph(nodes), that provides decomposed predictions into fragment-wise contributions. We demonstrate the usefulness of our model in predicting a correction to density functional theory (DFT)-calculated atomization energies using Δ-learning. Our model predicts G4(MP2)-quality thermochemistry with an accuracy of <1 kJ mol^-1 for the GDB9 dataset. Besides the high accuracy of our predictions, we observe trends in the fragment corrections which quantitatively describe the deficiencies of B3LYP. Node-wise predictions significantly outperform our previous model predictions from a global state vector. This effect is most pronounced as we explore the generality by predicting on more diverse test sets indicating node-wise predictions are less sensitive to extending machine learning models to larger molecules.

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.

An Overview of Basis Set Effects for Diatomic Boron Nitride Compounds (B2N(∓,0)): A Quantum Symmetry Breaking

The symmetry breaking (SB) of B2 not only exhibits an energy barrier for ionic or neutral forms dependent on various basis sets but it also exhibits a few SBs due to the asymmetry stretching and bending mode interactions. SB obeys the mechanical quantum theorem among discrete symmetries and their connection to the spin statistics in physical sciences. In this investigation, the unusual amount of energy barrier of SBs appeared upon the orbit–orbit coupling of BNB (both radical and ions) between transition states and the ground state. Our goal in this study is to understand the difference among the electromagnetic structures of the (B2N(∓,0)) variants due to effects of various basis sets and methods and also the quantum symmetry breaking phenomenon. In the D∞h point group of (B2N(∓,0)) variants, the unpaired electron is delocalized, while in the asymmetric C∞v point group, it is localized on either one of the B atoms. Structures with broken symmetry, C∞v, can be stable by interacting with the D∞h point group. In viewpoints of quantum chemistry, the second-order Jahn–Teller effect permits the unpaired electron to localize on boron atom, rather than being delocalized. In this study, we observed that the energy barrier of SB for BNB increases by post HF methods.

Interstellar hide and go seek: C3H4O. There and back (again)

The molecular species C₃H₄O represents a striking example of an astrochemical conundrum. With more than 60 structural isomers theoretically possible, to date only acrolein (CH₂CHCHO) has been identified in the Sgr B2(N) region of the interstellar medium (ISM). The topography of the singlet potential energy surface is complicated, with three low-lying minima predicted to be almost isoenergetic: cis and trans-acrolein, and methylketene (CH₃CHCO). Our CCSD(T)/cc-pVTZ calculations confirm that methylketene is energetically lower than cis-acrolein, lying only 1.9 kJ mol^-1 above the trans-isomer, which is the global minimum. In this respect, methylketene is a promising candidate for interstellar observation. Unfortunately, however, despite several searches its astronomical detection has been unsuccessful. To this end, the key question is whether in fact methylketene exists as a discrete chemical entity in the ISM at all? In this paper, we present a detailed examination of the C₃H₄O potential energy surface, with specific focus on formation pathways. CCSD(T)/cc-pVTZ calculations enable a more elaborate interpretation of reaction mechanisms than was published hitherto. Our results show that gauche-propargyl alcohol and syn and anti-allenol emerge as interesting new targets for observational astronomers in TMC-1: given the recent discovery of the propargyl radical in this region, barrierless product channels involving OH˙ lend support to their candidacy as possible interstellar species. Finally, this work provides accurate spectral data of these three potential molecules, to be used for searches in interstellar space.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.

BH9, a New Comprehensive Benchmark Data Set for Barrier Heights and Reaction Energies: Assessment of Density Functional Approximations and Basis Set Incompleteness Potentials

The calculation of accurate reaction energies and barrier heights is essential in computational studies of reaction mechanisms and thermochemistry. To assess methods regarding their ability to predict these two properties, high-quality benchmark sets are required that comprise a reasonably large and diverse set of organic reactions. Due to the time-consuming nature of both locating transition states and computing accurate reference energies for reactions involving large molecules, previous benchmark sets have been limited in scope, the number of reactions considered, and the size of the reactant and product molecules. Recent advances in coupled-cluster theory, in particular local correlation methods like DLPNO-CCSD(T), now allow the calculation of reaction energies and barrier heights for relatively large systems. In this work, we present a comprehensive and diverse benchmark set of barrier heights and reaction energies based on DLPNO-CCSD(T)/CBS called BH9. BH9 comprises 449 chemical reactions belonging to nine types common in organic chemistry and biochemistry. We examine the accuracy of DLPNO-CCSD(T) vis-a-vis canonical CCSD(T) for a subset of BH9 and conclude that, although there is a penalty in using the DLPNO approximation, the reference data are accurate enough to serve as a benchmark for density functional theory (DFT) methods. We then present two applications of the BH9 set. First, we examine the performance of several density functional approximations commonly used in thermochemical and mechanistic studies. Second, we assess our basis set incompleteness potentials regarding their ability to mitigate basis set incompleteness errors. The number of data points, the diversity of the reactions considered, and the relatively large size of the reactant molecules make BH9 the most comprehensive thermochemical benchmark set to date and a useful tool for the development and assessment of computational methods.

Strong π-stacking causes unusually large anisotropic thermal expansion and thermochromism

π-stacking in ground-state dimers/trimers/tetramers of N-butoxyphenyl(naphthalene)diimide (BNDI) exceeds 50 kcal ⋅ mol^-1 in strength, drastically surpassing that for the ^*3[pyrene]₂ excimer (∼30 kcal ⋅ mol^-1; formal bond order = 1) and similar to other weak-to-moderate classical covalent bonds. Cooperative π-stacking in triclinic (BNDI-T) and monoclinic (BNDI-M) polymorphs effects unusually large linear thermal expansion coefficients (α _a , α _b , α _c , β) of (452, -16.8, -154, 273) × 10^-6 ⋅ K^-1 and (70.1, -44.7, 163, 177) × 10^-6 ⋅ K^-1, respectively. BNDI-T exhibits highly reversible thermochromism over a 300-K range, manifest by color changes from orange (ambient temperature) toward red (cryogenic temperatures) or yellow (375 K), with repeated thermal cycling sustained for over at least 2 y.

CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions?

We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.

Explaining and Fixing DFT Failures for Torsional Barriers

Most torsional barriers are predicted with high accuracies (about 1 kJ/mol) by standard semilocal functionals, but a small subset was found to have much larger errors. We created a database of almost 300 carbon-carbon torsional barriers, including 12 poorly behaved barriers, that stem from the Y═C-X group, where Y is O or S and X is a halide. Functionals with enhanced exchange mixing (about 50%) worked well for all barriers. We found that poor actors have delocalization errors caused by hyperconjugation. These problematic calculations are density-sensitive (i.e., DFT predictions change noticeably with the density), and using HF densities (HF-DFT) fixes these issues. For example, conventional B3LYP performs as accurately as exchange-enhanced functionals if the HF density is used. For long-chain conjugated molecules, HF-DFT can be much better than exchange-enhanced functionals. We suggest that HF-PBE0 has the best overall performance.

Do Docking Sites Persist Upon Fluorination? The Diadamantyl Ether-Aromatics Challenge for Rotational Spectroscopy and Theory

Fluorinated derivatives of biological molecules have proven to be highly efficient at modifying the biological activity of a given protein through changes in the stability and the kind of docking interactions. These interactions can be hindered or facilitated based on the hydrophilic/hydrophobic character of a particular protein region. Diadamantyl ether (C₂₀H₃₀O) possesses both kinds of docking sites, serving as a good template to model these important contacts with aromatic fluorinated counterparts. In this work, an experimental study on the structures of several complexes between diadamantyl ether and benzene as well as a series of fluorinated benzenes is reported to analyze the effect of H→F substitution on the interaction and structure of the resulting molecular clusters using rotational spectroscopy. All experimentally observed complexes are largely dominated by London dispersion interactions with the hydrogen‐terminated surface areas of diadamantyl ether. Already single substitution of one hydrogen atom with fluorine changes the preferred docking site of the complexes. However, the overall contributions of the different intermolecular interactions are similar for the different complexes, contrary to previous studies focusing on the difference in interactions using fluorinated and non‐fluorinated molecules.

Assessing the Applicability of the Geometric Counterpoise Correction in B2PLYP/Double-ζ Calculations for Thermochemistry, Kinetics, and Noncovalent Interactions

We present a proof-of-concept study of the suitability of Kruse and Grimme’s geometric counterpoise correction (gCP) for basis set superposition errors (BSSEs) in double-hybrid density functional calculations with a double-ζ basis set. The gCP approach only requires geometrical information as an input and no orbital/density information is needed. Therefore, this correction is practically free of any additional cost. gCP is trained against the Boys and Bernardi counterpoise correction across a set of 528 noncovalently bound dimers. We investigate the suitability of the approach for the B2PLYP/def2-SVP level of theory, and reveal error compensation effects—missing London dispersion and the BSSE—associated with B2PLYP/def2-SVP calculations, and present B2PLYP-gCP-D3(BJ)/def2-SVP with the reparametrised DFT-D3(BJ) and gCP corrections as a more balanced alternative. Benchmarking results on the S66x8 benchmark set for noncovalent interactions and the GMTKN55 database for main-group thermochemistry, kinetics, and noncovalent interactions show a statistical improvement of the B2PLYP-gCP-D3(BJ) scheme over plain B2PLYP and B2PLYP-D3(BJ). B2PLYP-D3(BJ) shows significant overestimation of interaction energies, barrier heights with larger deviations from the reference values, and wrong relative stabilities in conformers, all of which can be associated with BSSE. We find that the gCP-corrected method represents a significant improvement over B2PLYP-D3(BJ), particularly for intramolecular noncovalent interactions. These findings encourage future developments of efficient double-hybrid DFT strategies that can be applied when double-hybrid calculations with large basis sets are not feasible due to system size.

Extrapolation in quantum chemistry: Insights on energetics and reaction dynamics

Since there is no exact solution for problems in physics and chemistry, extrapolation methods may assume a key role in quantitative quantum chemistry. Two topics where it bears considerable impact are addressed, both at the heart of computational quantum chemistry: electronic structure and reaction dynamics. In the first, the problem of extrapolating the energy obtained by solving the electronic Schrödinger equation to the limit of the complete one-electron basis set is addressed. With the uniform-singlet-and-triplet-extrapolation (USTE) scheme at the focal point, the emphasis is on recent updates covering from the energy itself to other molecular properties. The second topic refers to extrapolation of quantum mechanical reactive scattering probabilities from zero total angular momentum to any of the values that it may assume when running quasiclassical trajectories, QCT/QM-[Formula: see text]J. With the extrapolation guided in both cases by physically motivated asymptotic theories, realism is seeked by avoiding unsecure jumps into the unknown. Although, mostly review oriented, a few issues are addressed for the first time here and there. Prospects for future work conclude the overview.

Protobranching as repulsion-induced attraction: a prototype for geminal stabilization

Noncovalent interactions are traditionally defined within the context of their attractive components, such as electrostatics and dispersion. Sources of molecular strain are derived through the destabilization of Coulombic and exchange repulsion. Due to this binary designation, the underlying origin of geminal stability with respect to alkanes (referred to as protobranching) has been an active subject for debate between these competing perspectives. We recast this stabilization as a complementary (Gestalt) interaction between dispersion and exchange repulsion, each impacting the other. We use triplet hydrogen and argon dimer as foundational van der Waals adducts to develop a procedure for the visualization and quantification of both exchange repulsion, Δρ_SCF, and medium-range correlation, ΔΔρ, as perturbations in electron density. We use the framework of the DFT-D3 correction to reproduce the shape of the dispersion potential at medium range and successfully model the trend in stability for the eighteen isomers of octane with a diverse series of functionals: BLYP, B3LYP, BP86, PBE, and PBE0. Collectively, our findings show that protobranching is a manifestation of steric repulsion-reduction in vibrational enthalpy and medium-range electron correlation.

Difficulties of Popular Density Functionals to Describe the Conformational Isomerism in Iodoacetic Acid

Matrix isolation studies in solid argon and neon at 4.2 K reveal that iodoacetic acid initially only exists as its ground state (c,x) conformer with an almost perpendicular I-C-C═O dihedral angle, but UV irradiation in the 240-255 nm range leads to population of the 0.8 kcal mol^-1 less stable (c,c) isomer. The latter structure exhibits a close 3.23 Å contact of the iodine and carbonyl oxygen atoms decidedly below the sum of their van der Waals radii (3.50 Å). Increasing the matrix temperature by only a few Kelvin triggers the thermal back reaction of (c,c) to (c,x) and leads to an estimated upper limit of 0.38 kcal mol^-1 for the associated torsional barrier. While wave function methods including completely uncorrelated Hartree-Fock theory have no problem to identify (c,c) as a proper minimum, many popular density functionals fail to describe the C-C torsional potential in cis-iodoacetic acid qualitatively correct. We assessed the performance of 12 density functionals of different levels of sophistication, namely, the BLYP, PBE, TPSS, B3LYP, BHandHLYP, PBE0, M06-2X, CAM-B3LYP, ωB97X-D3, B2-PLYP, B2GP-PLYP, and DSD-PBEP86 methods, against accurate extrapolated CCSD(T)/CBS(T-Q)//MP2/def2-TZVPP energies and found that almost all of them yield acceptable relative energies. Still, even some of the best performers fail to find a reasonably deep minimum in the region of the (c,c) conformer, and addition of the empirical D3-dispersion correction does not remedy the qualitative shortcoming. Instead, inclusion of a sufficient amount of (long-range) exact exchange and likely a proper treatment of medium-range correlation effects all along the torsional coordinate play an important role in the proper description of the sub-van der Waals iodine-oxygen contact. More modern, recommended functionals do not suffer from the described shortcoming.

Double hybrid DFT calculations with Slater type orbitals

On a comprehensive database with 1,644 datapoints, covering several aspects of main-group as well as of transition metal chemistry, we assess the performance of 60 density functional approximations (DFA), among them 36 double hybrids (DH). All calculations are performed using a Slater type orbital (STO) basis set of triple-ζ (TZ) quality and the highly efficient pair atomic resolution of the identity approach for the exchange- and Coulomb-term of the KS matrix (PARI-K and PARI-J, respectively) and for the evaluation of the MP2 energy correction (PARI-MP2). Employing the quadratic scaling SOS-AO-PARI-MP2 algorithm, DHs based on the spin-opposite-scaled (SOS) MP2 approximation are benchmarked against a database of large molecules. We evaluate the accuracy of STO/PARI calculations for B3LYP as well as for the DH B2GP-PLYP and show that the combined basis set and PARI-error is comparable to the one obtained using the well-known def2-TZVPP Gaussian-type basis set in conjunction with global density fitting. While quadruple-ζ (QZ) calculations are currently not feasible for PARI-MP2 due to numerical issues, we show that, on the TZ level, Jacob's ladder for classifying DFAs is reproduced. However, while the best DHs are more accurate than the best hybrids, the improvements are less pronounced than the ones commonly found on the QZ level. For conformers of organic molecules and noncovalent interactions where very high accuracy is required for qualitatively correct results, DHs provide only small improvements over hybrids, while they still excel in thermochemistry, kinetics, transition metal chemistry and the description of strained organic systems.

Thermal stabilities and conformational behaviors of isocyanurates and cyclotrimerization energies of isocyanates: a computational study

Isocyanurates are cyclic trimers of isocyanate molecules. They are generally known as highly thermostable compounds. However, it is interesting how the thermal stabilities of the isocyanurate molecules will be altered depending on the substituents of their three nitrogen atoms. We performed computational investigations on the thermochemical behaviors of isocyanurate molecules with various alkyl and phenyl substituents. The cyclotrimerization processes of isocyanates are highly exothermic. Our best estimate of the enthalpy change for the cyclotrimerization of methyl isocyanate into trimethyl isocyanurate was -66.4 kcal mol-1. Additional negative cyclotrimerization enthalpy changes were observed for n-alkyl-substituted isocyanates. This trend was enhanced with an extension of n-alkyl chains. Conversely, low negative cyclotrimerization enthalpy changes were shown for secondary and tertiary alkyl-substituted isocyanates. The n-alkyl-substituted isocyanurates were shown to be stabilized due to attractive dispersion interactions between the substituents. Meanwhile, the branched alkyl-substituted isocyanurates were destabilized due to the deformation of their isocyanurate rings. For various alkyl-substituted isocyanates, the sum of the deformation energy of the isocyanurate ring and the intramolecular inter-substituent nonbonding interaction energies was found to be linearly correlated with their cyclotrimerization energies. The cyclotrimerization energy for phenyl isocyanate was shown to have significantly deviated from the linear relationship observed for the alkyl-substituted isocyanurates. This is probably attributable to a remarkable change in the orbital resonance interactions during the cyclotrimerization of phenyl isocyanate to triphenyl isocyanurate.

Determination of the absolute configuration of conformationally flexible molecules by simulation of chiro-optical spectra: a case study

The assignment of the absolute configuration (AC) of two conformational flexible organic molecules by means of TD-DFT simulation of the electronic circular dichroism (ECD) spectra is presented. The factors leading to a reliable assignment were evaluated in the various steps of the process. The effects of different functionals and basis sets in the geometry optimization step is very limited in terms of the resulting optimized geometries, whereas the inclusion of the solvent in the calculations has a much larger effect on the correct evaluation of the conformational ratio. B3LYP and M06-2x were found to be the most accurate functionals for geometry optimization. CAM-B3LYP and ωB97X-D provided the best results in the TD-DFT simulations.

Evaluating Transition Metal Barrier Heights with the Latest Density Functional Theory Exchange-Correlation Functionals: The MOBH35 Benchmark Database

A new database of transition metal reaction barrier heights (MOBH35) is presented. Benchmark energies (forward and reverse barriers and reaction energy) are calculated using DLPNO-CCSD(T) extrapolated to the complete basis set limit using a Weizmann-1-like scheme. Using these benchmark energies, the performance of a wide selection of density functional theory (DFT) exchange-correlation functionals, including the latest from the Martin, Truhlar, and Head-Gordon groups, is evaluated. It was found, using the def2-TZVPP basis set, that the ωB97M-V (MAD 1.7 kcal/mol), ωB97M-D3BJ (MAD 1.9 kcal/mol), ωB97X-V (MAD 2.0 kcal/mol), and revTPSS0-D4 (MAD 2.2 kcal/mol) hybrid functionals are recommended. The double-hybrid functionals B2K-PLYP (MAD 1.7 kcal/mol) and revDOD-PBEP86-D4 (MAD 1.8 kcal/mol) also performed well, but this has to be balanced by their increased computational cost.

Semiempirical Quantum-Chemical Methods with Orthogonalization and Dispersion Corrections

We present two new semiempirical quantum-chemical methods with orthogonalization and dispersion corrections: ODM2 and ODM3 (ODM x). They employ the same electronic structure model as the OM2 and OM3 (OM x) methods, respectively. In addition, they include Grimme's dispersion correction D3 with Becke-Johnson damping and three-body corrections E ABC for Axilrod-Teller-Muto dispersion interactions as integral parts. Heats of formation are determined by adding explicitly computed zero-point vibrational energy and thermal corrections, in contrast to standard MNDO-type and OM x methods. We report ODM x parameters for hydrogen, carbon, nitrogen, oxygen, and fluorine that are optimized with regard to a wide range of carefully chosen state-of-the-art reference data. Extensive benchmarks show that the ODM x methods generally perform better than the available MNDO-type and OM x methods for ground-state and excited-state properties, while they describe noncovalent interactions with similar accuracy as OM x methods with a posteriori dispersion corrections.

Effect of the exchange–correlation functional on the synchronicity/nonsynchronicity in bond formation in Diels–Alder reactions: a reaction force constant analysis

The performance of 24 KS-DFT-based methods (GGA, MGGA, HGGA, HMGGA, and DHGGA) was assessed, finding that M11 and M06-2X (HMGGA) predicting reliable TS geometries, energetics, and (a)synchronicities in Diels–Alder reactions.

Stable dicationic dioxoliums and fate of their dioxolyl radicals

A glimpse into uncharted territory: the synthesis and study of dicationic dioxolium salts allow for assessing the fate of the corresponding elusive dioxolyl radicals.

Thermochemistry of Guanine Tautomers Re-Examined by Means of High-Level CCSD(T) Composite Ab Initio Methods

We obtained accurate gas-phase tautomerization energies for a set of 14 guanine tautomers by means of high-level thermochemical procedures approximating the CCSD(T) energy at the complete basis set (CBS) limit. For the five low-lying tautomers, we use the computationally demanding W1-F12 composite method for obtaining the tautomerization energies. The relative W1-F12 tautomerization enthalpies at 298K are: 0.00 (1), 2.37 (2), 2.63 (3), 4.03 (3′), and 14.31 (4) kJmol−1. Thus, as many as four tautomers are found within a small energy window of less than 1.0kcalmol−1 (1kcalmol−1=4.184kJmol−1). We use these highly accurate W1-F12 tautomerization energies to evaluate the performance of a wide range of lower-level composite ab initio procedures. The Gn composite procedures (G4, G4(MP2), G4(MP2)-6X, G3, G3B3, G3(MP2), and G3(MP2)B3) predict that the enol tautomer (3) is more stable than the keto tautomer (2) by amounts ranging from 0.36 (G4) to 1.28 (G3(MP2)) kJmol−1. We also find that an approximated CCSD(T)/CBS energy calculated as HF/jul-cc-pV{D,T}Z+CCSD/jul-cc-pVTZ+(T)/jul-cc-pVDZ results in a root-mean-square deviation (RMSD) of merely 0.11kJmol−1 relative to the W1-F12 reference values. We use this approximated CCSD(T)/CBS method to obtain the tautomerization energies of 14 guanine tautomers. The relative tautomerization enthalpies at 298K are: 0.00 (1), 2.20 (2), 2.51 (3), 4.06 (3′), 14.30 (4), 25.65 (5), 43.78 (4′), 53.50 (6′), 61.58 (6), 77.37 (7), 82.52 (8′), 86.02 (9), 100.70 (10), and 121.01 (8) kJmol−1. Using these tautomerization enthalpies, we evaluate the performance of standard and composite methods for the entire set of 14 guanine tautomers. The best-performing procedures emerge as (RMSDs are given in parentheses): G4(MP2)-6X (0.51), CCSD(T)+ΔMP2/CBS (0.52), and G4(MP2) (0.64kJmol−1). The worst performers are CCSD(T)/AVDZ (1.05), CBS-QB3 (1.24), and CBS-APNO (1.38kJmol−1).

Semi-empirical or non-empirical double-hybrid density functionals: which are more robust?

The development of non-empirical double-hybrid density functionals (DHDFs) is a very active research area with the number of approaches in this field having increased rapidly. At the same time, there is a lack of published work that provides a fair assessment and comparison between non-empirical and semi-empirical DHDFs on an equal footing. Herein, we close this gap and present a thorough analysis of both classes of DHDFs on the large GMTKN55 benchmark database for general main-group thermochemistry, kinetics, and noncovalent interactions [Goerigk et al., Phys. Chem. Chem. Phys., 2017, 19, 32184-32215]. In total, 115 variations of dispersion-corrected and -uncorrected DHDFs are tested, which will be condensed to an in-depth assessment of 31 methods: 19 non-empirical and 12 semi-empirical DHDFs. As such, our study represents the largest DHDF study ever conducted and can serve as an important benchmark informing method developers and users alike. Our results show that semi-empirical DHDFs are the most robust density functional approximations and more reliable and accurate than non-empirical ones. In fact, some non-empirical approaches are even outperformed by hybrid approaches or even dispersion-corrected and -uncorrected MP2 and SCS-MP2. SOS0-PBE0-2-D3(BJ) is the only exception and the only non-empirical DHDF that we can safely recommend for general applicability. However, it is still outperformed by six semi-empirical DHDFs, of which we would like to particularly recommend the following five: ωB97X-2-D3(BJ), DSD-BLYP-D3(BJ), DSD-PBEP86-D3(BJ), B2NC-PLYP-D3(BJ), and B2GPPLYP-D3(BJ). Our findings seriously question current trends in the field and they highlight that novel strategies have to be found in order to outperform the currently best density functional theory methods on the market. We hope that our study can function as an important cornerstone inspiring such a change of direction in the field.

Adaptive responses of sterically confined intramolecular chalcogen bonds

The existence of intramolecular chalcogen bonds (IChBs) in 2,6-disubstituted arylchalcogen derivatives is determined by the substituents and the sigma hole donor behavior of the chalcogen atom in the molecule.

The responsive behavior of an entity towards its immediate surrounding is referred to as an adaptive response. The adaptive responses of a noncovalent interaction at the molecular scale are reflected from its structural and functional roles. Intramolecular chalcogen bonding (IChB), an attractive interaction between a heavy chalcogen E (E = Se or Te) centered sigma hole and an ortho-heteroatom Lewis base donor D (D = O or N), plays an adaptive role in defining the structure and reactivity of arylchalcogen compounds. In this perspective, we describe the adaptive roles of a chalcogen centered Lewis acid sigma hole and a proximal Lewis base (O or N) in accommodating built-in steric stress in 2,6-disubstituted arylchalcogen compounds. From our perspective, the IChB components (a sigma hole and the proximal Lewis base) act in synergism to accommodate the overwhelming steric force. The adaptive responses of the IChB components are inferred from the observed molecular structures and reactivity. These include (a) adaptation of a conformation without IChBs, (b) adaptation of a conformation with weak IChBs, (c) twisting the skeletal aryl ring while maintaining IChBs, (d) ionization of the E–X bond (e.g., X = Br) to relieve stress and (e) intramolecular cyclization to relieve steric stress. A comprehensive approach, involving X-ray data analysis, density functional theory (DFT) calculations, reaction pattern analysis and principal component analysis (PCA), has been employed to rationalize the adaptive behaviors of IChBs in arylchalcogen compounds. We believe that the perception of ChB as an adaptive/stimulus responsive interaction would profit the futuristic approaches that would utilise ChB as self-assembly and molecular recognition tools.

Perspective: Chemical Information Encoded in Electron Density

In this perspective, we review the chemical information encoded in electron density and other ingredients used in semilocal functionals. This information is usually looked at from the functional point of view: the exchange density or the enhancement factor are discussed in terms of the reduced density gradient. However, what parts of a molecule do these 3D functions represent? We look at these quantities in real space, aiming to understand the electronic structure information they encode and provide an insight from the quantum chemical topology (QCT). Generalized gradient approximations (GGAs) provide information about the presence of chemical interactions, whereas meta-GGAs can differentiate between the different bonding types. By merging these two techniques, we show new insight into the failures of semilocal functionals owing to three main errors: fractional charges, fractional spins, and non-covalent interactions. We build on simple models. We also analyze the delocalization error in hydrogen chains, showing the ability of QCT to reveal the delocalization error introduced by semilocal functionals. Then, we show how the analysis of localization can help understand the fractional spin error in alkali atoms, and how it can be used to correct it. Finally, we show that the poor description of GGAs of isodesmic reactions in alkanes is due to 1,3-interactions.

CBS extrapolation in electronic structure pushed to the end: a revival of minimal and sub-minimal basis sets

The complete basis set (CBS) limit is secluded in calculations of electronic structure, and hence CBS extrapolation draws immediate attention.

‘Diet GMTKN55’ offers accelerated benchmarking through a representative subset approach

The GMTKN55 benchmarking protocol allows comprehensive analysis and ranking of density functional approximations with diverse chemical behaviours. This work reports diet versions of GMTKN55 which reproduce key properties of the full protocol at substantially reduced numerical cost. ‘Diet GMTKN55’ can thus be used for benchmarking expensive methods, or in combination with solid state benchmarks.

Multilevel Quantum Mechanics Theories and Molecular Mechanics Calculations of the Cl- + CH3I Reaction in Water

The Cl^- + CH₃I → CH₃Cl + I^- reaction in water was studied using combined multilevel quantum mechanism theories and molecular mechanics with an explicit water solvent model. The study shows a significant influence of aqueous solution on the structures of the stationary points along the reaction pathway. A detailed, atomic-level evolution of the reaction mechanism shows a concerted one-bond-broken and one-bond-formed mechanism, as well as a synchronized charge-transfer process. The potentials of mean force calculated with the CCSD(T) and DFT treatments of the solute produce a free activation barrier at 24.5 and 19.0 kcal/mol, respectively, which agrees with the experimental one at 22.0 kcal/mol. The solvent effects have also been quantitatively analyzed: in total, the solvent effects raise the activation energy by 20.2 kcal/mol, which shows a significant impact on this reaction in water.

Benchmarking density functional tight binding models for barrier heights and reaction energetics of organic molecules

Density Functional Tight Binding (DFTB) models are two to three orders of magnitude faster than ab initio and Density Functional Theory (DFT) methods and therefore are particularly attractive in applications to large molecules and condensed phase systems. To establish the applicability of DFTB models to general chemical reactions, we conduct benchmark calculations for barrier heights and reaction energetics of organic molecules using existing databases and several new ones compiled in this study. Structures for the transition states and stable species have been fully optimized at the DFTB level, making it possible to characterize the reliability of DFTB models in a more thorough fashion compared to conducting single point energy calculations as done in previous benchmark studies. The encouraging results for the diverse sets of reactions studied here suggest that DFTB models, especially the most recent third-order version (DFTB3/3OB augmented with dispersion correction), in most cases provide satisfactory description of organic chemical reactions with accuracy almost comparable to popular DFT methods with large basis sets, although larger errors are also seen for certain cases. Therefore, DFTB models can be effective for mechanistic analysis (e.g., transition state search) of large (bio)molecules, especially when coupled with single point energy calculations at higher levels of theory. © 2017 Wiley Periodicals, Inc.