Kinetic study of thermal Z to E isomerization reactions of azobenzene and 4-dimethylamino-4'-nitroazobenzene in ionic liquids [1-R-3-methylimidazolium bis(trifluoromethylsulfonyl)imide with R = butyl, pentyl, and hexyl

Thermal Z to E isomerization reactions of azobenzene and 4-dimethylamino-4'-nitroazobenzene were examined in three ionic liquids of general formula 1-R-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (R = butyl, pentyl, and hexyl). The first-order rate constants and activation energies for the reactions of azobenzene measured in these ionic liquids were consistent with those measured in ordinary organic solvents, which indicated that the slow isomerization through the inversion mechanism with a nonpolar transition state was little influenced by the solvent properties, such as the viscosity and dielectric constant, of ionic liquids. On the other hand, the rate constants and the corresponding frequency factors of the Arrhenius plot were significantly reduced for the isomerization of 4-dimethylamino-4'-nitroazobenzene in ionic liquids compared with those for the isomerization in ordinary organic molecular solvents with similar dielectric properties. Although these ionic liquids are viscous, the apparent viscosity dependence of the rate constant could not be explained either by the Kramers-Grote-Hynes model or by the Agmon-Hopfield model for solution reactions. It is proposed that the positive and the negative charge centers of a highly polar rotational transition state are stabilized by the surrounding anions and cations, respectively, and that the ions must be rearranged so as to form highly ordered solvation shells around the charge centers of the reactant in the transition state. This requirement for the orderly solvation in the transition state results in unusually small frequency factors of 10(4)-10(7) s(-1).     

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 2. Application to transphosphorylation thio effects in solution

Density-functional and semiempirical quantum methods and continuum dielectric and explicit solvation models are applied to study the role of solvation on the stabilization of native and thio-substituted transphosphorylation reactions. Extensive comparison is made between results obtained from the different methods. For the semiempirical methods, explicit solvation was treated using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach and the implicit solvation was treated using a recently developed smooth solvation model implemented into a d-orbital semiempirical framework (MNDO/d-SCOSMO) within CHARMM. The different quantum and solvation methods were applied to the transesterification of a 3'-ribose,5'-methyl phosphodiester that serves as a nonenzymatic model for the self-cleavage reaction catalyzed by the hammerhead and hairpin ribozymes. Thio effects were studied for a double sulfur substitution at the nonbridging phosphoryl oxygen positions. The reaction profiles of both the native and double sulfur-substituted reactions from the MNDO/d-SCOSMO calculations were similar to the QM/MM results and consistent with the experimentally observed trends. These results underscore the need for a d-orbital semiempirical representation for phosphorus and sulfur for the study of experimentally observed thio effects in enzymatic and nonenzymatic phosphoryl transfer reactions. One of the major advantages of the present approach is that it can be applied to model chemical reactions at a significantly lower computational cost than either the density-functional calculations with implicit solvation or the semiempirical QM/MM simulations with explicit solvent.     

Computational Analysis of Solute–Solvent Coupling Magnitude in the Z/E Isomerization Reaction of Nitroazobenzene and Benzylideneanilines

The dynamic solvent effect often arises in solution reactions, where coupling between chemical reaction and solvent fluctuation plays a decisive role in the reaction kinetics. In this study, the Z/E isomerization reaction of nitoroazobenzene and benzylideneanilines in the ground state was computationally studied by molecular dynamics simulations. The non-equilibrium solvation effect was analyzed using two approaches: (1) metadynamics Gibbs energy surface exploration and (2) solvation Gibbs energy evaluation using a frozen solvation droplet model. The solute–solvent coupling parameter (C_coupled) was estimated by the ratio of the solvent fluctuation Gibbs energy over the corresponding isomerization activation Gibbs energy. The results were discussed in comparison with the ones estimated by means of the analytical models based on a reaction–diffusion equation with a sink term. The second approach using a frozen solvation droplet reached qualitative agreement with the analytical models, while the first metadynamics approach failed. This is because the second approach explicitly considers the non-equilibrium solvation in the droplet, which consists of a solute at the reactant geometry immersed in the pre-organized solvents fitted with the solute at the transition state geometry.     

Experimental and theoretical studies on the isomerization of allyl thiocyanate to allyl isothiocyanate

The mechanism of isomerization of allyl thiocyanate to allyl isothiocyanate has been investigated both experimentally and theoretically. The kinetic study indicates that the reaction is unimolecular and is not ionic. The entropy of activation suggests strongly that the mechanism involves a cyclic transition state. The rate of reaction was retarded to a small extent in polar solvents relative to that in nonpolar solvents. Ab initio MO calculations indicate, in agreement with the experimental results, that the reaction proceeds through a cyclic transition state, one in which the SCN moiety is almost linear. Thus, this is a [3,3] sigmatropic rearrangement. The charge separation in the transition state was substantial. The retardation of the reaction in polar solvents was attributed to the difference in solvation in the original state and in the transition state. © 1997 John Wiley & Sons, Inc.     

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 1. Implementation

The present paper describes the extension of a recently developed smooth conductor-like screening model for solvation to a d-orbital semiempirical framework (MNDO/d-SCOSMO) with analytic gradients that can be used for geometry optimizations, transition state searches, and molecular dynamics simulations. The methodology is tested on the potential energy surfaces for separating ions and the dissociative phosphoryl transfer mechanism of methyl phosphate. The convergence behavior of the smooth COSMO method with respect to discretization level is examined and the numerical stability of the energy and gradient are compared to that from conventional COSMO calculations. The present method is further tested in applications to energy minimum and transition state geometry optimizations of neutral and charged metaphosphates, phosphates, and phosphoranes that are models for stationary points in transphosphorylation reaction pathways of enzymes and ribozymes. The results indicate that the smooth COSMO method greatly enhances the stability of quantum mechanical geometry optimization and transition state search calculations that would routinely fail with conventional solvation methods. The present MNDO/d-SCOSMO method has considerable computational advantages over hybrid quantum mechanical/molecular mechanical methods with explicit solvation, and represents a potentially useful tool in the arsenal of multi-scale quantum models used to study biochemical reactions.     

Self-consistent combination of the three-dimensional RISM theory of molecular solvation with analytical gradients and the Amsterdam density functional package

The three-dimensional reference interaction site model with the closure relation by Kovalenko and Hirata (3D-RISM-KH) in combination with the density functional theory (DFT) method has been implemented in the Amsterdam density functional (ADF) software package. The analytical first derivatives of the free energy with respect to displacements of the solute nuclear coordinates have also been developed. This enables study of chemical reactions, including reaction coordinates and transition state search, with the molecular solvation described from the first principles. The method yields all of the features available by using other solvation approaches, for instance infrared spectra of solvated molecules. To evaluate the accuracy of the present method, test calculations have been carried out for a number of small molecules, including four glycine conformers, a set of small organic compounds, and carbon nanotubes of various lengths in aqueous solution. Our predictions for the solvation free energy agree well with other approaches as well as experiment. This new development makes it possible to calculate at modest computational cost the electronic properties and molecular solvation structure of a solute molecule in a given molecular liquid or mixture from the first principles.     

Aquation of the ruthenium-based anticancer drug NAMI-A: a density functional study

We carried out density functional theory (DFT) calculations to investigate the thermodynamics and the kinetics of the double aquation reaction of the anticancer drug NAMI-A. Three explicit water molecules were included in the calculations to improve the PB solvation energies. Our calculations show that the chloride substitution reactions on the considered Ru(III) octahedral complex follow a dissociative interchange mechanism, I(d), passing through a loose heptacoordinate transition state. We calculated an activation enthalpy and free energy for the first aquation step of 101.5 and 103.7 kJ mol(-1), respectively, values that are in good agreement with the available experimental results. The activation enthalpy and free energy for the second aquation step were found significantly higher, 118.7 and 125.0 kJ mol(-1), again in agreement with the experimental evidence indicating a slower rate for the second aquation.     

Nucleophilic substitution reactions at liquid/liquid interfaces: molecular dynamics simulation of a model S(N)1 dissociation reaction at the water/carbon tetrachloride interface

The ionic dissociation step of the nucleophilic substitution reaction t-BuCl --> t-Bu(+) + Cl(-) is studied at the water/carbon tetrachloride interface using molecular dynamics computer simulations. The empirical valence bond approach is used to couple two diabatic states, covalent and ionic, in the electronically adiabatic limit. The umbrella sampling technique is used to calculate the potential of mean force along the reaction coordinate (defined as the t-Bu to Cl distance) at several interface regions of varying distances from the Gibbs dividing surface. We find a significant increase of the ionic dissociation barrier height and of the reaction free energy at the interface relative to bulk water. This is shown to be due to the reduced polarity of the interface which causes a destabilization of the pure ionic state. However, deformation to the neat interface structure in the form of water protrusions into the organic phase may provide partial stabilization of the ionic species. The importance of these structural effects is examined by repeating the calculations with an artificially smooth interface. The destabilization of the ionic state at the interface also manifests itself with a rapid (picosecond time scale) recombination dynamics of the ions to form the parent molecule followed by a slow vibrational relaxation.     

Chemoselective nucleophilic fluorination induced by selective solvation of the SN2 transition state

Reaction of the fluoride ion with secondary alkyl halides leads to 90% of elimination reaction and only 10% of nucleophilic substitution in dipolar aprotic solvents. Adding water to the organic phase, the SN2 yield increases in the cost of decreased reactivity. Using ab initio calculations, we have shown that it is possible to increase the reaction rate and the selectivity toward the SN2 process through supramolecular organocatalysis. The catalytic concept is based on selective solvation of the transition state through two hydrogen bonds provided by the 1,4-benzenedimethanol. The two hydrogen bonds between the catalyst and the SN2 transition state favor this pathway while just one strong hydrogen bond between the catalyst and the fluoride ion leads to a lower stabilization of the nucleophile, resulting in a higher reaction rate. Our calculations predict that the substitution product increases to 40% yield because of the selective catalysis provided by the 1,4-benzenedimethanol.     

杂环反铂(Ⅱ)抗癌药物水解反应机理的DFT研究

用DFT-B3LYP方法和IEF-PCM溶剂化模型研究了反铂抗癌药物trans-[PtCl2(piperidine)(Am)](Am=2-picoline(1),3-picoline(2),4-picoline(3)),trans-[PtCl2(piperidine)(piperazine)](4),trans-[PtCl2(pipera-zine)2](5)and trans-[PtCl2(iminoether)2](6)的水解过程.水解反应是药物与DNA靶分子作用的关键活化步骤.全优化和表征了一水解和二水解反应经由一般的SN2路径过程所有物种的势能面稳定点.结果发现反应过程遵循已经建立的平面正方形配合物的配体取代反应理论,即取代反应通常通过一个三角双锥过渡态结构的铂配体交换反应发生.得到的过渡态结构与以前的相关工作一致,所有反应都是吸热反应;所有体系的二水解能垒都高于一水解.与顺铂相比,这些配合物都有更快的水解反应速率;并与以前类似的反铂配合物的研究做了比较.研究结果提供了这些配合物水解反应过程的详细能量变化,对理解药物与DNA靶分子的作用机理和新型反铂抗癌药物的设计有帮助.     

Hybrid QM/MM molecular dynamics simulations for an ionic SN2 reaction in the supercritical water: OH- + CH3Cl --> CH3OH + Cl-

A hybrid real space quantum mechanical/molecular mechanical (RS-QM/MM) method has been applied to an ionic S(N)2 reaction (OH- + CH3Cl --> CH3OH + Cl-) in water solution to investigate dynamic solvation effects of the supercritical water (SCW) on the reaction. It has been demonstrated that the approaching process of OH- to methyl group is prevented by water molecules in the ambient water (AW), while the reaction takes place easily in the gas phase. Almost the same solvation effect on the dynamics of OH- is observed in the SCW, though the bulk density of water is substantially reduced compared with that of the AW. It has been shown that the solvation of the SCW around the OH anion is locally identical to that of the AW due to the strong ion-dipole interactions between OH- and water molecules. At the transition state, the QM/MM simulations have revealed that the excess electron is quite flexible, and the charge volume, as well as the fractional charges on atoms, vary seriously depending on the instantaneous solvent configurations. However, it has been found that the solvation energy in the SCW can be qualitatively related to the HOMO volume of the system by Born's equation.     

Quantum mechanics investigation on initial decomposition of ammonia borane in glyme

The chemical kinetics of ammonia borane (AB) in glyme solution is studied using quantum mechanics (QM) based calculations along with experimental results available in the literature. The primary objective of this study is to propose a detailed reaction mechanism that explains the formation of species observed during AB decomposition for temperatures ranging from 323 to 368 K. The quantum mechanics investigation uses transition state theory to identify the relevant reaction pathways. Intrinsic reaction coordinate calculations use the identified transition‐state structure to link the reactants to the products. These calculations were performed using the Gaussian 09 program package, including the solvation model based on density (SMD) with acetonitrile as the solvent. Thermodynamic properties of species at equilibrium or at transition states were computed using the G4(MP2) compound method. Sensitivity analysis was performed using a species conservation model to identify reactions and species that play a critical role. This study confirms the previous experimental observation regarding the initiation of decomposition of AB in glyme. It also elucidates the role of DADB, ammonium borohydride salt ([BH₄]⁻[NH₄]⁺) and BH₂NH₂ in hydrogen release and intermediates formed during initial phase of AB decomposition. This work shows how QM calculations along with experimental results can contribute to our understanding of the complex chemical kinetics involved during AB dehydrogenation.     

Uncatalyzed transfer hydrogenation of quinones and related systems: a theoretical mechanistic study

Quantum chemistry calculations have been used to study the uncatalyzed transfer hydrogenation between a range of hydrogen donors and acceptors, in the gas phase and in solution. Our study shows in the first place that in order to obtain reliable condensed-phase transition structures, it is necessary to perform geometry optimization in the presence of a continuum. In addition, the use of a free energy of solvation obtained with the UB3-LYP/6-31+G(d,p)/IEF-PCM/UA0 combination, in conjunction with UMPWB1K/6-311+G(3df,2p)//B3-LYP/6-31+G(d,p) gas-phase energies, gives the best agreement with experimental barriers. In condensed phases, the geometries and energies of the transition structures are found to relate to one another in a manner consistent with the Hammond postulate. There is also a correlation between the barriers and the energies of the radical intermediates in accord with the Bell-Evans-Polanyi principle. We find that in the gas phase, all the transfer-hydrogenation reactions examined proceed via a radical pathway. In condensed phases, some of the reactions follow a radical mechanism regardless of the solvent. However, for some reactions there is a change from a radical mechanism to an ionic mechanism as the solvent becomes more polar. Our calculations indicate that the detection of radical adducts by EPR does not necessarily indicate a predominant radical mechanism, because of the possibility of a concurrent ionic reaction. We also find that the transition structures for these reactions do not necessarily have a strong resemblance to the intermediates, and therefore one should be cautious in utilizing the influence of polar effects on the rate of reaction as a means of determining the mechanism.     

Theory of solvation in polar nematics

We develop a linear response theory of solvation of ionic and dipolar solutes in anisotropic, axially symmetric polar solvents. The theory is applied to solvation in polar nematic liquid crystals. The formal theory constructs the solvation response function from projections of the solvent dipolar susceptibility on rotational invariants. These projections are obtained from Monte Carlo simulations of a fluid of dipolar spherocylinders which can exist both in the isotropic and nematic phases. Based on the properties of the solvent susceptibility from simulations and the formal solution, we have obtained a formula for the solvation free energy which incorporates the experimentally available properties of nematics and the length of correlation between the dipoles in the liquid crystal. The theory provides a quantitative framework for analyzing the steady-state and time-resolved optical spectra and makes several experimentally testable predictions. The equilibrium free energy of solvation, anisotropic in the nematic phase, is given by a quadratic function of cosine of the angle between the solute dipole and the solvent nematic director. The sign of solvation anisotropy is determined by the sign of dielectric anisotropy of the solvent: solvation anisotropy is negative in solvents with positive dielectric anisotropy and vice versa. The solvation free energy is discontinuous at the point of isotropic-nematic phase transition. The amplitude of this discontinuity is strongly affected by the size of the solute becoming less pronounced for larger solutes. The discontinuity itself and the magnitude of the splitting of the solvation free energy in the nematic phase are mostly affected by microscopic dipolar correlations in the nematic solvent. Illustrative calculations are presented for the equilibrium Stokes shift and the Stokes shift time correlation function of coumarin-153 in 4-n-pentyl-4'-cyanobiphenyl and 4,4-n-heptyl-cyanopiphenyl solvents as a function of temperature in both the nematic and isotropic phases.     

Solvent-dependent transition states for decarboxylations

The rate constants and kinetic isotope effects for decarboxylation of 4-pyridylacetic acid depend strongly on whether the solvent is water or dioxane, and the present paper interprets this finding. We calculate the solvent dependence of the free energy barrier and of the (13)C and (18)O kinetic isotope effects using a quantum mechanical solvation model based on class IV charges and semiempirical atomic surface tensions. The calculations provide a consistent interpretation of the experimental results, which provides a striking confirmation of the soundness of the solvation modeling. Even more significantly, the agreement of theory and experiment gives us confidence in the physical picture of the reaction provided by the model. This indicates that the location of the transition state, as measured by the length of the breaking C--C bond, is 0.24 A later than the gas phase in dioxane and 0.37 A later than the gas phase in water. Charge development at the transition state also depends strongly on the solvent; in particular the CO(2) moiety is 0.07 electronic charge units more negative at the transition state in dioxane than in water.     

Hybrid Solvation Model with First Solvation Shell for Calculation of Solvation Free Energy

We present a hybrid solvation model with first solvation shell to calculate solvation free energies. This hybrid model combines the quantum mechanics and molecular mechanics methods with the analytical expression based on the Born solvation model to calculate solvation free energies. Based on calculated free energies of solvation and reaction profiles in gas phase, we set up a unified scheme to predict reaction profiles in solution. The predicted solvation free energies and reaction barriers are compared with experimental results for twenty bimolecular nucleophilic substitution reactions. These comparisons show that our hybrid solvation model can predict reliable solvation free energies and reaction barriers for chemical reactions of small molecules in aqueous solution.     

Two-state model based on the block-localized wave function method

The block-localized wave function (BLW) method is a variant of ab initio valence bond method but retains the efficiency of molecular orbital methods. It can derive the wave function for a diabatic (resonance) state self-consistently and is available at the Hartree-Fock (HF) and density functional theory (DFT) levels. In this work we present a two-state model based on the BLW method. Although numerous empirical and semiempirical two-state models, such as the Marcus-Hush two-state model, have been proposed to describe a chemical reaction process, the advantage of this BLW-based two-state model is that no empirical parameter is required. Important quantities such as the electronic coupling energy, structural weights of two diabatic states, and excitation energy can be uniquely derived from the energies of two diabatic states and the adiabatic state at the same HF or DFT level. Two simple examples of formamide and thioformamide in the gas phase and aqueous solution were presented and discussed. The solvation of formamide and thioformamide was studied with the combined ab initio quantum mechanical and molecular mechanical Monte Carlo simulations, together with the BLW-DFT calculations and analyses. Due to the favorable solute-solvent electrostatic interaction, the contribution of the ionic resonance structure to the ground state of formamide and thioformamide significantly increases, and for thioformamide the ionic form is even more stable than the covalent form. Thus, thioformamide in aqueous solution is essentially ionic rather than covalent. Although our two-state model in general underestimates the electronic excitation energies, it can predict relative solvatochromic shifts well. For instance, the intense pi-->pi* transition for formamide upon solvation undergoes a redshift of 0.3 eV, compared with the experimental data (0.40-0.5 eV).     

Polar and solvation effects in reactions of oxygen atoms and hydroxyl and alkoxyl radicals with oxygen-containing compounds

Experimental data on the activation energies of reactions of H-abstraction from oxygencontaining compounds by oxygen atoms and hydroxyl and alkoxyl radicals in the gas and liquid phases have been analyzed by means of the parabolic model of the transition state. The contribution of polar interaction to the activation energies of the reactions has been calculated. The contribution of solvation to the activation energy has been calculated by comparison of the reaction parameters of the respective reaction in the liquid and gas phases.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 38–42, January, 1994.     

溶液中甲醇和二氯亚砜的化学反应

用B3LYP方法和SCIPCM模型(模拟溶剂效应)研究了甲醇和二氯亚砜在两种非极性(ε<15)和两种极性(ε>15)溶剂中的反应(最终产物为氯代甲烷和二氧化硫). 反应过程由反应(1)和反应(2)组成, CH3OS(O)Cl是反应(1)的主要产物和反应(2)的反应物. 反应(2)有“前面取代”(经过渡态TS3f)和“背后取代”(先经CH3OS(O)Cl的电离, 再经过渡态TS3b)两种机理. 计算表明, 在气相和四种溶剂中反应(1)和(2)都是放热反应, 反应(1)具有相同的反应途径(经过渡态→中间体→过渡态), 溶剂的极性对反应(2)有很大的影响. 在气相和非极性溶剂中, TS3f的能量比(CH3OSO++Cl-)离子对(中间体IM2)的能量低, 反应(2)应为前面取代机理; 在极性溶剂中, IM2和TS3b的能量都比TS3f低, 反应(2)应为背后取代机理.     

Ionic liquid acceleration of solid-phase suzuki-miyaura cross-coupling reactions

[reaction: see text] Room-temperature ionic liquids promote various transition metal-catalyzed reactions in the solution phase. Here, for the first time, we show that these effects are translatable to solid-phase reactions. The Suzuki-Miyaura cross-coupling of 4-iodophenol immobilized on polystyrene-Wang resin with various arylboronic acids was significantly accelerated by the ionic liquid 1-butyl-3-methylimidazolium tetrafluoroborate ([bmim][BF(4)(-)]).