Empirical valence bond model of an SN2 reaction in polar and nonpolar solvents

Enhanced ClNO2 Formation at the Interface of Sea-Salt Aerosol

The reactive uptake of N2O5 on sea-spray aerosol plays a key role in regulating the NOx concentration in the troposphere. Despite numerous field and laboratory studies, a microscopic understanding of its heterogeneous reactivity remains unclear. Here, we use molecular simulation and theory to elucidate the chlorination of N2O5 to form ClNO2, the primary reactive channel within sea-spray aerosol. We find that the formation of ClNO2 is markedly enhanced at the air-water interface due to the stabilization of the charge-delocalized transition state, as evident from the formulation of bimolecular rate theory in heterogeneous environments. We explore the consequences of the enhanced interfacial reactivity in the uptake of N2O5 using numerical solutions of molecular reaction-diffusion equations as well as their analytical approximations. Our results suggest that the current interpretation of aerosol branching ratios needs to be revisited.

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.

Alkyl hydrogen atom abstraction reactions of the CN radical with ethanol

We present a study of the abstraction of alkyl hydrogen atoms from the β and α positions of ethanol by the CN radical in solution using the Empirical Valence Bond (EVB) method. We have built separate 2 × 2 EVB models for the H_β and H_α reactions, where the atom transfer is parameterized using ab initio calculations. The intra- and intermolecular potentials of the reactant and product molecules were modelled with the General AMBER Force Field, with some modifications. We have carried out the dynamics in water and chloroform, which are solvents of contrasting polarity. We have computed the potential of mean force for both abstractions in each of the solvents. They are found to have a small and early barrier along the reaction coordinate with a large energy release. Analyzing the solvent structure around the reaction system, we have found two solvents to have little effect on either reaction. Simulating the dynamics from the transition state, we also study the fate of the energies in the HCN vibrational modes. The HCN molecule is born vibrationally hot in the CH stretch in both reactions and additionally in the HCN bends for the H_α abstraction reaction. In the early stage of the dynamics, we find that the CN stretch mode gains energy at the expense of the energy in CH stretch mode.

Reaction dynamics at liquid interfaces

The liquid interface is a narrow, highly anisotropic region, characterized by rapidly varying density, polarity, and molecular structure. I review several aspects of interfacial solvation and show how these affect reactivity at liquid/liquid interfaces. I specifically consider ion transfer, electron transfer, and SN2 reactions, showing that solvent effects on these reactions can be understood by examining the unique structure and dynamics of the liquid interface region.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

On unjustifiably misrepresenting the EVB approach while simultaneously adopting it

In recent years, the EVB has become a widely used tool in the QM/MM modeling of reactions in condensed phases and in biological systems, with ever increasing popularity. However, despite the fact that its power and validity have been repeatedly established since 1980, a recent work (Valero, R.; et al. J. Chem. Theory Comput. 2009, 5, 1) has strongly criticized this approach, while not discussing the fact that one of the authors is effectively using it himself for both gas-phase and solution studies. Here, we have responded to the most serious unjustified assertions of that paper, covering both the more problematic aspects of that work and the more complex scientific aspects. Additionally, we have demonstrated that the poor EVB results shown in Valero et al. which where presented as verification of the unreliability of the EVB model were in fact obtained by the use of incorrect parameters, without comparing to the correct surface obtained by our program.