Pericyclic reaction benchmarks: hierarchical computations targeting CCSDT(Q)/CBS and analysis of DFT performance

Radical Addition Reactions: Hierarchical Ab Initio Benchmark and DFT Performance Study

We performed a hierarchical ab initio benchmark study of the gas-phase radical addition reactions of X⋅+C2H2 and X⋅+C2H4 (X⋅ = CH3⋅, NH2⋅, OH⋅, SH⋅). The hierarchical series of ab initio methods (HF, MP2, CCSD, CCSD(T)) were paired with a hierarchal series of Dunning basis sets with and without diffuse functions ((aug)-cc-pVDZ, (aug)-cc-pVTZ, (aug)-cc-pVQZ). The HF ground-state wavefunctions were transformed into quasi-restricted orbital (QRO) reference wavefunctions to address spin contamination. Following extrapolation to the CBS limit, the energies from our highest- QRO-CCSD(T)/CBS+ level converged within 0.0-3.4 kcal mol-1 and 0.0-1.0 kcal mol-1 concerning the ab initio method and basis set, respectively. Our QRO-CCSD(T)/CBS+ reference data was used to evaluate the performance of 98 density functional theory (DFT) approximations. The MAE of the best functionals for reaction barriers and energies were: OLYP (1.9 kcal mol-1), BMK (1.0 kcal mol-1), M06-2X (0.9 kcal mol-1), MN12-SX (0.8 kcal mol-1) and CAM-B3LYP (0.8 kcal mol-1). These functionals also accurately reproduce key geometrical parameters of the stationary points within an average 2 % deviation from the reference QRO-CCSD(T)/cc-pVTZ level.

Assessing Computational Methods to Calculate the Binding Energies of Dimers of Five-Membered Heterocyclic Molecules

Computational electronic structure methods, including ab initio and density functional theory (DFT), have been assessed in calculating the binding energies of 14 five-membered heterocyclic dimers. The configurations were generated using classical molecular dynamics before optimization at the MP2/aug-cc-pVTZ. Benchmark binding energies are calculated at the CCSD(T)/CBS level of theory. Among the ab initio methods, the DLPNO-CCSD(T)/CBS method has the best performance, reproducing CCSD(T)/CBS with a mean absolute deviation (MAD) of 0.17 kcal/mol. In addition, a schematic CCSD(T)/CBS approach perfectly reproduces the canonical CCSD(T)/CBS with a mean absolute error of 0.08 kcal/mol. Regarding DFT functionals, it has been found that counterpoise corrections have negligible effects on the accuracy of the functionals. Furthermore, including the D3 empirical dispersion considerably enhances the accuracy of the DFT functionals. As a result, outstanding performance is noted for the double hybrid functional B2K-PLYP, with a mean absolute error of 0.25 kcal/mol. In addition to the B2K-PLYP double hybrid functional, M05-D3, B97D, M05-2X-D3, M05-2X, M06-HF, M08-HX, M11, TPSSh-D3, and RSX-0DH-D3(BJ) have MAD values lower than 0.5 kcal/mol. These functionals are recommended for further investigations of five-membered heterocyclic clusters.

Computational Insights into Tunable Reversible Network Materials: Accelerated ReaxFF Kinetics of Furan-Maleimide Diels-Alder Reactions for Self-Healing and Recyclability

In this study, ReaxFF molecular dynamics simulations were benchmarked and used to study the relative kinetics of the retro Diels-Alder reaction between furan and N-methylmaleimide. This reaction is very important for the creation of polymer networks with self-healing and recyclable properties, since they can be used as reversible linkers in the network. So far, the reversible Diels-Alder reaction has not yet been studied by using reactive molecular dynamics simulations. This work is, thus, the first step in simulating a covalent adaptable network (CAN) using Diels-Alder reactions as reversible linkers. For both endo and exo, the bond breaking in 40 product molecules was simulated using the bond boost method and the endo/exo ratio was evaluated. This ratio was benchmarked against density functional theory (DFT) and experimental results for a changing set of bond boost parameters. Given their importance to understand how the CAN performs, the effect of the addition of a polymer backbone and the effect of temperature were successfully simulated using our newly parametrized reactive force field.

Fulminic acid: a quasibent spectacle

Fulminic acid (HCNO) played a critical role in the early development of organic chemistry, and chemists have sought to discern the structure and characteristics of this molecule and its isomers for over 200 years. The mercurial nature of the extremely flat H-C-N bending potential of fulminic acid, with a nearly vanishing harmonic vibrational frequency at linearity, remains enigmatic and refractory to electronic structure theory, as dramatic variation with both orbital basis set and electron correlation method is witnessed. To solve this problem using rigorous electronic wavefunction methods, we have employed focal point analyses (FPA) to ascertain the ab initio limit of optimized linear and bent geometries, corresponding vibrational frequencies, and the HCN + O(3P) → HCNO reaction energy. Electron correlation treatments as extensive as CCSDT(Q), CCSDTQ(P), and even CCSDTQP(H) were employed, and complete basis set (CBS) extrapolations were performed using the cc-pCVXZ (X = 4-6) basis sets. Core electron correlation, scalar relativistic effects (MVD1), and diagonal Born-Oppenheimer corrections (DBOC) were all included and found to contribute significantly in determining whether vibrationless HCNO is linear or bent. At the all-electron (AE) CCSD(T)/CBS level, HCNO is a linear molecule with ω5(H-C-N bend) = 120 cm-1. However, composite AE-CCSDT(Q)/CBS computations give an imaginary frequency (51i cm-1) at the linear optimized geometry, leading to an equilibrium structure with an H-C-N angle of 173.9°. Finally, at the AE-CCSDTQ(P)/CBS level, HCNO is once again linear with ω5 = 45 cm-1, and inclusion of both MVD1 and DBOC effects yields ω5 = 32 cm-1. A host of other topics has also been investigated for fulminic acid, including the dependence of re and ωi predictions on a variety of CBS extrapolation formulas, the question of multireference character, the N-O bond energy and enthalpy of formation, and issues that give rise to the quasibent appellation.

Retro-Cope elimination of cyclic alkynes: reactivity trends and rational design of next-generation bioorthogonal reagents

The retro-Cope elimination reaction between dimethylhydroxylamine (DMHA) and various cyclic alkynes has been quantum chemically explored using DFT at ZORA-BP86/TZ2P. The purpose of this study is to understand the role of the following three unique activation modes on the overall reactivity, that is (i) additional cycloalkyne predistortion via fused cycles, (ii) exocyclic heteroatom substitution on the cycloalkyne, and (iii) endocyclic heteroatom substitution on the cycloalkyne. Trends in reactivity are analyzed and explained by using the activation strain model (ASM) of chemical reactivity. Based on our newly formulated design principles, we constructed a priori a suite of novel bioorthogonal reagents that are highly reactive towards the retro-Cope elimination reaction with DMHA. Our findings offer valuable insights into the design principles for highly reactive bioorthogonal reagents in chemical synthesis.

Controlling Diastereoselectivity in Dearomatizing Diels-Alder Reactions of Nitroarenes with 2-Trimethylsilyloxycyclohexadiene

Dearomative Diels-Alder cycloadditions between nitroarenes and 2-trimethylsilyloxycyclohexadiene are carried out under high pressure at room temperature in the absence of any chemical promoter. Reactions are performed with different arenes, including the highly aromatic naphthalenes and quinolines. They lead to 3D-scaffolds with exquisite exo-diastereoselectivity. The exo approach is characterized by lower distortion of the substrates in a late TS and by more favorable orbital interactions presumably between the nitro group and the dienic part, explaining the stereoselectivity.

Explainable chemical artificial intelligence from accurate machine learning of real-space chemical descriptors

Machine-learned computational chemistry has led to a paradoxical situation in which molecular properties can be accurately predicted, but they are difficult to interpret. Explainable AI (XAI) tools can be used to analyze complex models, but they are highly dependent on the AI technique and the origin of the reference data. Alternatively, interpretable real-space tools can be employed directly, but they are often expensive to compute. To address this dilemma between explainability and accuracy, we developed SchNet4AIM, a SchNet-based architecture capable of dealing with local one-body (atomic) and two-body (interatomic) descriptors. The performance of SchNet4AIM is tested by predicting a wide collection of real-space quantities ranging from atomic charges and delocalization indices to pairwise interaction energies. The accuracy and speed of SchNet4AIM breaks the bottleneck that has prevented the use of real-space chemical descriptors in complex systems. We show that the group delocalization indices, arising from our physically rigorous atomistic predictions, provide reliable indicators of supramolecular binding events, thus contributing to the development of Explainable Chemical Artificial Intelligence (XCAI) models.

Quantum Tunneling Instability in Pericyclic Reactions

Several cycloreversion reactions of the retro-Diels-Alder type were computationally assessed to understand their quantum tunneling (QT) reactivity. N2, CO, and other leaving groups were considered based on their strong exothermicity, as it reduces their thermodynamic and kinetic stabilities. Our results indicate that for many of these reactions, it is essential to take into account their QT decomposition rate, which can massively weaken their molecular stability and shorten their half-lives even at deep cryogenic temperatures. In practical terms, this indicates that many supposedly stable molecules will actually be unsynthesizable or unisolable, and therefore trying to prepare or detect them would be a futile attempt. In addition, we discuss the importance of tunneling to correctly understand the enthalpy of activation and the collective atomic effect on the tunneling kinetic isotope effects to test if third-row atoms can tunnel in a chemical reaction. This project raises the question of the importance of in silico chemistry to guide in vitro chemistry, especially in cases where the latter cannot solve its own uncertainties.

Enantiodetermining processes in the synthesis of alanine, serine, and isovaline

In this study, quantum chemical calculations were used to explore the synthesis of three chiral α-amino acids, specifically alanine, serine, and isovaline, from reactants found in interstellar space. Our focus is on the crucial step in the synthesis pathway that determines the chirality of the amino acids. The results indicate that in the case of alanine, the determination of enantiomer is primarily influenced by the direction of the collision of molecules or functional groups, which leads to the formation of a chirality center in a crucial intermediate. However, contrary to chemical expectations, the enantiodetermining/enantioselection step for serine and isovaline synthesis occurs prior to the creation of a chirality center.

(3+2)-Cycloadditions of Levoglucosenone (LGO) with Fluorinated Nitrile Imines Derived from Trifluoroacetonitrile: An Experimental and Computational Study

The in situ-generated N-aryl nitrile imines derived from trifluoroacetonitrile smoothly undergo (3+2)-cycloadditions onto the enone fragment of the levoglucosenone molecule, yielding the corresponding, five-membered cycloadducts. In contrast to the 'classic' C(Ph),N(Ph) nitrile imine, reactions with fluorinated C(CF3),N(Ar) analogues lead to stable pyrazolines in a chemo- and stereoselective manner. Based on the result of X-ray single crystal diffraction analysis, their structures were established as exo-cycloadducts with the location of the N-Ar terminus of the 1,3-dipole at the α-position of the enone moiety. The DFT computation demonstrated that the observed reaction pathway results from the strong dominance of kinetic control over thermodynamic control.

Net Intermolecular Silyloxypyrone-Based (5+2) Cycloadditions Utilizing Amides as Enabling and Cleavable Tethers

Silyloxypyrone-based (5+2) cycloadditions were facilitated by amides that allowed for increased reactivity and a pathway for cleaving the tether to afford net intermolecular cycloadducts. Various amides underwent facile cycloaddition, and several experiments revealed steric and electronic factors that accelerate the reaction. tert-Butyl amides reacted faster than less hindered variants in multiple cases. In the case of dearomative oxidopyrylium-indole (5+2) cycloadditions, an amine-based tether was ineffective, whereas amides enabled this powerful transformation. Theoretical calculations evidenced a concerted asynchronous reaction in which the amide facilitates a conformational driving force enabling cycloaddition. Finally, a one-pot acylation/(5+2) cycloaddition/nucleophilic lactam opening and other examples of tosyl lactam opening of a modified cycloadduct were demonstrated.

Not antiaromaticity gain, but increased asynchronicity enhances the Diels-Alder reactivity of tropone

Tropone is an unreactive diene in normal electron demand Diels-Alder reactions, but it can be activated via carbonyl umpolung by using hydrazone ion analogs. Recently, the higher reactivity of hydrazone ion analogs was ascribed to a raised HOMO energy induced by antiaromaticity (L. J. Karas, A. T. Campbell, I. V. Alabugin and J. I. Wu, Org. Lett., 2020, 22, 7083). We show that this is incorrect, and that the activation barrier is lowered by increased asynchronicity.

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.

Reaction profiles for quantum chemistry-computed [3 + 2] cycloaddition reactions

Bio-orthogonal click chemistry based on [3 + 2] dipolar cycloadditions has had a profound impact on the field of biochemistry and significant effort has been devoted to identify promising new candidate reactions for this purpose. To gauge whether a prospective reaction could be a suitable bio-orthogonal click reaction, information about both on- and off-target activation and reaction energies is highly valuable. Here, we use an automated workflow, based on the autodE program, to compute over 5000 reaction profiles for [3 + 2] cycloadditions involving both synthetic dipolarophiles and a set of biologically-inspired structural motifs. Based on a succinct benchmarking study, the B3LYP-D3(BJ)/def2-TZVP//B3LYP-D3(BJ)/def2-SVP level of theory was selected for the DFT calculations, and standard conditions and an (aqueous) SMD model were imposed to mimic physiological conditions. We believe that this data, as well as the presented workflow for high-throughput reaction profile computation, will be useful to screen for new bio-orthogonal reactions, as well as for the development of novel machine learning models for the prediction of chemical reactivity more broadly.

How Bases Catalyze Diels-Alder Reactions

We have quantum chemically studied the base-catalyzed Diels-Alder (DA) reaction between 3-hydroxy-2-pyrone and N-methylmaleimide using dispersion-corrected density functional theory. The uncatalyzed reaction is slow and is preceded by the extrusion of CO2 via a retro-DA reaction. Base catalysis, for example, by triethylamine, lowers the reaction barrier up to 10 kcal mol-1 , causing the reaction to proceed smoothly at low temperature, which quenches the expulsion of CO2 , yielding efficient access to polyoxygenated natural compounds. Our activation strain analyses reveal that the base accelerates the DA reaction via two distinct electronic mechanisms: i) by the HOMO-raising effect, which enhances the normal electron demand orbital interaction; and ii) by donating charge into 3-hydroxy-2-pyrone which accumulates in its reactive region and promotes strongly stabilizing secondary electrostatic interactions with N-methylmaleimide.

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.

Automated Variable Electric-Field DFT Application for Evaluation of Optimally Oriented Electric Fields on Chemical Reactivity

Recent theoretical work and experiments at molecular junctions have provided a strong conceptualization for the effects of oriented electric fields (OEFs) on organic reactions. Depending on the axis of application, OEFs can increase (or decrease) the reaction rate or distinguish between isomeric pathways. Despite the conceptual elegance of OEFs, which may be applied externally or induced locally, as tools for catalyzing organic reactions, implementation in synthetically relevant systems has been hampered by inefficiencies in evaluating reaction sensitivity to field effects. Herein, we describe the development of the Automated Variable Electric-Field DFT Application (A.V.E.D.A.) for streamlined evaluation of a reaction's susceptibility to OEFs. This open-source software was designed to be accessible for nonexpert users of computational and programming tools. Following initiation by a single command (and with no subsequent intervention) the Linux workflow manages a series of density functional theory calculations and mathematical manipulations to optimize local-minimum and transition-state structures in oriented electric fields of increasing magnitude. The resulting molecular and reaction dipole moments, field-perturbed geometries, and net effective activation energies are compiled for user interpretation. Ten representative pericyclic reactions that showcase the development and evaluation of A.V.E.D.A. are described.

Mechanism and the Origins of Periselectivity in Cycloaddition Reactions of Benzyne with Dienes

Density functional theory calculations have been used to explore the reaction mechanism of (4 + 2) and (2 + 2) cycloadditions of benzyne with classical dienes. The results indicate the following: (1) (4 + 2) products arise via concerted pathways, (2) (2 + 2) products arise via stepwise pathways with diradical intermediates, and (3) these diradical intermediates are formed via isomerization of carbene intermediates. The origins of periselectivity in these reactions are analyzed using distortion/interaction analysis for the key steps, and they indicate that the tiny distortion in the very early [4 + 2] transition structure, coupled with an entropic favorability, controls selective (4 + 2) cycloaddition.