Optimizing the structures of minimum and transition state on the free energy surface

Presented here is the application of a scheme for optimizing the structures of minima and transition states on the free energy surface (FES) for a path along a fixed reaction coordinate with the aid of ab initio molecular dynamics (AIMD) simulation. In the direction of the reaction coordinate, the values corresponding to the stationary points were optimized using the quasi-Newton method, in which the gradient of the free energy along the reaction coordinate was obtained by a constraint AIMD method, and the Bofill Hessian update scheme was used. The equilibrium values for the other directions were taken as the corresponding averages in the dynamic simulation. This scheme was applied to several elementary bimolecular addition reactions: (A) BH(3) + H(2)O --> H(2)O.BH(3); (B) BF(3) + NH(3) --> FB(3).NH(3); (C) SO(3) + NH(3) --> O(3)S.NH(3); (D) C(2)H(4) + CCl(2) --> H(4)C(2).CCl(2); (E) Ni(NH(2))(2) + PH(3) --> (NH(2))(2)Ni.PH(3); (F) W(CO)(5) + CO --> W(CO)(6). For reactions A, B, C, and F, no transition state (TS) exists on the potential energy surface (PES). However there is a TS on the FES. This stems from the curvature difference of the PES and -TDeltaS as a function of the reaction coordinate. For all reactions, it is found that the TS shifts toward the complexation product with increasing temperature because of the curvature increase of -TDeltaS. The equilibrium bond distances for the inactive coordinates perpendicular to the reaction coordinate always increase with temperature, which is due to the thermal excitation and anharmonicity of the PES.     

气相中Nb＋活化CS2中的C-S键

以Nb＋与CS2反应作为第二前过渡金属离子与CS2反应生成金属硫化物离子和CS的范例体系．采用密度泛函UB3LYP方法,对于Nb＋采用Stuttgart赝势基组,对于C和S采用6—311＋G（2d）基组,计算研究了Nb＋在基态和激发态时与CS2气相反应的机理．全参数优化了反应势能面上各驻点的儿何构型,并且用频率分析方法和内禀反应坐标方法对过渡态进行了验证．结果表明Nb＋与CS2的反应是插入-消去反应,在反应过程中会发生系间窜越,并且找到了两个势能面的能量最低交叉点．     

A new mechanism for the production of highly vibrationally excited OH in the mesosphere: an ab initio study of the reactions of O2(A 3Sigmau+ and A' 3Deltau)+H

In an attempt to explain the observed nightglow emission from OH(v=10) in the mesosphere that has the energy greater than the exothermicity of the H+O(3) reaction, potential energy surfaces were calculated for reactions of high lying electronic states of O(2)(A (3)Sigma(u) (+) and A' (3)Delta(u)) with atomic hydrogen H((2)S) to produce the ground state products OH((2)Pi)+O((3)P). From collinear two-dimensional scans, several adiabatic and nonadiabatic pathways have been identified. Multiconfigurational single and double excitation configuration interaction calculations show that the adiabatic pathways on a (4)Delta potential surface from O(2)(A' (3)Delta)+H and a (4)Sigma(+) potential surface from O(2)(A (3)Sigma(u) (+))+H are the most favorable, with the zero-point corrected barrier heights of as low as 0.191 and 0.182 eV, respectively, and the reactions are fast. The transition states for these pathways are collinear and early, and the reaction coordinate suggests that the potential energy release of ca. 3.8 eV (larger than the energy required to excite OH to v=10) is likely to favor high vibrational excitation.     

Guided ion beam and theoretical study of the reactions of Hf+ with H2, D2, and HD

The kinetic energy dependences of reactions of the third-row transition metal cation Hf(+) with H(2), D(2), and HD were determined using a guided ion beam tandem mass spectrometer. A flow tube ion source produces Hf(+) in its (2)D (6s(2)5d(1)) electronic ground state level. Corresponding state-specific reaction cross sections are obtained. The kinetic energy dependences of the cross sections for the endothermic formation of HfH(+) and HfD(+) are analyzed to give a 0 K bond dissociation energy of D(0)(Hf(+)-H)=2.11±0.08 eV. Quantum chemical calculations at several levels of theory performed here generally overestimate the experimental bond energy but results obtained using the Becke-half-and-half-LYP functional show good agreement. Theory also provides the electronic structures of these species and the reactive potential energy surfaces. Results from the reactions with HD provide insight into the reaction mechanisms and indicates that Hf(+) reacts via a statistical mechanism. We also compare this third-row transition metal system with the first-row and second-row congeners, Ti(+) and Zr(+), and find that Hf(+) has a weaker M(+)-H bond. As most third-row transition metal hydride cation bonds exceed their lighter congeners, this trend is unusual but can be understood using promotion energy arguments.     

气相中Cu＋和Zn＋与SCO反应的理论研究

为更清晰地揭示M^＋与SCO基元反应的机理, 采用密度泛函B3LYP方法, 在6-311＋＋G**基组水平上研究了Cu^＋＋SCO和Zn^＋＋SCO反应体系. 对反应势能面上各驻点的几何构型进行了全优化, 用频率分析方法和内禀反应坐标(IRC)方法对过渡态进行了验证. 在Cu^＋与SCO的反应中, 对影响反应机理和反应速率的势能面交叉现象进行了讨论, 运用Hammond假设和Yoshizawa等的内禀反应坐标垂直激发的计算方法找到了势能面交叉点. 计算结果表明, C—S和C—O键的活化都是通过插入消去机理, 但C—S键的活化在能量上更占优势. 计算确认了标题反应的主通道, 所有的计算结果与实验吻合.     

Semiclassical nonadiabatic dynamics based on quantum trajectories for the O(3P,1D) + H2 system

The O(3P,1D) + H2 --> OH + H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states.     

激发态Ti＋(3d14s2)与丙炔醇气相反应理论研究

用量子化学密度泛函(DFT)方法研究了激发态Ti^＋(3d¹4s²)与丙炔醇(PPA)气相反应的机理. 在B3LYP/DZVP水平上, 优化了反应的两个通道的反应物、中间体、过渡态和产物的几何构型, 并在MP4/[6-311＋G**(C,H,O)＋Lanl2dz (Ti)]水平上计算了各驻点的单点能量. 为了确证过渡态的真实性, 在B3LYP/DZVP水平上进行了内禀坐标(IRC)计算和频率分析, 获得了二重态反应势能面, 确定了反应机理. 研究结果表明生成产物为[C₃H₃O]^＋和Ti—H的通道是主要反应途径.     

Reaction of bare VO+ and FeO+ with ammonia: a theoretical point of view

The potential energy surfaces corresponding to the dehydration reaction of NH(3) by VO(+) ((3)Sigma, (1)Delta, (5)Sigma) and FeO(+) ((6)Sigma, (4)Delta) metal oxide cations have been investigated within the framework of the density functional theory in its B3LYP formulation and by employing new optimized basis sets for iron and vanadium. The reaction is proposed to occur through two hydrogen shifts from the nitrogen to the oxygen atom giving rise to multicentered transition states. Possible spin crossing between surfaces at different spin multiplicities has been considered. The energy profiles are compared with the corresponding ones for the insertion of bare cations to investigate the influence on reactivity of the presence of the oxygen ligand. The topological analysis of the gradient field of the electron localization function has been used to characterize the nature of the bonds for all the minima and transition states along the paths.     

Probes of spin conservation in heavy metal reactions: experimental and theoretical studies of the reactions of Re+ with H2, D2, and HD

A guided ion beam tandem mass spectrometer is used to examine the kinetic energy dependence of reactions of the third-row transition metal cation, Re(+), with molecular hydrogen and its isotopologues. A flow tube ion source produces Re(+) in its (7)S(3) electronic ground state. Reaction with H(2), D(2), and HD forms Re H(+)(Re D(+)) in endothermic processes. Modeling of the endothermic reaction cross sections yields the 0 K bond dissociation energy of D(0)(Re(+)-H)=2.29+/-0.07 eV (221+/-6 kJ/mol). The experimental thermochemistry is consistent with ab initio calculations, performed here and in the literature. Theory also provides the electronic structures of these species and is used to examine the reactive potential energy surfaces. Results from reactions with HD provide insight into the reaction mechanisms and indicate that the late metal ion, Re(+), reacts largely via a statistical mechanism. This is consistent with the potential energy surfaces which locate a stable Re H(2) (+)((5)B(2)) complex. Results for this third-row transition metal system are compared with the first-row congener (Mn(+)) and found to have much higher reactivity towards dihydrogen and stronger M(+)-H bonds. These differences can be attributed to efficient coupling among surfaces of different spin along with lanthanide contraction and relativistic effects.     

Theoretical Study on Photoisomerization of 2-Methylfuran to 3-Methylfuran

1 INTRODUCTION 2-Methylfuran belongs to the basic heteroaromatic compounds relevant to many fields of modern che- mistry, ranging from the study of natural products and biologically active substances to the develop- ment of building blocks for organic synthesis and conducting polymers[1]. Since the photochemistry ofR-furan was gradually recognized in 1960s[2～7], lots of interest has been aroused. Herein we only study one branch of photoche- mistry of R-furan: the isomerization of 2-methy…     

Electronic structure and solvation of copper and silver ions: a theoretical picture of a model aqueous redox reaction

Electronic states and solvation of Cu and Ag aqua ions are investigated by comparing the Cu(2+) + e(-)--> Cu(+) and Ag(2+) + e(-)--> Ag(+) redox reactions using density functional-based computational methods. The coordination number of aqueous Cu(2+) is found to fluctuate between 5 and 6 and reduces to 2 for Cu(+), which forms a tightly bound linear dihydrate. Reduction of Ag(2+) changes the coordination number from 5 to 4. The energetics of the oxidation reactions is analyzed by comparing vertical ionization potentials, relaxation energies, and vertical electron affinities. The model is validated by a computation of the free energy of the full redox reaction Ag(2+) + Cu(+) --> Ag(+) + Cu(2+). Investigation of the one-electron states shows that the redox active frontier orbitals are confined to the energy gap between occupied and empty states of the pure solvent and localized on the metal ion hydration complex. The effect of solvent fluctuations on the electronic states is highlighted in a computation of the UV absorption spectrum of Cu(+) and Ag(+).     

CH3NO2+X(X=H,OH, CH3, CH2[^3B1], O[^3P])→CH2NO2+XH吸氢反应过渡态 结 构和反应位垒的DFT计算研究

采用DFT(B3LYP)方法,分别在6-311g(d,p),6-311++g(d,p)和自洽相关基组cc-pVIZ水平上优化了基态硝基甲烷和自由基H,OH,CH3,CH2[^3B1]以及O[^3P]等发生吸氢反应时的过渡态结构,并计算了反应的位垒。研究表明,对同一反应,不同基组下优化得到的过渡态几何结构基本一致;反应位垒数值的大小也基本接近,经校正,硝基甲烷同自由基反应位垒的理论计算值同实验结果基本吻合。     

Experimental and theoretical study of the photodissociation reaction of thiophenol at 243 nm: intramolecular orbital alignment of the phenylthiyl radical

The photoinduced hydrogen (or deuterium) detachment reaction of thiophenol (C(6)H(5)SH) or thiophenol-d(1) (C(6)H(5)SD) pumped at 243 nm has been investigated using the H (D) ion velocity map imaging technique. Photodissociation products, corresponding to the two distinct and anisotropic rings observed in the H (or D) ion images, are identified as the two lowest electronic states of phenylthiyl radical (C(6)H(5)S). Ab initio calculations show that the singly occupied molecular orbital of the phenylthiyl radical is localized on the sulfur atom and it is oriented either perpendicular or parallel to the molecular plane for the ground (B(1)) and the first excited state (B(2)) species, respectively. The experimental energy separation between these two states is 2600+/-200 cm(-1) in excellent agreement with the authors' theoretical prediction of 2674 cm(-1) at the CASPT2 level. The experimental anisotropy parameter (beta) of -1.0+/-0.05 at the large translational energy of D from the C(6)H(5)SD dissociation indicates that the transition dipole moment associated with this optical transition at 243 nm is perpendicular to the dissociating S-D bond, which in turn suggests an ultrafast D+C(6)H(5)S(B(1)) dissociation channel on a repulsive potential energy surface. The reduced anisotropy parameter of -0.76+/-0.04 observed at the smaller translational energy of D suggests that the D+C(6)H(5)S(B(2)) channel may proceed on adiabatic reaction paths resulting from the coupling of the initially excited state to other low-lying electronic states encountered along the reaction coordinate. Detailed high level ab initio calculations adopting multireference wave functions reveal that the C(6)H(5)S(B(1)) channel may be directly accessed via a (1)(n(pi),sigma(*)) photoexcitation at 243 nm while the key feature of the photodissociation dynamics of the C(6)H(5)S(B(2)) channel is the involvement of the (3)(n(pi),pi(*))-->(3)(n(sigma),sigma(*)) profile as well as the spin-orbit induced avoided crossing between the ground and the (3)(n(pi),sigma(*)) state. The S-D bond dissociation energy of thiophenol-d(1) is accurately estimated to be D(0)=79.6+/-0.3 kcalmol. The S-H bond dissociation energy is also estimated to give D(0)=76.8+/-0.3 kcalmol, which is smaller than previously reported ones by at least 2 kcalmol. The C-H bond of the benzene moiety is found to give rise to the H fragment. Ring opening reactions induced by the pi-pi(*)n(pi)-pi(*) transitions followed by internal conversion may be responsible for the isotropic broad translational energy distribution of fragments.     

Au+(1S, 3D)与N2O(1Σ+)反应机理的理论研究

采用密度泛函B3LYP方法, O和N用6-311+G*基组, Au+用赝势基组(8s7p6d)/[6s5p3d], 研究了Au+(1S, 3D)离子和N2O(1Σ+)分子的反应机理. 报道了在基态单重态和激发三重态势能面上各反应物、中间体和过渡态的构型特征及能量. 结果表明, 两个主反应通道Au+(1S)+ N2O(1Σ+)→1NA-Complex-1→1NA-TS1→1NA-Complex-2→1NA-Crossing→[3OAuNN]+和Au+(1S)+ N2O(1Σ+)→1NB-Complex→1NB-Crossing→[AuNN(1Σ+)]++O(3P)都需经过反应交叉势能面, 出现“系间窜越”. 用内禀坐标单点计算垂直激发态的方法确定了势能面交叉点, 并用含时密度泛函TD-B3LYP方法进一步探讨了自旋翻转机理.     

几个硫-叶立德反应机理的量子拓扑研究

采用MP2(FC)/6-311++G(d,p)对硫叶立德和类硫叶立德自由基反应机理进行了探讨。优化了中间体、过滤态和产物的几何构型。本文侧重从量子拓扑学的角度,对IRC(内禀反应坐标)反应进程中各点进行电子密度拓扑分析,讨论了反应过程中化学键的断裂、生成和化学键的变化规律。找到了这类反应的能量过渡态和结构过渡态,上述两个反应都是先经历一个没有形成三元环拓扑结构的能量过渡态,再经历一个形成了三元环拓扑结构的结构过渡态,最后到达产物。     

Theoretical Study on the Mechanism of the Gas-phase Reaction of Sc＋ with Propargyl Alcohol

In order to elucidate the reaction mechanisms of reaction Sc^＋ with propargyl alcohol (PPA), the triplet potential energy surface for the reactions has been theoretically investigated using a DFT method. The geometries for the reactants, intermediates, transition states and products were completely optimized at B3LYP/DZVP level. The single point energy of each stationary point was calculated at MP4/(6-311＋G** for C, H, O and Lanl2dz for Sc^＋) level. All the transition states were verified by the vibrational analysis and the internal reaction coordinate (IRC) calculations. The present results show that the reaction takes an insertion-elimination mechanism both along the O—H and C—O bond activation branches, but the C—O bond activation is much more favorable in energy than the O—H bond activation. All theoretical results not only support the existing conclusions inferred from early experiment, but also complement the pathway and mechanism for this reaction.     

Mechanisms of the reactions of W AND W+ with H2O: computational studies

The mechanism of the reactions of W and W(+) with the water molecule have been studied for several lower-lying electronic states of tungsten centers at the CCSD(T)/6-311G(d,p)+SDD and B3LYP/6-31G(d,p)+SDD levels of theory. It is shown that these reactions are essentially multistate processes, during which lower-lying electronic states of the systems cross several times. They start with the formation of initial prereaction M(H(2)O) complexes with M-H(2)O bonding energies of 9.6 and 48.2 kcal/mol for M = W and W(+), followed by insertion of the metal center into an O-H bond with 20.0 and 53.3 kcal/mol barriers for neutral and cationic systems, respectively. The overall process of M + H(2)O --> t-HM(OH) is calculated to be highly exothermic, 48.4 and 48.8 kcal/mol for M = W and W(+). From the HM(OH) intermediate the reaction may proceed via several different channels, among which the stepwise HM(OH) --> HMO + H --> (H)(2)MO and concerted HM(OH) --> (H)(2)MO pathways are more favorable and can compete (energetically) with each other. For the neutral system (M = W), the concerted process is the most favorable, whereas for the charged system (M = W(+)), the stepwise pathway is slightly more favorable. From the energetically most favorable intermediate (H)(2)MO the reactions proceed via H(2)-molecule formation with a 53.1 kcal/mol activation barrier for the neutral system. For the cationic system, H-H formation and dissociation is an almost barrierless process. The overall reaction of W and W(+) with the water molecule leading to H(2) + MO formation is found to be exothermic by 48.2 and 39.8 kcal/mol, respectively. In the gas phase with the collision-less conditions the reactions W((7)S) + H(2)O --> H(2) + WO((3)Sigma(+)), and W(+)((6)D) + H(2)O --> H(2) + WO(+)((4)Sigma(+)) are expected to proceed via a 10.4 and 5.1 kcal/mol overall energy barrier corresponding to the first O-H dissociation at the TS1. On the basis of these PESs, we predict kinetic rate constants for the reactions of W and W(+) with H(2)O.     

邻位吡啶自由基多通道分解反应的动力学和反应机理研究

Utilizing Gaussian94 program package, all species involved in decomposition reactions of o-pyridyl radical were optimized fully at B3LYP/6-311＋＋G^＊＊ level. Intrinsic reaction coordinate calculations were employed to confirm the connections of the transition states and products, and transition states were ascertained by the number of imaginary frequency （0 or 1）. The reaction mechanism was elucidated by the vibrational mode analysis and electronic population analysis, and the reaction rate constants were calculated with transition state theory.