A reaction density functional theory study of solvent effect in the nucleophilic addition reactions in aqueous solution

Whereas the proper choice of reaction solvent constitutes the cornerstone of the green solvent concept,solvent effects on chemical reactions are not mechanistically well understood due to the lack of feasible molecular models.Herein,by taking the case study of nucleophilic addition reaction in aqueous solution,we extend the proposed multiscale reaction density functional theory(RxDFT)method to investigate the intrinsic free energy profile and total free energy profile,and study the solvent effect on the activation and reaction free energy for the nucleophilic addition reactions of hydroxide anion with methanal and carbon dioxide in aqueous solution.The predictions of the free energy profile in aqueous solution for these two nucleophilic addition reactions from RxDFT have a satisfactory agreement with the results from the RISM and MD-FEP simulation.Meanwhile,the solvent effect is successfully addressed by examining the difference of the free energy profile between the gas phase and aqueous phase.In addition,we investigate the solvent effect on the reactions occurred near solid-liquid interfaces.It is shown that the activation free energy is significantly depressed when reaction takes place in the region within 10A distance to the substrate surface owing to the decrease of hydration free energy at the solid-liquid interface.

A reaction density functional theory study of solvent effects on keto-enol tautomerism and isomerization in pyruvic acid

It is important to study the solvent effect on keto-enol tautomerism that has applications in many areas of chemical engineering.In this work,we use a multiscale reaction density functional theory(Rx DFT)to study the keto-enol tautomerism and isomerization of pyruvic acid.The results show that both effects of solvation and water assistance could reduce the reaction barriers.The water molecule participates the reaction as a catalyst to accept/give the protons with forming a hexagonal ring in the transition state.As a result of this temporary and intermediate hexagonal ring,the solute configuration undergoes a small variation during the reaction,giving a diminished contribution to the intrinsic reaction free energy.The solvent distribution shows a local ordering behavior near the solute that also reduces the contribution of solvation effect to the reaction barrier.Water assistance plays a major role in both pre-reaction and postreaction process.In terms of the driving force for the reaction,the effects of both solvation and water assistance are important.

Nonequilibrium Solvent Free Energy Curve from Molecular Theory in Electron Transfer Reaction

The microscopic moleeular theory for electron transfer in a model solvent ishahr developed. The nonlinear response of the solvent molecules is be computedquanitatively in a new way. Adopting computer simulation daa and choosingappropriate reaction coordinae, a reasonable free energy dinram is constructed and thercorganhaion energy for the product state is calculated.

Solvent effects on Diels-Alder reaction in ionic liquids:A reaction density functional study

Extensive experimental studies have been performed on the Diels-Alder(DA)reactions in ionic liquids(ILs),which demonstrate that the IL environment can significantly influence the reaction rates and selectivity.However,the underlying microscopic mechanism remains ambiguous.In this work,the multiscale reaction density functional theory is applied to explore the effect of 1-butyl-3-methylimidazolium hexafluorophosphate([BMIM][PF_(6)])solvent on the reaction of cyclopentadiene(CP)with acrolein,methyl acrylate,or acrylonitrile.By analyzing the free energy landscape during the reaction,it is found that the polarization effect has a relatively small influence,while the solvation effect makes both the activation free energy and reaction free energy decrease.In addition,the rearrangement of local solvent structure shows that the cation spatial distribution responds more evidently to the reaction than the anion,and this indicates that the cation plays a dominant role in the solvation effect and so as to affect the reaction rates and selectivity of the DA reactions.

Theoretical analysis of the double-differential cross-sections of neutron,proton,deuteron,^(3)He,andαfor the p+^(6) Li reaction

Based on the unified Hauser–Feshbach and exciton model,which can describe the particle emission processes between discrete energy levels with energy,angular momentum,and parity conservations,a statistical theory of light nucleus reaction(STLN)is developed to calculate the double-differential cross-sections of the outgoing neutron and light charged particles for the proton-induced^(6) Li reaction.A significant difference is observed between the p+^(6) Li and p+^(7) Li reactions owing to the discrepancies in the energy-level structures of the targets.The reaction channels,including sequential and simultaneous emission processes,are analyzed in detail.Taking the double-differential cross-sections of the outgoing proton as an example,the influence of contaminations(such as^(1) H,^(7)Li,^(12)C,and^(16)O)on the target is identified in terms of the kinetic energy of the first emitted particles.The optical potential parameters of the proton are obtained by fitting the elastic scattering differential cross-sections.The calculated total double-differential cross-sections of the outgoing proton and deuteron at E_(p)=14 MeV agree well with the experimental data for different outgoing angles.Simultaneously,the mixed double differential cross-sections of^(3) He andαare in good agreement with the measurements.The agreement between the measured data and calculated results indicates that the two-body and three-body breakup reactions need to be considered,and the pre-equilibrium reaction mechanism dominates the reaction processes.Based on the STLN model,a PLUNF code for the p+^(6) Li reaction is developed to obtain an ENDF-6-formatted file of the double-differential cross-sections of the nucleon and light composite charged particles.

Energy diffusion controlled reaction rate in dissipative Hamiltonian systems

In this paper the energy diffusion controlled reaction rate in dissipative Hamiltonian systems is investigated by using the stochastic averaging method for quasi Hamiltonian systems. The boundary value problem of mean first- passage time （MFPT） of averaged system is formulated and the energy diffusion controlled reaction rate is obtained as the inverse of MFPT. The energy diffusion controlled reaction rate in the classical Kramers bistable potential and in a two-dimensional bistable potential with a heat bath are obtained by using the proposed approach respectively. The obtained results are then compared with those from Monte Carlo simulation of original systems and from the classical Kraraers theory. It is shown that the reaction rate obtained by using the proposed approach agrees well with that from Monte Carlo simulation and is more accurate than the classical Kramers rate.

Mixed Quantum Classical Reaction Rates based on the Phase Space Formulation of the Hierarchical Equations of Motion

In a previous work[J.Chem.Phys.140,174105(2014)],we have shown that a mixed quantum classical(MQC)rate theory can be derived to investigate the quantum tunneling effects in the proton transfer reactions.However,the method is based on the high temperature approximation of the hierarchical equation of motion(HEOM)with the Debye-Drude spectral density,and results in a multistate Zusman type of equation.We now extend this theory to include quantum effects of the bath degrees of freedom.By writing the full HEOM into a multidimensional partial differential equation in phase space,we can define a new reaction coordinate,and the previous method can be generalized to the full quantum regime.The validity of the new method is demonstrated by using numerical examples,including the spin-Boson model,and the double well model for proton transfer reaction.The new method is found to resolve some key problems of the previous theory based on high temperature approximation,including possible numerical instability in long time simulation and wrong rate constant at low temperatures.

Dissolution Precipitation Wave Structure of Hydrothermal Ore Zoning

Hydrothermal ore zoning is a transport-reaction problem in which infiltration is the principal Prcness of transport and dissolution/Precipitation is the Principal process of chemical reactions.Neglecting diffusion and ion exchange/adsorption would not affect the basic attributes of hydrothermal ore zoning. Hydrothermal ore zoning belongs essentially to infiltration metasomatic zoning, it results from the formation and propagation of dissolution/precipitation waves through Permeable media. The authors apply the theory of coupled infiltration and dissolution/precipitation reactions in Physicochemical hydrodynamics to studying the structural characteristics of dissolution/precipitation waves, and apply furthermore the coherence principle in dynamic theory of multicomponent coupled systems to revealing the dynamic mechanisms of their formation. The results of investigation verify and develop . C. 's theory of infiltration metasomatic zoning,on the one hand, raising it from the qualitative, equilibrium thermodynamic basis to the quantitative dynamic level;on the other hand, and more importantly, applying theories of Physicochemical hydrodynamics and dynamics of multicomponent coupled systems to bringing to light the dynamic mechanisms of formation of the structure of hydrothermal ore zoning, and advancing a theory of hydrothermal ore zoning, putting forward new ideas on the nature of the problem of hydrothermal ore zoning, the essence of hydrothermal ore zoning and the structural characteristics and mechanisms of formation of hydrothermal ore zoning.

Variational transition state theory study of the reactions Li+HF and Li+HCl on the BO potential energy surfaces

Variational transition state theory is used to study the kinetics of the reactions Li+HF and Li+HCl on the BO potential energy surfaces, fitted to ab initio values by Lagana et al. The saddle point, the location of generalized transition state and other kinetic features on the bottleneck of the po- tential energy surface are obtained. We point out that the location of generalized transition state tends to the direction in which the two atoms with larger reduced mass separate off. It can be determined whether a potential energy surface is accurate or not by a fine analysis of its bottleneck region and a comparison of kinetic calculation with experimental result. This is helpful for building a more accurate potential energy surface.

Hydrogenation of graphene nanoflakes and C-H bond dissociation of hydrogenated graphene nanoflakes： a density functional theory study

The Gibbs free energy change for the hydro- genation of graphene nanoflakes Cn （n = 24, 28, 30 and 32） and the C-H bond dissociation energy of hydrogenated graphene nanoflakes CnHm （n = 24, 28, 30 and 32; and m = 1, 2 and 3） are evaluated using density functional theory calculations. It is concluded that the graphene nanoflakes and hydrogenated graphene nanoflakes accept the orth- aryne structure with peripheral carbon atoms bonded via the most triple bonds and leaving the least unpaired dan- gling electrons. Five-membered rings are formed at the deep bay sites attributing to the stabilization effect from the pairing of dangling electrons. The hydrogenation reactions which eliminate one unpaired dangling electron and thus decrease the overall multiplicity of the graphene nanoflakes or hydrogenated graphene nanoflakes are spontaneous with negative or near zero Gibbs free energy change. And the resulting C-H bonds are stable with bond dissociation energy in the same range as those of aromatic compounds. The other C-H bonds are not as stable attributing to the excessive unpaired dangling electrons being filled into the C-H anti-bond orbital.