Analysis of transition state stabilization by non-covalent interactions in the Houk-List model of organocatalyzed intermolecular Aldol additions using functional-group symmetry-adapted perturbation theory

An efficient method by combining different basis sets and SAPT levels

In symmetry-adapted perturbation theory (SAPT), accurate calculations on non-covalent interaction (NCI) for large complexes with more than 50 atoms are time-consuming using large basis sets. More efficient ones with smaller basis sets usually result in poor prediction in terms of dispersion and overall energies. In this study, we propose two composite methods with baseline calculated at SAPT2/aug-cc-pVDZ and SAPT2/aug-cc-pVTZ with dispersion term corrected at SAPT2+ level using bond functions and smaller basis set with δ MP2 corrections respectively. Benchmark results on representative NCI data sets, such as S22, S66, and so forth, show significant improvements on the accuracy compared to the original SAPT Silver standard and comparable to SAPT Gold standard in some cases with much less computational cost.

Nanoscale solvent organization in metal-organic framework ZIF-8 probed by EPR of flexible β-phosphorylated nitroxides

Metal-organic frameworks (MOFs) draw increasing attention as nanoenvironments for chemical reactions, especially in the field of catalysis. Knowing the specifics of MOF cavities is decisive in many of these cases; yet, obtaining them in situ remains very challenging. We report the first direct assessment of the apparent polarity and solvent organization inside MOF cavities using a dedicated structurally flexible spin probe. A stable β-phosphorylated nitroxide radical was incorporated into the cavities of a prospective MOF ZIF-8 in trace amounts. The spectroscopic properties of this probe depend on local polarity, structuredness, stiffness and cohesive pressure and can be precisely monitored by Electron Paramagnetic Resonance (EPR) spectroscopy. Using this approach, we have demonstrated experimentally that the cavities of bare ZIF-8 are sensed by guest molecules as highly non-polar inside. When various alcohols fill the cavities, remarkable self-organization of solvent molecules is observed leading to a higher apparent polarity in MOFs compared to the corresponding bulk alcohols. Accounting for such nanoorganization phenomena can be crucial for optimization of chemical reactions in MOFs, and the proposed methodology provides unique routes to study MOF cavities inside in situ, thus aiding in their various applications.

Effect of Solvents on Proline Modified at the Secondary Sphere: A Multivariate Exploration

The critical influence of solvent effects on proline-catalyzed aldol reactions has been extensively described. Herein, we apply multivariate regression strategies to probe the influence of different solvents on an aldol reaction catalyzed by proline modified at its secondary sphere with boronic acids. In this system, both in situ binding of the boronic acid to proline and the outcome of the aldol reaction are impacted by the solvent-controlled microenvironment. Thus, with the aim of uncovering mechanistic insight and an ancillary aim of identifying methodological improvements, we designed a set of experiments, spanning 15 boronic acids in five different solvents. Based on hypothesized intermediates or interactions that could be responsible for the selectivity in these reactions, we proposed several structural configurations for the library of boronic acids. Subsequently, we compared the statistical models correlating the outcome of the reaction in different solvents with molecular descriptors produced for each of these proposed configurations. The models allude to the importance of different interactions in controlling selectivity in each of the studied solvents. As a proof-of-concept for the practicality of our approach, the models in chloroform ultimately led to lowering the ketone loading to only two equivalents while retaining excellent yield and enantio- and diastereo-selectivity.

Modulation of π character upon complexation captured by molecular rotation spectra

Two configurations of the furan–CF3Cl complex have been observed by high-resolution rotational spectroscopy. One is characterized by a dominant Cl lone pairs∙∙∙π*aromatic interaction and the other is stabilized by a...

On the question of steric repulsion versus noncovalent attractive interactions in chiral phosphoric acid catalyzed asymmetric reactions

The origin of enantioselectivity in asymmetric catalysis is often built around the differential steric interaction in the enantiocontrolling transition states (TSs). A closer perusal of enantiocontrolling TSs in an increasingly diverse range of reactions has revealed that the cumulative effect of weak noncovalent interactions could even outweigh the steric effects. While enunciating this balance is conspicuously important, quantification of such intramolecular forces within a TS continues to remain scarce and challenging. Herein, we demonstrate the utility of the fragment molecular orbital method in establishing the relative contributions of various attractive and repulsive contributions in the total interaction energy between the suitably chosen fragments in enantiocontrolling TSs. Three types of reactions of high contemporary importance, namely, axially chiral phosphoric acid (CPA) catalyzed kinetic resolution of rac-α-methyl-γ-hydroxy ester (reaction I), asymmetric dearomative amination of β-naphthols by dimethyl azodicarboxylate (IIa and IIb), and intramolecular desymmetrization of β,β-disubstituted methyl oxetanes (IIIa) and hydroxyl oxetane (IIIb), bearing a tethered alcohol (-OCH₂CH₂OH or -(CH₂)₂CH₂OH), are considered. In all the five reactions, the differences in the stabilizing contributions arising due to electrostatic, charge-transfer, and dispersion interactions between the catalyst and the reacting partners in the enantiocontrolling transition states are weighed against the destabilizing exchange interaction. The balancing interactions are found to be between dispersion and exchange repulsion in reaction I, a combination of charge transfer and dispersion energies offsets the repulsive energy in reaction IIb involving the electron rich anthryl groups in the catalyst, whereas the -(CF₃)₂C₆H₄ 3,3'-substituent in the catalyst (reaction IIa) leads to a trade-off between dispersion and exchange energies. In reactions IIIa and IIIb, however, electrostatic and dispersion energies help compensate the repulsive interactions. These quantitative insights on the intramolecular interactions in the stereocontrolling TSs could help in the rational design of asymmetric catalysis.

Isomerization of Functionalized Olefins Using the Dinuclear Catalyst [PdI(μ-Br)(PtBu3)]2: A Mechanistic Study

In a combined experimental and computational study, the isomerization activity of the dinuclear palladium(I) complex [Pd^I(μ‐Br)(P^tBu₃)]₂ towards allyl arenes, esters, amides, ethers, and alcohols has been investigated. The calculated energy profiles for catalyst activation for two alternative dinuclear and mononuclear catalytic cycles, and for catalyst deactivation are in good agreement with the experimental results. Comparison of experimentally observed E/Z ratios at incomplete conversion with calculated kinetic selectivities revealed that a substantial amount of product must form via the dinuclear pathway, in which the isomerization is promoted cooperatively by two palladium centers. The dissociation barrier towards mononuclear Pd species is relatively high, and once the catalyst enters the energetically more favorable mononuclear pathway, only a low barrier has to be overcome towards irreversible deactivation.

Optimized damping parameters for empirical dispersion corrections to symmetry-adapted perturbation theory

Symmetry-adapted perturbation theory (SAPT) has become an invaluable tool for studying the fundamental nature of non-covalent interactions by directly computing the electrostatics, exchange (steric) repulsion, induction (polarization), and London dispersion contributions to the interaction energy using quantum mechanics. Further application of SAPT is primarily limited by its computational expense, where even its most affordable variant (SAPT0) scales as the fifth power of system size [O(N⁵)] due to the dispersion terms. The algorithmic scaling of SAPT0 is reduced from O(N⁵)→O(N⁴) by replacing these terms with the empirical D3 dispersion correction of Grimme and co-workers, forming a method that may be termed SAPT0-D3. Here, we optimize the damping parameters for the -D3 terms in SAPT0-D3 using a much larger training set than has previously been considered, namely, 8299 interaction energies computed at the complete-basis-set limit of coupled cluster through perturbative triples [CCSD(T)/CBS]. Perhaps surprisingly, with only three fitted parameters, SAPT0-D3 improves on the accuracy of SAPT0, reducing mean absolute errors from 0.61 to 0.49 kcal mol^-1 over the full set of complexes. Additionally, SAPT0-D3 exhibits a nearly 2.5× speedup over conventional SAPT0 for systems with ∼300 atoms and is applied here to systems with up to 459 atoms. Finally, we have also implemented a functional group partitioning of the approach (F-SAPT0-D3) and applied it to determine important contacts in the binding of salbutamol to G-protein coupled β₁-adrenergic receptor in both active and inactive forms. SAPT0-D3 capabilities have been added to the open-source Psi4 software.

Ab Initio Study of Chiral Discrimination in the Glycidol Dimer

Chiral discrimination, the ability of a chiral molecule to exhibit different weak intermolecular interactions than its mirror image, is investigated for dimers of oxiranemethanol (glycidol). In this regard, high-level ab initio calculations were performed to study the chiral recognition effects in the homochiral and heterochiral dimers of glycidol. Fourteen dimer structures, seven homochiral and seven heterochiral, were studied: they all feature two intermolecular O-H···O hydrogen bonds. These structures have been determined with the second-order Møller-Plesset perturbation theory (MP2) using the aug-cc-pVTZ basis set and verified to pertain to actual local minima. The benchmark interaction energy values were computed using MP2 extrapolated from the aug-cc-pVQZ and aug-cc-pV5Z bases with a higher-level correction from a coupled-cluster calculation in the aug-cc-pVTZ basis. The global minimum structure is a homochiral one, with the two hydrogen bonds forming a part of a ring with eight heavy atoms. A similar heterochiral structure has a binding energy smaller by about 0.6 kcal/mol. The largest diastereomeric energy difference is about 1.0 kcal/mol. Further insight into the origins of chiral discrimination was provided by symmetry-adapted perturbation theory (SAPT) and a functional-group SAPT (F-SAPT) difference analysis to investigate the direct and indirect effects of two -H/-CH₂OH substitutions leading from an achiral ethylene oxide dimer to the chiral glycidol dimer. Last but not least, harmonic frequency shifts relative to a noninteracting glycidol molecule were calculated and analyzed for all conformations to get insight into the origins of chiral discrimination. It is found that the largest frequency shifts are related to the effect of hydrogen bonding on the O-H stretch mode, the stability of the ring involving both hydrogen bonds, and the transition between two nonequivalent minima of the glycidol molecule.

Computationally-Guided Development of a Chelated NHC-P Iridium(I) Complex for the Directed Hydrogen Isotope Exchange of Aryl Sulfones

Herein, we report the rational, computationally-guided design of an iridium(I) catalyst system capable of enabling directed hydrogen isotope exchange (HIE) with the challenging sulfone directing group. Substrate binding energy was used as a parameter to guide rational ligand design via an in silico catalyst screen, resulting in a lead series of chelated iridium(I) NHC-phosphine complexes. Subsequent preparative studies show that the optimal catalyst system displays high levels of activity in HIE, and we demonstrate the labeling of a broad scope of substituted aryl sulfones. We also show that the activity of the catalyst is maintained at low pressures of deuterium gas and apply these conditions to tritium radiolabeling, including the expedient synthesis of a tritium-labeled drug molecule.

Non-covalent interactions and their impact on the complexation thermodynamics of noble gases with methanol

Accurate ab initio calculations provide the reliable information needed to study the potential energy surfaces that control the non-covalent interactions (NCIs) responsible for the formation of weak van der Waals complexes.

Symmetry-Adapted Perturbation Theory Decomposition of the Reaction Force: Insights into Substituent Effects Involved in Hemiacetal Formation Mechanisms

The decomposition of the reaction force based on symmetry-adapted perturbation theory (SAPT) has been proposed. This approach was used to investigate the substituent effects along the reaction coordinate pathway for the hemiacetal formation mechanism between methanol and substituted aldehydes of the form CX₃CHO (X = H, F, Cl, and Br), providing a quantitative evaluation of the reaction-driving and reaction-retarding force components. Our results highlight the importance of more favorable electrostatic and induction effects in the reactions involving halogenated aldehydes that leads to lower activation energy barriers. These substituent effects are further elucidated by applying the functional-group partition of symmetry-adapted perturbation theory (F-SAPT). The results show that the reaction is largely driven by favorable direct noncovalent interactions between the CX₃ group on the aldehyde and the OH group on methanol.

Unprecedented E-stereoselectivity on the sigmatropic Hurd-Claisen rearrangement of Morita-Baylis-Hillman adducts: a joint experimental-theoretical study

Herein we report the first systematic investigation of the tandem mercury(ii) catalysed transvinylation/Hurd-Claisen rearrangement of MBH adducts derived from alkyl acrylates. This is the first report of E-selectivity for MBH adducts with alkyl side chains and is complementary to the previously reported Johnson-Claisen and Eschenmoser-Claisen rearrangements. The rearrangement products were obtained in good yields and could be readily converted to 2-alkenyl δ-valerolactones. Combined DFT and F-SAPT studies demonstrate that reaction rates are primarily governed by non-covalent interactions dictating the relative stability of the transition states. Our F-SAPT calculations revealed that the hyperconjugative effects are not so significant, but that electrostatic interactions, instead, are the driving forces for the relative E : Z stereoselectivity.

Substituent effects in the so-called cationπ interaction of benzene and its boron-nitrogen doped analogues: overlooked role of σ-skeleton

Despite massive efforts to pinpoint the substituent effects in the so-called cationπ systems, no consensus has been yet reached on how substituents exercise their effects in the interaction of the aromatic molecule with the metal ion. The π-polarization (the Hunter model) and the direct local effect (the Wheeler-Houk model) are two lines of thought applied to this problem, but the justification of both approaches is based on insufficiently proven assumptions and approximations. In order to shed more light on this issue we propose a new approach which enables us to gauge directly the energetic trends resulting from the interaction of the ring with the cation. In our method we add one more partitioning level to the interacting quantum atoms (IQA) scheme and decompose the IQA interaction energies into contributions resulting from σ and π electron densities of the aromatic ring. The new approach, which is named partitioned-IQA, abbreviated as p-IQA, has been applied to complexes of derivatives of benzene or azaborine interacting with a sodium cation. The p-IQA approach reveals that in these systems both σ and π electronic moieties are polarized. Interestingly, for the majority of cases the σ-polarization outweighs the π one, contrary to the Hunter model. However, the Wheeler-Houk model is not precise, either, since the σ-polarization shows some degree of non-locality. In addition, the substituents are found to have a negligible influence on the ring orbital-overlapping capability, i.e. the covalency. Therefore, the substituent effect in the cationπ interaction is a nonlocal classical effect, indicating that neither Hunter model nor Wheeler-Houk model is able to fully describe all the aspects of the substituent effects. The p-IQA conclusions for the considered systems have been compared with the results from the functional-group SAPT (F-SAPT) method. We believe that the presented partitioning in the IQA framework will provide a deeper insight into the substituent effects in the cationπ interactions, which is beyond the σ-π atomic charge population separation.

A quantification scheme for non-covalent interactions in the enantio-controlling transition states in asymmetric catalysis

The origin of enantioselectivity in asymmetric catalysis is attributed to the energy difference between lower and higher energy diastereomeric transition states, which are respectively responsible for the formation of major and minor enantiomers. Although the increase in the number of transition state models emphasizes the role of weak non-covalent interactions in asymmetric induction, the strength of such interactions is seldom quantified. Through this article, we propose a simple and effective method of quantifying the total non-covalent interaction in stereocontrolling transition states belonging to a group of three representative asymmetric catalytic reactions involving chiral phosphoric acids. Our method relies on rational partitioning of a given transition state into two (or three) sub-units, such that the complex network of intramolecular interactions can be ameliorated to a set of intermolecular interactions between two sub-units. The computed strength of interaction obtained using the counterpoise (CP) method on suitably partitioned transition states provides improved estimates of non-covalent interactions, which are also devoid of basis set superposition error (BSSE). It has been noted that catalysts decorated with larger aromatic arms provide cumulative non-covalent interactions (C-Hπ, N-Hπ and ππ) to the tune of 10 to 15 kcal mol-1. Fine-tuning of the magnitude and nature of these interactions can provide valuable avenues in the design of asymmetric catalysts.

Analysis of transition state stabilization by non-covalent interactions in organocatalysis: application of atomic and functional-group partitioned symmetry-adapted perturbation theory to the addition of organoboron reagents to fluoroketones

This work seeks to apply symmetry-adapted perturbation theory (SAPT) to the recent study of Hoveyda and co-workers [K. A. Lee et al., Nat. Chem. 2016, 8, 768] where an allyl addition to a ketone became enantioselective when the ketone was fluorinated. Through the application of atomic SAPT (A-SAPT) and functional-group SAPT (F-SAPT), the non-covalent interactions between specific atoms and functional groups in the transition states associated with the fluoroketone reactions can be quantified. Our A-SAPT analysis confirms that a HF contact thought to enhance stereoselectivity shows a strong preference for one of the transition states leading to the experimentally observed product enantiomer. Other key atom-atom contacts invoked to rationalize relative transition state energies are also found to behave as expected based on chemical intuition and contact distances. On the other hand, hypothesized steric clashes between substrate phenyl or ortho-methyl phenyl groups and the catalyst are not supported by F-SAPT computations, and indeed, these are actually favorable π-π interactions.

Assessment of Density Functional Methods for Geometry Optimization of Bimolecular van der Waals Complexes

We explore the suitability of three popular density functionals (B97-D3, B3LYP-D3, M05-2X) for producing accurate equilibrium geometries of van der Waals (vdW) complexes with diverse binding motifs. For these functionals, optimizations using Dunning's aug-cc-pVDZ basis set best combine accuracy and a reasonable computational expense. Each DFT/aug-cc-pVDZ combination produces optimized equilibrium geometries for 21 small vdW complexes of organic molecules (up to four non-hydrogen atoms total) that agree with high-level CCSD(T)/CBS reference geometries to within ±0.1 Å for the averages of the center-of-mass displacement and the mean least root-mean-squared displacement. The DFT/aug-cc-pVDZ combinations are also able to reproduce the optimal center-of-mass displacements interpolated from CCSD(T)/CBS radial potential energy surfaces in both NBC7x and HBC6 test sets to within ±0.1 Å. We therefore conclude that each of these denisty functional methods, together with the aug-cc-pVDZ basis set, is suitable for producing equilibrium geometries of generic nonbonded complexes.

C-H···O Hydrogen Bonding. The Prototypical Methane-Formaldehyde System: A Critical Assessment

Distinguishing the functionality of C-H···O hydrogen bonds (HBs) remains challenging, because their properties are difficult to quantify reliably. Herein, we present a study of the model methane-formaldehyde complex (MFC). Six stationary points on the MFC potential energy surface (PES) were obtained at the CCSD(T)/ANO2 level. The CCSDT(Q)/CBS interaction energies of the conformers range from only -1.12 kcal mol^-1 to -0.33 kcal mol^-1, denoting a very flat PES. Notably, only the lowest energy stationary point (MFC1) corresponds to a genuine minimum, whereas all other stationary points-including the previously studied ideal case of a_e(C-H···O) = 180°-exhibit some degree of freedom that leads to MFC1. Despite the flat PES, we clearly see that the HB properties of MFC1 align with those of the prototypical water dimer O-H···O HB. Each HB property generally becomes less prominent in the higher-energy conformers. Only the MFC1 conformer prominently exhibits (1) elongated C-H donor bonds, (2) attractive C-H···O═C interactions, (3) n(O) → σ*(C-H) hyperconjugation, (4) critical points in the electron density from Bader's method and from the noncovalent interactions method, (5) positively charged donor hydrogen, and (6) downfield NMR chemical shifts and nonzero ²J(C_M-H_M···O_F) coupling constants. Based on this research, some issues merit further study. The flat PES hinders reliable determinations of the HB-induced shifts of the C-H stretches; a similarly difficult challenge is observed for the experiment. The role of charge transfer in HBs remains an intriguing open question, although our BLW and NBO computations suggest that it is relevant to the C-H···O HB geometries. These issues notwithstanding, the prominence of the HB properties in MFC1 serves as clear evidence that the MFC is predominantly bound by a C-H···O HB.

Psi4 1.1: An Open-Source Electronic Structure Program Emphasizing Automation, Advanced Libraries, and Interoperability

Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements.