Steric and Solvation Effects in Ionic SN2 Reactions

Spontaneous Oxidation in Aqueous Microdroplets: Water Radical Cation as Primary Oxidizing Agent

Exploration of the unique chemical properties of interfaces can unlock new understanding. A striking example is the finding of accelerated reactions, particularly spontaneous oxidation reactions, that occur without assistance of catalysts or external oxidants at the air interface of both aqueous and organic solutions (provided they contain some water). This finding opened a new area of interfacial chemistry but also caused heated debate regarding the primary chemical species responsible for the observed oxidation. An overview of the literature covering oxidation in microdroplets with air interfaces is provided, together with a critical examination of previous findings and hypotheses. The water radical cation/radical anion pair, formed spontaneously and responsible for the electric field at or near the droplet/air interface, is suggested to constitute the primary redox species. Mechanisms of accelerated microdroplet reactions are critically discussed and it is shown that hydroxyl radical/hydrogen peroxide formation in microdroplets does not require that these species be the primary oxidant. Instead, we suggest that hydroxyl radical and hydrogen peroxide are the products of water radical cation decay in water. The importance of microdroplet chemistry in the prebiotic environment is sketched briefly and the role of partial solvation in reaction acceleration is noted.

Solvent-induced dual nucleophiles and the α-effect in the SN2 versus E2 competition

We have quantum chemically investigated how microsolvation affects the various E2 and SN2 pathways, their mutual competition, and the α-effect of the model reaction system HOO-(H2O)n + CH3CH2Cl, at the CCSD(T) level. Interestingly, we identify the dual nature of the α-nucleophile HOO- which, upon solvation, is in equilibrium with HO-. This solvent-induced dual appearance gives rise to a rich network of competing reaction channels. Among both nucleophiles, SN2 is always favored over E2, and this preference increases upon increasing microsolvation. Furthermore, we found a pronounced α-effect, not only for SN2 substitution but also for E2 elimination, i.e., HOO- is more reactive than HO- in both cases. Our activation strain and quantitative molecular orbital analyses reveal the physical mechanisms behind the various computed trends. In particular, we demonstrate that two recently proposed criteria, required for solvent-free nucleophiles to display the α-effect, must also be satisfied by microsolvated HOO-(H2O)n nucleophiles.

Interfacial chemical reactivity enhancement

Interfacial enhancements of chemical reaction equilibria and rates in liquid droplets are predicted using a combined theoretical and experimental analysis strategy. Self-consistent solutions of reaction and adsorption equilibria indicate that interfacial reactivity enhancement is driven primarily by the adsorption free energy of the product (or activated complex). Reactant surface activity has a smaller indirect influence on reactivity due to compensating reactant interfacial concentration and adsorption free energy changes, as well as adsorption-induced depletion of the droplet core. Experimental air-water interfacial adsorption free energies and critical micelle concentration correlations provide quantitative surface activity estimates as a function of molecular structure, predicting an increase in interfacial reactivity with increasing product size and decreasing product polarity, aromaticity, and charge (but less so for anions than cations). Reactions with small, neutral, or charged products are predicted to have little reactivity enhancement at an air-water interface unless the product is rendered sufficiently surface active by, for example, interactions with interfacial water dangling OH groups, charge transfer, or voltage fluctuations.

Quantum Simulation of Molecules in Solution

Quantum chemical calculations on quantum computers have been focused mostly on simulating molecules in the gas phase. Molecules in liquid solution are, however, most relevant for chemistry. Continuum solvation models represent a good compromise between computational affordability and accuracy in describing solvation effects within a quantum chemical description of solute molecules. In this work, we extend the variational quantum eigensolver to simulate solvated systems using the polarizable continuum model. To account for the state dependent solute-solvent interaction we generalize the variational quantum eigensolver algorithm to treat non-linear molecular Hamiltonians. We show that including solvation effects does not impact the algorithmic efficiency. Numerical results of noiseless simulations for molecular systems with up to 12 spin-orbitals (qubits) are presented. Furthermore, calculations performed on a simulated noisy quantum hardware (IBM Q, Mumbai) yield computed solvation free energies in fair agreement with the classical calculations.

Predicting Solvent Effects on SN2 Reaction Rates: Comparison of QM/MM, Implicit, and MM Explicit Solvent Models

Solvents are one of the key variables in the optimization of a synthesis yield or properties of a synthesis product. In this paper, contemporary solvent models are applied to predict the rates of S_N2 reactions in a range of aqueous and non-aqueous solvents. High-level CCSD(T)/CBS//M06-2X/6-31+G(d) gas phase energies were combined with solvation free energies from SMD, SM12, and ADF-COSMO-RS continuum solvent models, as well as molecular mechanics (MM) explicit solvent models with different atomic charge schemes to predict the rate constants of three S_N2 reactions in eight protic and aprotic solvents. It is revealed that the prediction of rate constants in organic solvents is not necessarily less challenging than in water and popular solvent models struggle to predict their rate constants to within 3 log units of experimental values. Among the continuum solvent models, the ADF-COSMO-RS model performed the best in predicting absolute rate contants while the SM12 model was best at predicting relative rate constants with an average accuracy of about 1.5 and 0.8 log units, respectively. The use of computationally more demanding MM explicit solvent models did not translate to improvements in absolute rate constants but was quite effective at predicting relative rate constants due to systematic error cancellation. Free energy barriers obtained from umbrella sampling with explicit solvent QM/MM simulations led to excellent agreement with experimental values, provided that a validated level of theory is used to treat the QM region.

Spectroscopic Study of the Br- +CH3 I→I- +CH3 Br SN 2 Reaction

Mass spectrometry and anion photoelectron spectroscopy have been used to study the gas-phase <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>S</mml:mi> <mml:mi>N</mml:mi></mml:msub> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>${{{\rm S}}_{{\rm N}}2}$</mml:annotation> </mml:semantics> </mml:math> reaction involving <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:msup><mml:mrow><mml:mi>B</mml:mi> <mml:mi>r</mml:mi></mml:mrow> <mml:mo>-</mml:mo></mml:msup> <mml:annotation>${{{\rm B}{\rm r}}^{-}}$</mml:annotation> </mml:semantics> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow> <mml:msub><mml:mrow><mml:mi>C</mml:mi> <mml:mi>H</mml:mi></mml:mrow> <mml:mn>3</mml:mn></mml:msub> <mml:mi>I</mml:mi></mml:mrow> <mml:annotation>${{{\rm C}{\rm H}}_{3}{\rm I}}$</mml:annotation> </mml:semantics> </mml:math> . The anion photoelectron spectra associated with the reaction intermediates of this <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>S</mml:mi> <mml:mi>N</mml:mi></mml:msub> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>${{{\rm S}}_{{\rm N}}2}$</mml:annotation> </mml:semantics> </mml:math> reaction are presented. High-level CCSD(T) calculations have been utilised to investigate the reaction intermediates that may form as a result of the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>S</mml:mi> <mml:mi>N</mml:mi></mml:msub> <mml:mn>2</mml:mn></mml:mrow> <mml:annotation>${{{\rm S}}_{{\rm N}}2}$</mml:annotation> </mml:semantics> </mml:math> reaction along various different reaction pathways, including back-side attack and front-side attack. In addition, simulated vertical detachment energies of each reaction intermediate have been calculated to rationalise the photoelectron spectra.

Unexpected steric hindrance failure in the gas phase F- + (CH3)3CI SN2 reaction

Base-induced elimination (E2) and bimolecular nucleophilic substitution (S_N2) reactions are of significant importance in physical organic chemistry. The textbook example of the retardation of S_N2 reactivity by bulky alkyl substitution is widely accepted based on the static analysis of molecular structure and steric environment. However, the direct dynamical evidence of the steric hindrance of S_N2 from experiment or theory remains rare. Here, we report an unprecedented full-dimensional (39-dimensional) machine learning-based potential energy surface for the 15-atom F⁻ + (CH₃)₃CI reaction, facilitating the reliable and efficient reaction dynamics simulations that can reproduce well the experimental outcomes and examine associated atomic-molecular level mechanisms. Moreover, we found surprisingly high “intrinsic” reactivity of S_N2 when the E2 pathway is completely blocked, indicating the reaction that intends to proceed via E2 transits to S_N2 instead, due to a shared pre-reaction minimum. This finding indicates that the competing factor of E2 but not the steric hindrance determines the small reactivity of S_N2 for the F⁻ + (CH₃)₃CI reaction. Our study provides new insight into the dynamical origin that determines the intrinsic reactivity in gas-phase organic chemistry.

Base-induced elimination (E2) and bimolecular nucleophilic substitution (SN2) are of significant importance in physical organic chemistry. Here, the authors show that the competing factor of E2 as opposed to steric hindrance determines the low reactivity of SN2 in the F⁻ + (CH₃)₃CI reaction.

Evaluation of Deep Learning Architectures for Aqueous Solubility Prediction

Determining the aqueous solubility of molecules is a vital step in many pharmaceutical, environmental, and energy storage applications. Despite efforts made over decades, there are still challenges associated with developing a solubility prediction model with satisfactory accuracy for many of these applications. The goals of this study are to assess current deep learning methods for solubility prediction, develop a general model capable of predicting the solubility of a broad range of organic molecules, and to understand the impact of data properties, molecular representation, and modeling architecture on predictive performance. Using the largest currently available solubility data set, we implement deep learning-based models to predict solubility from the molecular structure and explore several different molecular representations including molecular descriptors, simplified molecular-input line-entry system strings, molecular graphs, and three-dimensional atomic coordinates using four different neural network architectures-fully connected neural networks, recurrent neural networks, graph neural networks (GNNs), and SchNet. We find that models using molecular descriptors achieve the best performance, with GNN models also achieving good performance. We perform extensive error analysis to understand the molecular properties that influence model performance, perform feature analysis to understand which information about the molecular structure is most valuable for prediction, and perform a transfer learning and data size study to understand the impact of data availability on model performance.

Proton transfer-induced competing product channels of microsolvated Y-(H2O)n + CH3I (Y = F, Cl, Br, I) reactions

The potential energy profiles of three proton transfer-involved product channels for the reactions of Y^-(H₂O)_1,2 + CH₃I (Y = F, Cl, Br, I) were characterized using the B97-1/ECP/d method. These three channels include the (1) PT_CH3 product channel that transfers a proton from methyl to nucleophile, (2) HO^--induced nucleophilic substitution (HO^--S_N2) product channel, and (3) oxide ion substitution (OIS) product channel that gives CH₃O^- and HY products. The reaction enthalpies and barrier heights follow the order OIS > PT_CH3 > HO^--S_N2 > Y^--S_N2, and thus HO^--S_N2 can compete with the most favored Y^--S_N2 product channel under singly-/doubly-hydrated conditions, while the PT_CH3 channel only occurs under high collision energy and the OIS channel is the least probable. All product channels share the same pre-reaction complex, Y^-(H₂O)_n-CH₃I, in the entrance of the potential energy profile, signifying the importance of the pre-reaction complex. For HO^-/Y^--S_N2 channels, we considered front-side attack, back-side attack, and halogen-bonded complex mechanisms. Incremental hydration increases the barriers of both HO^-/Y^--S_N2 channels as well as their barrier difference, implying that the HO^--S_N2 channel becomes less important when further hydrated. Varying the nucleophile Y^- from F^- to I^- also increases the barrier heights and barrier difference, which correlates with the proton affinity of the nucleophiles. Energy decomposition analyses show that both the orbital interaction energy and structural deformation energy of the transition states determine the S_N2 barrier change trend with incremental hydration and varying Y^-. In brief, this work computes the comprehensive potential energy surfaces of the HO^--S_N2 and PT_CH3 channels and shows how proton transfer affects the microsolvated Y^-(H₂O)_1,2 + CH₃I reaction by competing with the traditional Y^--S_N2 channel.

Structural basis for UFM1 transfer from UBA5 to UFC1

Ufmylation is a post-translational modification essential for regulating key cellular processes. A three-enzyme cascade involving E1, E2 and E3 is required for UFM1 attachment to target proteins. How UBA5 (E1) and UFC1 (E2) cooperatively activate and transfer UFM1 is still unclear. Here, we present the crystal structure of UFC1 bound to the C-terminus of UBA5, revealing how UBA5 interacts with UFC1 via a short linear sequence, not observed in other E1-E2 complexes. We find that UBA5 has a region outside the adenylation domain that is dispensable for UFC1 binding but critical for UFM1 transfer. This region moves next to UFC1’s active site Cys and compensates for a missing loop in UFC1, which exists in other E2s and is needed for the transfer. Overall, our findings advance the understanding of UFM1’s conjugation machinery and may serve as a basis for the development of ufmylation inhibitors.

Ufmylation is a well-established ubiquitin-like protein modification, but its mechanism is largely unclear. Here, the authors present a crystal structure of the ufmylation-specific E1-E2 complex, revealing differences to the ubiquitination machinery and mechanistic details of the ufmylation process.

Energetic Descriptors of Steric Hindrance in Real Space: An Improved IQA Picture*

Steric hindrance (SH) plays a central role in the modern chemical narrative, lying at the core of chemical intuition. As it however happens with many successful chemical concepts, SH lacks an underlying physically sound root, and multiple mutually inconsistent approximations have been devised to relate this fuzzy concept to computationally derivable descriptors. We here argue that being SH related to spatial as well as energetic features of interacting systems, SH can be properly handled if we chose a real space energetic stance like the Interacting Quantum Atoms (IQA) approach. Drawing on previous work by Popelier and coworkers (ChemistryOpen 8, 560, 2019) we build an energetic estimator of SH, referred to as E_ST . We show that the rise in the self-energy of a fragment that accompanies steric congestion is a faithful proxy for the chemist's SH concept if we remove the effect of charge transfer. This can be done rigorously, and the E_ST here defined provides correct sterics even for hydrogen atoms, where the plain use of deformation energies leads to non-chemical results. The applicability of E_ST is validated in several chemical scenarios, going from atomic compressions to archetypal S_N2 reactions. E_ST is shown to be a robust steric hindrance descriptor.

Transition-State Expansion: A Quantitative Model for Counterion Effects in Ionic Reactions

Ionic reactions are the most common reactions used in chemical synthesis. In relatively low dielectric constant solvents (e.g., dichloromethane, toluene), ions usually exist as ion pairs. Despite the importance of counterions, a quantitative description of how the paired 'counterion' affects the reaction kinetic is still elusive. We introduce a general and quantitative model, namely transition-state expansion (TSE), that describes how the size of a counterion affects the transition-state structure and the kinetics of an ionic reaction. This model could rationalize the counterion effects in nucleophilic substitutions and gold-catalyzed enyne cycloisomerizations.

•

A quantitative model for counterion effects was introduced

•

'Transition State Expansion' was used to describes how the size of a counterion affects the kinetics

Solvent Effects on the Symmetric and Asymmetric SN2 Reactions in the Acetonitrile Solution: A Reaction Density Functional Theory Study

Bimolecular nucleophilic substitution (S_N2) reactions are of great importance in chemistry and biochemistry due to their capability of constructing functional groups. In this work, we investigate the solvent effect on the free energy profiles of symmetric and asymmetric S_N2 reactions in the acetonitrile solution using the proposed reaction density functional theory (RxDFT) method. This multiscale method utilizes quantum density functional theory for calculating intrinsic reaction free energy coupled with classical density functional theory for addressing solvation contribution. We find that the presence of acetonitrile brings both the polarization effect and solvation effect on the reaction pathways. For the eight selected symmetric S_N2 reactions, the predicated reaction pathways agree well with the results from the direct and thermodynamic cycle (TC) methods with the SMD-M062X solvation model. In addition, the polarization effect reduces the free energy barriers by about 6 kcal/mol, while the solvation effect increases the barriers by about 18 kcal/mol. For the four selected asymmetric S_N2 reactions, the predicted reaction pathways agree well with the results from the Monte Carlo simulations and experiments. The polarization effect and the solvation effect mutually reduce the free energy barriers, and the solvation effect plays a dominant role.

Ambident Nucleophilic Substitution: Understanding Non-HSAB Behavior through Activation Strain and Conceptual DFT Analyses

The ability to understand and predict ambident reactivity is key to the rational design of organic syntheses. An approach to understand trends in ambident reactivity is the hard and soft acids and bases (HSAB) principle. The recent controversy over the general validity of this principle prompted us to investigate the competing gas-phase SN 2 reaction channels of archetypal ambident nucleophiles CN- , OCN- , and SCN- with CH3 Cl (SN 2@C) and SiH3 Cl (SN 2@Si), using DFT calculations. Our combined analyses highlight the inability of the HSAB principle to correctly predict the reactivity trends of these simple, model reactions. Instead, we have successfully traced reactivity trends to the canonical orbital-interaction mechanism and the resulting nucleophile-substrate interaction energy. The HOMO-LUMO orbital interactions set the trend in both SN 2@C and SN 2@Si reactions. We provide simple rules for predicting the ambident reactivity of nucleophiles based on our Kohn-Sham molecular orbital analysis.

QM/MM Calculations for the Cl- + CH3Cl SN2 Reaction in Water Using CM5 Charges and Density Functional Theory

The prototypical S_N2 reaction of chloride ion with methyl chloride has been reinvestigated in aqueous solution using QM/MM methodology featuring MO6-2X/6-31+G(d) calculations with the TIP4P water model, and partial charges were computed with the CM5 method. Though the DFT method yields excellent gas-phase energetics for the reaction, the QM/MM approach is found to yield overestimation of the activation barrier by ca. 12 kcal/mol. The discrepancy is traced to underestimate of the magnitude of the partial charges on the chlorine atoms in the transition structure. When CM1 or CM3 charges based on semiempirical wave functions are used instead, the agreement with experiment is much improved. The findings emphasize the sensitivity of the results of QM/MM calculations to the choice of QM method, the MM force field, and implementation of the QM/MM interface.

A reaction density functional theory study of the solvent effect in prototype SN2 reactions in aqueous solution

Reaction density functional theory (RxDFT), combining quantum DFT with classical DFT, has been employed to investigate the solvent effect and free energy profiles of S_N2 reactions in aqueous solution.

Nucleophilic Substitution (SN 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent

The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.

Ion-Pair SN 2 Reaction of OH- and CH3 Cl: Activation Strain Analyses of Counterion and Solvent Effects

We have theoretically studied the non-identity S_N 2 reactions of M_n OH^(n-1) +CH₃ Cl (M⁺ =Li⁺ , Na⁺ , K⁺ , and MgCl⁺ ; n=0, 1) in the gas phase and in THF solution at the OLYP/6-31++G(d,p) level using polarizable continuum model (PCM) implicit solvation. We want to explore and understand the effect of the metal counterion M⁺ and solvation on the reaction profile and the stereoselectivity of these processes. To this end, we have explored the potential energy surfaces of the backside (S_N 2-b) and frontside (S_N 2-f) pathways. To explain the computed trends, we have carried out analyses with an extended activation strain model (ASM) of chemical reactivity that includes the treatment of solvation effects.

Microsolvation effects on the reactivity of oxy-nucleophiles: the case of gas-phase SN2 reactions of YO-(CH3OH) n=1,2 towards CH3Cl

The modified G4(MP2) method was applied to explore microsolvation effects on the reactivity of four solvated normal oxy-nucleophiles YO^-(CH₃OH) _n=1,2 (Y = CH₃, C₂H₅, FC₂H₄, ClC₂H₄), and five α-oxy-nucleophiles YO^-(CH₃OH) _n=1,2 (Y = HO, CH₃O, F, Cl, Br), in gas-phase S_N2 reactions towards the substrate CH₃Cl. Based on a Brønsted-type plot, our calculations reveal that the overall activation barriers of five microsolvated α-oxy-nucleophiles are obviously smaller than the prediction from the correlation line constructed by four normal microsolvated ones to different degrees, and clearly demonstrate the existence of an α-effect in the presence of one or two methanol molecule(s). Moreover, it was found that the α-effect of the mono-methanol microsolvated α-nucleophile is stronger than that of the monohydrated α-nucleophile. However, the α-effect of YO^-(CH₃OH)₂ becomes weaker for Y = HO and CH₃O, whereas it becomes stronger for Y = F, Cl, Br than that of YO^-(H₂O)₂, which can be explained by analyses of the activation strain model in the two cases. It was also found that the rationale about the low ionization energy of α-nucleophile inducing the α-effect was not widely significant. Graphical abstract Variation of alpha-effect in the gas-phase S_N2 reaction with the microsolvation.

Biological Activity of Ionic Liquids and Their Application in Pharmaceutics and Medicine

Ionic liquids are remarkable chemical compounds, which find applications in many areas of modern science. Because of their highly tunable nature and exceptional properties, ionic liquids have become essential players in the fields of synthesis and catalysis, extraction, electrochemistry, analytics, biotechnology, etc. Apart from physical and chemical features of ionic liquids, their high biological activity has been attracting significant attention from biochemists, ecologists, and medical scientists. This Review is dedicated to biological activities of ionic liquids, with a special emphasis on their potential employment in pharmaceutics and medicine. The accumulated data on the biological activity of ionic liquids, including their antimicrobial and cytotoxic properties, are discussed in view of possible applications in drug synthesis and drug delivery systems. Dedicated attention is given to a novel active pharmaceutical ingredient-ionic liquid (API-IL) concept, which suggests using traditional drugs in the form of ionic liquid species. The main aim of this Review is to attract a broad audience of chemical, biological, and medical scientists to study advantages of ionic liquid pharmaceutics. Overall, the discussed data highlight the importance of the research direction defined as "Ioliomics", studies of ions in liquids in modern chemistry, biology, and medicine.

Ion-Pair SN 2 Substitution: Activation Strain Analyses of Counter-Ion and Solvent Effects

The ion-pair SN 2 reactions of model systems MnF(n-1) +CH3Cl(M(+) =Li(+), Na(+), K(+), and MgCl(+); n=0, 1) have been quantum chemically explored by using DFT at the OLYP/6-31++G(d,p) level. The purpose of this study is threefold: 1) to elucidate how the counterion M(+) modifies ion-pair SN 2 reactivity relative to the parent reaction F(-) +CH3Cl; 2) to determine how this influences stereochemical competition between the backside and frontside attacks; and 3) to examine the effect of solvation on these ion-pair SN2 pathways. Trends in reactivity are analyzed and explained by using the activation strain model (ASM) of chemical reactivity. The ASM has been extended to treat reactivity in solution. These findings contribute to a more rational design of tailor-made substitution reactions.

Calculating Free Energy Changes in Continuum Solvation Models

We recently showed for a large data set of pKas and reduction potentials that free energies calculated directly within the SMD continuum model compares very well with corresponding thermodynamic cycle calculations in both aqueous and organic solvents [ Phys. Chem. Chem. Phys. 2015 , 17 , 2859 ]. In this paper, we significantly expand the scope of our study to examine the suitability of this approach for calculating general solution phase kinetics and thermodynamics, in conjunction with several commonly used solvation models (SMD-M062X, SMD-HF, CPCM-UAKS, and CPCM-UAHF) for a broad range of systems. This includes cluster-continuum schemes for pKa calculations as well as various neutral, radical, and ionic reactions such as enolization, cycloaddition, hydrogen and chlorine atom transfer, and SN2 and E2 reactions. On the basis of this benchmarking study, we conclude that the accuracies of both approaches are generally very similar-the mean errors for Gibbs free energy changes of neutral and ionic reactions are approximately 5 and 25 kJ mol(-1), respectively. In systems where there are significant structural changes due to solvation, as is the case for certain ionic transition states and amino acids, the direct approach generally afford free energy changes that are in better agreement with experiment.

New insights into hydrosilylation of unsaturated carbon-heteroatom (C═O, C═N) bonds by rhenium(V)-dioxo complexes

The hydrosilylation of unsaturated carbon-heteroatom (C═O, C═N) bonds catalyzed by high-valent rhenium(V)-dioxo complex ReO2I(PPh3)2 (1) were studied computationally to determine the underlying mechanism. Our calculations revealed that the ionic outer-sphere pathway in which the organic substrate attacks the Si center in an η(1)-silane rhenium adduct to prompt the heterolytic cleavage of the Si-H bond is the most energetically favorable process for rhenium(V)-dioxo complex 1 catalyzed hydrosilylation of imines. The activation energy of the turnover-limiting step was calculated to be 22.8 kcal/mol with phenylmethanimine. This value is energetically more favorable than the [2 + 2] addition pathway by as much as 10.0 kcal/mol. Moreover, the ionic outer-sphere pathway competes with the [2 + 2] addition mechanism for rhenium(V)-dioxo complex 1 catalyzing the hydrosilylation of carbonyl compounds. Furthermore, the electron-donating group on the organic substrates would induce a better activity favoring the ionic outer-sphere mechanistic pathway. These findings highlight the unique features of high-valent transition-metal complexes as Lewis acids in activating the Si-H bond and catalyzing the reduction reactions.

Perspective: Bimolecular chemical reaction dynamics in liquids

Bimolecular reactions in the gas phase exhibit rich and varied dynamical behaviour, but whether a profound knowledge of the mechanisms of isolated reactive collisions can usefully inform our understanding of reactions in liquid solutions remains an open question. The fluctuating environment in a liquid may significantly alter the motions of the reacting particles and the flow of energy into the reaction products after a transition state has been crossed. Recent experimental and computational studies of exothermic reactions of CN radicals with organic molecules indicate that many features of the gas-phase dynamics are retained in solution. However, observed differences may also provide information on the ways in which a solvent modifies fundamental chemical mechanisms. This perspective examines progress in the use of time-resolved infra-red spectroscopy to study reaction dynamics in liquids, discusses how existing theories can guide the interpretation of experimental data, and suggests future challenges for this field of research.

Quantum and Molecular Mechanical (QM/MM) Monte Carlo Techniques for Modeling Condensed-Phase Reactions

A recent review (Acc. Chem. Res. 2010, 43:142-151) examined our use and development of a combined quantum and molecular mechanical (QM/MM) technique for modelling organic and enzymatic reactions. Advances included the PDDG/PM3 semiempirical QM (SQM) method, computation of multi-dimensional potentials of mean force (PMF), incorporation of on-the-fly QM in Monte Carlo simulations, and a polynomial quadrature method for rapidly treating proton-transfer reactions. The current article serves as a follow up on our progress. Highlights include new reactions, alternative SQM methods, a polarizable OPLS force field, and novel solvent environments, e.g., "on water" and room temperature ionic liquids. The methodology is strikingly accurate across a wide range of condensed-phase and antibody-catalyzed reactions including substitution, decarboxylation, elimination, isomerization, and pericyclic classes. Comparisons are made to systems treated with continuum-based solvents and ab initio or density functional theory (DFT) methods. Overall, the QM/MM methodology provides detailed characterization of reaction paths, proper configurational sampling, several advantages over implicit solvent models, and a reasonable computational cost.

How alkyl halide structure affects E2 and SN2 reaction barriers: E2 reactions are as sensitive as SN2 reactions

High-level electronic structure calculations, including a continuum treatment of solvent, are employed to elucidate and quantify the effects of alkyl halide structure on the barriers of SN2 and E2 reactions. In cases where such comparisons are available, the results of these calculations show close agreement with solution experimental data. Structural factors investigated include α- and β-methylation, adjacency to unsaturated functionality (allyl, benzyl, propargyl, α to carbonyl), ring size, and α-halogenation and cyanation. While the influence of these factors on SN2 reactivity is mostly well-known, the present study attempts to provide a broad comparison of both SN2 and E2 reactivity across many cases using a single methodology, so as to quantify relative reactivity trends. Despite the fact that most organic chemistry textbooks say far more about how structure affects SN2 reactions than about how it affects E2 reactions, the latter are just as sensitive to structural variation as are the former. This sensitivity of E2 reactions to structure is often underappreciated.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

SN2 regioselectivity in the esterification of 5- and 7-membered azacycloalkane quaternary salts: a DFT study to reveal the transition state ring conformation prevailing over the ground state ring strain

S_N2 regioselectivity in 5- and 7-membered azacycloalkanes quaternary salts is directed by the transition state ring conformation.

Investigating the α-effect in gas-phase S(N)2 reactions of microsolvated anions

The α-effect-enhanced reactivity of nucleophiles with a lone-pair adjacent to the attacking center-was recently demonstrated for gas-phase S(N)2 reactions of HOO(-), supporting an intrinsic component of the α-effect. In the present work we explore the gas-phase reactivity of microsolvated nucleophiles in order to investigate in detail how the α-effect is influenced by solvent. We compare the gas-phase reactivity of the microsolvated α-nucleophile HOO(-)(H2O) to that of microsolvated normal alkoxy nucleophiles, RO(-)(H2O), in reaction with CH3Cl using a flowing afterglow-selected ion flow tube instrument. The results reveal enhanced reactivity of HOO(-)(H2O) and clearly demonstrate the presence of an α-effect for the microsolvated α-nucleophile. The association of the nucleophile with a single water molecule results in a larger Brønsted βnuc value than is the case for the unsolvated nucleophiles. Accordingly, the reactions of the microsolvated nucleophiles proceed through later transition states in which bond formation has progressed further. Calculations show a significant difference in solvent interaction for HOO(-) relative to the normal nucleophiles at the transition states, indicating that differential solvation may well contribute to the α-effect. The reactions of the microsolvated anions with CH3Cl can lead to formation of either the bare Cl(-) anion or the Cl(-)(H2O) cluster. The product distributions show preferential formation of the Cl(-) anion even though the formation of Cl(-)(H2O) would be favored thermodynamically. Although the structure of the HOO(-)(H2O) cluster resembles HO(-)(HOOH), we demonstrate that HOO(-) is the active nucleophile when the cluster reacts.

Impact of alkyl substituents on the gas-phase competition between substitution and elimination

The S(N)2 and E2 reactions of a series of alkyl bromides with varying substitution patterns at the α- and β-carbons have been studied in the gas phase using naphthoate and phenoxide-based nucleophiles. The experimental work is supported by calculations at the MP2/6-31+G(d,p)//MP2/6-31+G(d) level. The results parallel reactivity patterns observed in the condensed phase, but offer new insights into steric factors in S(N)2 processes. In the gas phase, polarizability is more important, and the highest S(N)2 reactivity is observed when the β-carbon is 2°. In addition, the data confirm that alkyl substituents at the β-carbon have a greater accelerating effect on E2 reactions than those at the α-carbon. Finally, computed data based on lowest enthalpy pathways provide poor descriptions of the reactions of the larger alkyl bromides and are skewed toward crowded systems that offer stabilizing, nonbonded interactions at the expense of conformational freedom.

Benchmark study on the smallest bimolecular nucleophilic substitution reaction: H⁻+CH₄-->CH₄+H⁻

We report here a benchmark study on the bimolecular nucleophilic substitution (S(N)2) reaction between hydride and methane, for which we have obtained reference energies at the coupled cluster toward full configuration-interaction limit (CC-cf/CBS). Several wavefunction (HF, MP2, coupled cluster) and density functional methods are compared for their reliability regarding these reference data.

THE MECHANISM OF ACID-CATALYZED DECARBOXYLATION OF PYRROLE-2-CARBOXYLIC ACID: INSIGHTS FROM CLUSTER-CONTINUUM MODEL CALCULATIONS

The decarboxylation of pyrrole-2-carboxylic acid comprises the addition of water to the carboxyl group and the C–C bond cleavage leading to the protonated carbonic acid. Herein possible concerted and stepwise mechanisms for the C-protonated and O-protonated pathways were extensively investigated by using the cluster-continuum model. The calculated results indicate that the initial hydration or the nucleophilic attack of water at the carbonyl group of both C- and O-protonated derivatives is the rate-determining step for the overall reaction, and the O-protonated pathway will dominate the whole reaction. The predicted activation Gibbs energies for the overall reaction initialized by the O-protonated species fall in the range of 83.3 ∼ 123.0 kJ/mol, showing good agreement with experimental values of 91.6 ∼ 101.3 kJ/mol. On the basis of extensive calculations, the remarkable dependence of the predicted mechanisms and thermodynamic values on the number of explicit water molecules in the cluster-continuum model was discussed.