含有一个非平面杂环胺配体的新型反铂抗癌药物的水解机理

用B3LYP杂化泛函和等电子聚焦极化连续模型(IEF-PCM)研究了trans-[PtCl2(NH3)(Am)](Am: 非平面哌啶或哌嗪)新型反铂抗癌药物的水解反应机理. 对经由一般的SN2机理的第一步和第二步水解反应势能面上的稳定点进行了全优化和表征. 在水解中, 最显著的结构变化发生在反应过渡态和中间体的五配位三角双锥的赤道面上. 与经典顺铂(cisplatin)比较, 反式[PtCl2(NH3)(piperazine)]的第一步和第二步水解活化能均低于顺铂, 而反式[PtCl2(NH3)(piperidine)]的第一步水解活化能稍高于顺铂, 第二步水解活化能稍低于顺铂. 计算表明, 这些含有非平面杂环胺反铂的配合物减小了赤道面上的立体效应和水解势垒.     

非对称反铂抗癌药物反-异丙胺·间羟甲基吡啶·二氯铂水解机理的理论研究

应用密度泛函理论PBE1PBE方法及CPCM模型计算具有空间位阻的非对称反铂抗癌药反-异丙胺·间羟甲基吡啶·二氯铂的水解反应机理.结果表明,由于空间效应,水分子从垂直于Pt平面四边形配位的方向进攻,其水解反应为水的H,O原子分别与Cl,Pt原子形成平面四边形结构的协同作用结果,Pt的5d电子和Cl的3p电子分别向水的H—O反键轨道离域,O的孤对电子向Pt的价层空轨道或Pt—Cl反键轨道离域,速率决定步骤经过一个近似三角双锥的过渡态完成.随着反应的进行,离域效应增强,Pt与O作用增强,而Pt—Cl键减弱.溶剂化效应使两步水解反应的各反应物、生成物和过渡态的能量较气相时低63.6～386.3kJ/mol,单点能垒较气相反应低约17.1～36.2kJ/mol.从空间位阻较小的异丙基相反方向进攻的反应通道更易进行,其中1B和2B通道活化焓(分别为79.7和87.8kJ/mol)最小,是第一、二步水解反应的主要通道.第二步水解活化能垒远高于顺铂,两步水解活化能垒均高于对称的反铂trans-[PtCl2(i-pra)2].     

DFT studies of the phosphinidene derivative HPNaF and its insertion reaction with R–H(R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>)

The formations of the phosphinidene derivative HPNaF and its insertion reactions with R–H (R=F, OH, NH₂, CH₃) have been systematically investigated employing the density functional theory (DFT), such as the B3LYP and MPW1PW91 methods. A comparison with the results of MP2 calculations shows that MPW1PW91 underestimates the barrier heights for the four reactions considered. Similarly, the same is also true for the B3LYP method depending on the selected reactions, but by much less than MPW1PW91, where the barrier heights of the four reactions are 25.2, 85.7, 119.0, and 142.4 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. All the mechanisms of the four reactions are identical to each other, i.e., an intermediate has been located during the insertion reaction. Then, the intermediate could dissociate to substituted phosphinidane(H₂RP) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are −92.2, −68.1, −57.2, and −44.3 kJ/mol at the B3LYP/6-311+G* level of theory, respectively, where both the B3LYP and MPW1PW91 methods underestimate the reaction energies compared with the MP2 results. The linear correlations between the calculated barrier heights and the reaction energies have also been observed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Probing the imine silylenoid HN=SiNaF and its insertions reaction with R–H (R=F,OH, NH<Subscript>2</Subscript>, CH<Subscript>3</Subscript>) using DFT

The geometries and isomerization of the imine silylenoid HN=SiNaF as well as its insertion reactions with some R–H molecules have been systematically investigated theoretically, where R=F, OH, NH₂, and CH₃, respectively. The barrier heights for the four insertion reactions are 67.7, 115.6, 153.5, and 271.5 kJ/mol at the B3LYP/6-311+G* level of theory, respectively. Here, all the mechanisms of the four reactions are identical to each other, i.e., a stable intermediate has been formed during the insertion reaction. Then, the intermediate could dissociate into the substituted silylene (HN=SiHR) and NaF with a barrier corresponding to their respective dissociation energies. Correspondingly, the reaction energies for the four reactions are 71.8, 95.5, 123.3, and 207.6 kJ/mol, respectively, which are linearly correlated with the calculated barrier heights. Furthermore, the effects of halogen substitutions (F, Cl, and Br) on the reaction activity have also been discussed. As a result, the relative reactivity among the four insertion reactions should be as follows: H–F > H–OH > H–NH₂ > H–CH₃.     

Direct determination of the ionization energies of PtC, PtO, and PtO2 with VUV radiation

Photoionization efficiency curves were measured for gas-phase PtC, PtO, and PtO2 using tunable vacuum ultraviolet (VUV) radiation at the Advanced Light Source. The molecules were prepared by laser ablation of a platinum tube, followed by reaction with CH4 or N2O and supersonic expansion. These measurements provide the first directly measured ionization energy for PtC, IE(PtC) = 9.45 +/- 0.05 eV. The direct measurement also gives greatly improved ionization energies for the platinum oxides, IE(PtO) = 10.0 +/- 0.1 eV and IE(PtO2) = 11.35 +/- 0.05 eV. The ionization energy connects the dissociation energies of the neutral and cation, leading to greatly improved 0 K bond dissociation energies for the neutrals: D0(Pt-C) = 5.95 +/- 0.07 eV, D0(Pt-O) = 4.30 +/- 0.12 eV, and D0(OPt-O) = 4.41 +/- 0.13 eV, as well as enthalpies of formation for the gas-phase molecules DeltaH(0)(f,0)(PtC(g)) = 701 +/- 7 kJ/mol, DeltaH(0)(f,0)(PtO(g)) = 396 +/- 12 kJ/mol, and DeltaH(0)(f,0)(PtO2(g)) = 218 +/- 11 kJ/mol. Much of the error in previous Knudsen cell measurements of platinum oxide bond dissociation energies is due to the use of thermodynamic second law extrapolations. Third law values calculated using statistical mechanical thermodynamic functions are in much better agreement with values obtained from ionization energies and ion energetics. These experiments demonstrate that laser ablation production with direct VUV ionization measurements is a versatile tool to measure ionization energies and bond dissociation energies for catalytically interesting species such as metal oxides and carbides.     

A DFT study of H-isomerisation in alkoxy-, alkylperoxy- and alkyl radicals: Some implications for radical chain reactions in polymer systems

Intramolecular 1-n H-shift (n = 2, 3… 7) reactions in alkoxy, alkyl and peroxy radicals were studied by density functional theory (DFT) at the B3LYP/6-311+G∗∗ level and compared with respective intermolecular H-transfers. It was found that starting from 1 to 3 H-shift the barrier heights stepwise decrease with increasing n reaching minimum for 1-5 and 1-6 H-shifts. This dependence can be ascribed to the decrease of the strain with increasing transition state (TS) ring size, which is minimal in six- and seven-member rings. The barrier heights of H-shifts in alkyl radicals are systematically larger than those in alkoxy radicals: the respective activation energies (E_a) of 1-5 and 1-6 H-shifts are about 59-67 kJ/mol for alkyl radical and 21-34 kJ/mol for alkoxy radicals. Further increase of the TS ring size in 1-7 H-shifts leads to the increase of the barrier to 44 kJ/mol in the hexyloxy radical and 84 kJ/mol for n-heptyl radical. We have also found that intermolecular H-transfer reactions in all three types of free radicals have smaller barriers than respective intramolecular 1-5 or 1-6 H-shifts by 4-25 kJ/mol. The mentioned difference can be explained in terms of enhanced nonbonding repulsion interaction in the cyclic TS structures compared to respective intermolecular TS. B3LYP/6-311+G∗∗ geometric parameters and imaginary frequencies for 1-n H-shifts TS are consistent with respective calculated barrier heights. Reactivity of some other radicals compared to alkoxy, peroxy and alkyl radicals as well as other factors influencing their reactivity (π-conjugation, steric effect and ring strain in cyclic TS, etc.) are also briefly discussed in relation to free radical reactions in polymer systems.     

Mechanistic Study on the Monofunctional Binding of Hydrolysis Products of Non-planar Heterocyclic Complexes Trans-[PtCl_2NH_3（piperidine）] and Trans-[PtCl_2（piperidine）_2] to DNA Purine Bases

The monofunctional substitution reactions between trans-[PtCl(H2O)(NH3)(pip)]+,trans-[Pt(H2O)2(NH3)(pip)]2+,trans-[PtCl(H2O)(pip)2]+,trans-[Pt(H2O)2(pip)2]2+ (pip = piperidine) and adenine/guanine nucleotides are explored by using B3LYP hybrid functional and IEF-PCM salvation models. For the trans-[Pt(H2O)2(NH3)(pip)]2+ and trans-[PtCl(H2O)(NH3)(pip)]+ complexes,the computed barrier heights in aqueous solution are 13.5/13.5 and 11.6/11.6 kcal/mol from trans-Pt-chloroaqua complex to trans/cis-monoadduct for adenine and guanine,and the corresponding values are 20.7/20.7 and 18.8/18.8 kcal/mol from trans-Pt-diaqua complex to trans/cis-monoadduct for adenine and guanine,respectively. For trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+,the corresponding values are 21.5/21.3 and 19.4/19.4 kcal/mol,and 26.0/26.0 and 20.7/20.8 kal/mol for adenine and guanine,respectively. Our calculations demonstrate that the barrier heights of chloroaqua are lower than the corresponding values of diaqua for adenine and guanine. In addition,the free energies of activation for guanine in aqueous solution are all smaller than that for adenine,which predicts a preference of 1.9 kcal/mol when trans-[PtCl(H2O)(NH3)(pip)]+ and trans-[Pt(H2O)2(NH3)(pip)]2+ are the active agents and ～1.9 and ～ 5.3 kcal/mol when trans-[PtCl(H2O)(pip)2]+ and trans-[Pt(H2O)2(pip)2]2+ are the active agents,respectively. For the reaction of trans-Pt-chloroaqua (or diaqua) to cis-monoadduct,we obtain the same transition-state structure as from the reaction of trans-Pt-chloroaqua (or diaqua) to trans-monoadduct,which seems that the trans-Pt-chloroaqua (or diaqua) complex can generate trans-or cis-monoadduct via the same transition-state.     

Iridium-imine and -amine complexes relevant to the (S)-metolachlor process: structures, exchange kinetics, and C-H activation by Iri causing racemization

Iridium complexes of DMA-imine [2,6-dimethylphenyl-1'-methyl-2'-methoxyethylimine, 1 a) and (R)-DMA-amine [(1'R)-2,6-dimethylphenyl-1'-methyl-2'-methoxyethylamine, 2 a] that are relevant to the catalytic imine hydrogenation step of the Syngenta (S)-Metolachlor process were synthesized: metathetical exchange of [Ir2Cl2(cod)2] (cod=1,5-cyclooctadiene) with [Ag(1 a)2]BF4 and [Ag((R)-2 a)2]BF4 afforded [Ir(cod)(kappa2- -1 a)]BF4 (11) and [Ir(cod)(kappa2-(R)-2 a)]BF4 ((R)-19)), respectively. These complexes were then used in stopped-flow experiments to study the displacement of amine 2 a from complex 19 by imine 1 a to form the imine complex 11, thus modeling the product/substrate exchange step in the catalytic cycle. The data suggest a two-step associative mechanism characterized by k1=(2.6+/-0.3) x 10(2) M(-1) s(-1) and k2=(4.3+/-0.6) x 10(-2) s(-1) with the respective activation energies EA1=(7.5+/-0.6) kJ mol(-1) and EA2=(37+/-3) kJ mol(-1). Furthermore, complex 11 reacted with H2O to afford the hydrolysis product [Ir(cod)(eta(6-)-2,6-dimethylaniline)]BF4 (12), and with I2 to liberate quantitatively the DMA-iminium salt 14. On the other hand, the chiral amine complex (R)-19 formed the optically inactive eta6-bound compound [Ir(cod)(eta6-rac-2 a)]BF4 (rac-18) upon dissolution in THF at room temperature, presumably via intramolecular C-H activation. This racemization was found to be a two-step event with k'1=9.0 x 10(-4) s(-1) and k2=2.89 x 10(-5) s(-1), featuring an optically active intermediate prior to sp3 C-H activation. Compounds 11, 12, rac-18, and (R)-19 were structurally characterized by single-crystal X-ray analyses.     

Structure and bonding in solution of dioxouranium(VI) oxalate complexes: isomers and intramolecular ligand exchange

Structural isomers of [UO(2)(oxalate)(3)](4-), [UO(2)(oxalate)F(3)](3-), [UO(2)(oxalate)(2)F](3-), and [UO(2)(oxalate)(2)(H(2)O)](2-) have been studied by using EXAFS and quantum chemical ab initio methods. Theoretical structures and their relative energies were determined in the gas phase and in water using the CPCM model. The most stable isomers according to the quantum chemical calculations have geometries consistent with the EXAFS data, and the difference between measured and calculated bond distances is generally less than 0.05 A. The complex [UO(2)(oxalate)(3)](4-) contains two oxalate ligands forming five-membered chelate rings, while the third is bonded end-on to a single carboxylate oxygen. The most stable isomer of the other two complexes also contains the same type of chelate-bonded oxalate ligands. The activation energy for ring opening in [UO(2)(oxalate)F(3)](3-), deltaU++ = 63 kJ/mol, is in fair agreement with the experimental activation enthalpy, deltaH++ = 45 +/- 5 kJ/mol, for different [UO(2)(picolinate)F(3)](2-) complexes, indicating similar ring-opening mechanisms. No direct experimental information is available on intramolecular exchange in [UO(3)(oxalate)(3)](4-). The theoretical results indicate that it takes place via the tris-chelated intermediate with an activation energy of deltaU++ = 38 kJ/mol; the other pathways involve multiple steps and have much higher activation energies. The geometries and energies of dioxouranium(VI) complexes in the gas phase and solvent models differ slightly, with differences in bond distance and energy of typically less than 0.06 A and 10 kJ/mol, respectively. However, there might be a significant difference in the distance between uranium and the leaving/entering group in the transition state, resulting in a systematic error when the gas-phase geometry is used to estimate the activation energy in solution. This systematic error is about 10 kJ/mol and tends to cancel when comparing different pathways.     

Kinetics and Thermodynamics of Reactions Involving Criegee Intermediates: An Assessment of Density Functional Theory and Ab Initio Methods Through Comparison with CCSDT(Q)/CBS Data

Reactions involving Criegee intermediates (CIs, R₁R₂COO) are important in atmospheric ozonolysis models. In recent years, density functional theory (DFT) and CCSD(T)-based ab initio methods are increasingly being used for modeling reaction profiles involving CIs. We obtain highly accurate CCSDT(Q)/CBS reaction energies and barrier heights for ring-closing reactions involving atmospherically important CIs (R₁/R₂ = H, Me, OH, OMe, F, CN, cyclopropene, ethylene, acetaldehyde, and acrolein). We use this benchmark data to evaluate the performance of DFT, double-hybrid DFT (DHDFT), and ab initio methods for the kinetics and thermodynamics of these reactions. We find that reaction energies are more challenging for approximate theoretical procedures than barrier heights. Overall, taking both reaction energies and barrier heights into account, only one of the 58 considered DFT methods (the meta-GGA MN12-L) attains near chemical accuracy, with root-mean-square deviations (RMSDs) of 3.5 (barrier heights) and 4.7 (reaction energies) kJ mol⁻¹. Therefore, MN12-L is recommended for investigations where CCSD(T)-based methods are not computationally feasible. For reaction barrier heights performance does not strictly follow Jacob's Ladder, for example, DHDFT methods do not perform better than conventional DFT methods. Of the ab initio methods, the cost-effective CCSD(T)/CBS(MP2) approach gives the best performance for both reaction energies and barrier heights, with RMSDs of 1.7 and 1.4 kJ mol⁻¹, respectively. All the considered Gaussian-n methods show good performance with RMSDs below the threshold of chemical accuracy for both reaction energies and barrier heights, where G4(MP2) shows the best overall performance with RMSDs of 2.9 and 1.5 kJ mol⁻¹, respectively. © 2019 Wiley Periodicals, Inc.     

(H_2O)_n(n=0～2)催化HS和HOCl反应机理的理论研究

在CCSD(T)/6-311+G(3df,2p)//M06-2X/6-311+G(3df,2p)水平上研究了(H_2O)n(n=0~2)催化HS和HOCl的反应机理.结果表明,HS与HOCl反应中HS夺取HOCl上的H原子形成产物H_2S和ClO.在无水催化时,该反应存在2种不同的路径(分别经过过渡态TS1和TS2,二者互为顺反结构),对应的能垒分别为100.28和100.91kJ/mol,到达产物(H_2S+ClO)需吸收18.99kJ/mol能量,反应不易发生;在单个水分子参与时,水分子可通过形成弱相互作用或者作为H原子转移桥梁影响反应机理,获得了4种水催化路径,能垒(间于53.97~92.39kJ/mol之间)均低于无水催化过程.同时发现,在反应到达产物前,水分子可以与产物形成中间体IM,IM相对能仅为0.46kJ/mol,有利于产物形成;有2个水分子参与反应时,找到了3条催化路径,最优反应路径过渡态TS7的能垒为45.05kJ/mol,低于无水催化过程,相比单个水分子最优路径能垒(53.97kJ/mol)并无显著降低.     

Modified Gaussian-2 level investigation of the identity ion-pair SN2 reactions of lithium halide and methyl halide with inversion and retention mechanisms

Identity ion-pair S(N)2 reactions LiX + CH(3)X --> XCH(3) + LiX (X = F, Cl, Br, and I) have been investigated in the gas phase and in solution at the level of the modified Gaussian-2 theory. Two possible reaction mechanisms, inversion and retention, are discussed. The reaction barriers relative to the complexes for the inversion mechanism [DeltaH(cent) ( not equal )(inv)] are found to be much higher than the corresponding values for the gas phase anionic S(N)2 reactions, decreasing in the following order: F (263.6 kJ mol(-1)) > Cl (203.3 kJ mol(-1)) > Br (174.7 kJ mol(-1)) > I (150.7 kJ mol(-1)). The barrier gaps between the two mechanisms [DeltaH(cent) ( not equal ) (ret) - DeltaH(cent) ( not equal ) (inv)] increase in the order F (-62.7 kJ mol(-1)) < Cl (4.4 kJ mol(-1)) < Br (24.9 kJ mol(-1)) < I (45.1 kJ mol(-1)). Thus, the retention mechanism is energetically favorable for fluorine and the inversion mechanism is favored for other halogens, in contrast to the anionic S(N)2 reactions at carbon where the inversion reaction channel is much more favorable for all of the halogens. The stabilization energies for the dipole-dipole complexes CH(3)X. LiX (DeltaH(comp)) are found to be similar for the entire set of systems with X = F, Cl, Br, and I, ranging from 53.4 kJ mol(-1) for I up to 58.9 kJ mol(-1) for F. The polarizable continuum model (PCM) has been used to evaluate the direct solvent effects on the energetics of the anionic and ion-pair S(N)2 reactions. The energetic profiles are found to be still double-well shaped for most of the ion-pair S(N)2 reactions in the solution, but the potential profile for reaction LiI + CH(3)I is predicted to be unimodal in the protic solvent. Good correlations between central barriers [DeltaH(cent) ( not equal ) (inv)] with the geometric looseness of the inversion transition state %C-X( not equal ), the dissociation energies of the C-X bond (D(C-X)) and Li-X bond (D(Li-X)) are observed, respectively.     

Enthalpies of formation, bond dissociation energies and reaction paths for the decomposition of model biofuels: ethyl propanoate and methyl butanoate

The complete basis set method CBS-QB3 has been used to study the thermochemistry and kinetics of the esters ethyl propanoate (EP) and methyl butanoate (MB) to evaluate initiation reactions and intermediate products from unimolecular decomposition reactions. Using isodesmic and isogeitonic equations and atomization energies, we have estimated chemically accurate enthalpies of formation and bond dissociation energies for the esters and species derived from them. In addition it is shown that controversial literature values may be resolved by adopting, for the acetate radical, CH3C(O)O(.-), DeltaH(o)(f)298.15K) = -197.8 kJ mol(-1) and for the trans-hydrocarboxyl radical, C(.-)(O)OH, -181.6 +/- 2.9 kJ mol(-1). For EP, the lowest energy decomposition path encounters an energy barrier of approximately 210 kJ mol(-1) (approximately 50 kcal mol(-1)), which proceeds through a six-membered ring transition state (retro-ene reaction) via transfer of the primary methyl H atom from the ethyl group to the carbonyl oxygen, while cleaving the carbon-ether oxygen to form ethene and propanoic acid. On the other hand, the lowest energy path for MB has a barrier of approximately 285 kJ mol(-1), producing ethene. Other routes leading to the formation of aldehydes, alcohols, ketene, and propene are also discussed. Most of these intramolecular hydrogen transfers have energy barriers lower than that needed for homolytic bond fission (the lowest of which is 353 kJ mol(-1) for the C(alpha)-C(beta) bond in MB). Propene formation is a much higher energy demanding process, 402 kJ mol(-1), and it should be competitive with some C-C, C-O, and C-H bond cleavage processes.     

Activation barriers for addition of methyl radicals to oxygen-stabilized carbocations

Activation barriers (DeltaHMe(double dagger)) for adding methyl radicals to ions of the general formula CH3CR=OCH3+ have been measured by looking at the threshold energies for the reverse reaction, dissociative photoionization of ethers of the general formula RC(CH3)2OCH3. Dissociation by loss of a methyl radical has more favorable thermochemistry than loss of R*, yet the onset of R* loss occurs at lower energies than loss of CH3*. In other words, the more endothermic dissociation exhibits a lower appearance energy. Contrathermodynamic ordering of appearance energies is observed for R = Et, nPr, iPr, tBu, and neopentyl. The sum of the appearance energy difference, DeltaAE, and the thermochemical difference (DeltaDeltaH, calculated using G3 theory) gives a lower bound for the barrier for adding methyl radical to CH3CR=OCH3+. More specifically, the difference between that activation barrier and the one for adding R* to (CH3)2C=OCH3+, DeltaHMe(double dagger)-DeltaHR(double dagger), equals DeltaAE + DeltaDeltaH and has values in the range 20-24 kJ mol(-1) for the homologous series investigated. There is no systematic trend with the steric bulk of R, and available evidence suggests that DeltaHR(double dagger) does not have a value >5 kJ mol(-1). The difference in barrier heights, DeltaHMe(double dagger)-DeltaHiPr(double dagger) for CH3* plus iPrC(CH3)=OX+ vs iPr* + (CH3)2C=OX+, has the same value, regardless of whether X = H or CH3. Mixing of higher energy electronic configurations provides a qualitative theoretical explanation for some (but not all) observed trends in barrier heights.     

CS2和O(3P)及N(4S)原子化学反应动力学

含硫化合物的氧化是大气化学中特别感兴趣的领域之一。研究较多的含硫化合物有SO_2、OCS和CS_2。其中O(~3P)+CS_2化学反应动力学研究在文献中已见过报导。通常认为O(~3P)和CS_2反应有三个反应通道:     

含有羧基配体的蝎型钒氧配合物的合成、结构及其热分解动力学

设计合成了两种新型的以聚吡唑硼酸盐、氨基酸为配体的钒氧配合物VO[phCH2CH(NH2)COO][HB(pz)3](1)和VO(3,5-Me2pz)[HB(3,5-Me2pz)3](CH3COO)(2). 通过元素分析、红外光谱对配合物进行了表征, 并利用单晶X射线衍射技术解析了它们的结构. 非等温热分解动力学研究表明, 配合物1和2的热分解反应都是分两步进行的. 通过计算, 配合物1热分解的第一步反应的可能机理为成核与生长(n=1/4); 第二步反应的可能机理为化学反应. 其非等温动力学方程分别为, dα/dT=(A/β)e-E/RT(1/4)(1-α)[-ln(1-α)]-3 和dα/dT=(A/β)e-E/RT(1-α)2. 分解反应的表观活化能分别是223.52 和331.94 kJ·mol-1; 指前因子ln(A/s-1)分别是49.67 和57.50. 配合物2 热分解的第一步反应的可能机理为化学反应; 第二步反应的可能机理为成核与生长(n=1/2). 其非等温动力学方程分别为, dα/dT=(A/β)e-E/RT(1-α)2, 和dα/dT=(A/β)e-E/RT(1/2)(1-α)[-ln(1-α)]-1. 分解反应的表观活化能分别是300.56 和444.72 kJ·mol-1; 指前因子ln(A/s-1)分别是75.53 和92.50.     

Formic acid catalyzed hydrolysis of SO3 in the gas phase: a barrierless mechanism for sulfuric acid production of potential atmospheric importance

Computational studies at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels are performed to explore the changes in reaction barrier height for the gas phase hydrolysis of SO(3) to form H(2)SO(4) in the presence of a single formic acid (FA) molecule. For comparison, we have also performed calculations for the reference reaction involving water assisted hydrolysis of SO(3) at the same level. Our results show that the FA assisted hydrolysis of SO(3) to form H(2)SO(4) is effectively a barrierless process. The barrier heights for the isomerization of the SO(3)···H(2)O···FA prereactive collision complex, which is the rate limiting step in the FA assisted hydrolysis, are found to be respectively 0.59 and 0.08 kcal/mol at the B3LYP/6-311++G(3df,3pd) and MP2/6-311++G(3df,3pd) levels. This is substantially lower than the ~7 kcal/mol barrier for the corresponding step in the hydrolysis of SO(3) by two water molecules--which is currently the accepted mechanism for atmospheric sulfuric acid production. Simple kinetic analysis of the relative rates suggests that the reduction in barrier height facilitated by FA, combined with the greater stability of the prereactive SO(3)···H(2)O···FA collision complex compared to SO(3)···H(2)O···H(2)O and the rather plentiful atmospheric abundance of FA, makes the formic acid mediated hydrolysis reaction a potentially important pathway for atmospheric sulfuric acid production.     

Photochemical reactions of (eta 5-cyclopentadienyl)bis(t-butylacrylate) rhodium with silanes: dynamics of isomer interconversion via Rh(eta 2-silane) species

The compound CpRh(C(2)H(3)CO(2)(t)Bu)(2) 1 has been synthesised as a mixture of two pairs of interconverting isomers which differ in the relative orientations of the alkene substituents. The four isomers have been fully characterised by NMR spectroscopy. When complex 1 is photolysed in the presence of a silane, HSiR(2)R'R(2)R'= Et(3), Me(3), HEt(2), (OMe)(3) and Me(2)Cl] the corresponding Si-H oxidative addition products CpRh(SiR(2)R')(H)(C(2)H(3)CO(2)(t)Bu) and CpRh(H)(2)(SiR(2)R')(2) are formed. The Rh(III) complexes CpRh(SiR(2)R')(H)(C(2)H(3)CO(2)(t)Bu) exist in two isomeric forms of comparable energy which interconvert in an intramolecular process that does not involve a reversible [1,3] hydride or [1,3] silyl migration. The hydride (1)H NMR resonances for these species consequently broaden before coalescing into a single peak. For R(2)R'= Et(3), the activation parameters for interchange from the major to minor isomer were Delta H++= 60.2 +/- 2 kJ mol(-1) and Delta S++= 8 +/- 9 J mol(-1) K(-1), while for R(2)R'= Me(3) and Et(2)H, Delta H++= 61.5 +/- 1 kJ mol(-1), Delta S++= 6 +/- 5 J mol(-1) K(-1), and Delta H++= 61.8 +/- 3 kJ mol(-1), Delta S++= 12 +/- 9 J mol(-1) K(-1) respectively for conversion from the major isomer to the minor. For these complexes an eta(2)-Rh-H-Si transition state or intermediate is consistent with the evidence. When R(2)R'=(OMe)(3) and Me(2)Cl the change in appearance of the hydride resonances is more complex, with the activation parameters for interchange from the major to minor isomer for the former species being Delta H++= 78.3 +/- 2 kJ mol(-1) and Delta S++= 30 +/- 7 J mol(-1) K(-1) while for Me(2)Cl the barrier proved too high to measure before decomposition occurred. The complex spectral changes could be simulated when a discrete eta(2)-Rh-H-Si intermediate was involved in the isomer interconversion process and hence silane rotation in all these systems is proposed to involve two isomers of CpRh(eta(2)-HSiR(2)R')(C(2)H(3)CO(2)(t)Bu).     

Highly accurate first-principles benchmark data sets for the parametrization and validation of density functional and other approximate methods. Derivation of a robust, generally applicable, double-hybrid functional for thermochemistry and thermochemical kinetics

We present a number of near-exact, nonrelativistic, Born-Oppenheimer reference data sets for the parametrization of more approximate methods (such as DFT functionals). The data were obtained by means of the W4 ab initio computational thermochemistry protocol, which has a 95% confidence interval well below 1 kJ/mol. Our data sets include W4-08, which are total atomization energies of over 100 small molecules that cover varying degrees of nondynamical correlations, and DBH24-W4, which are W4 theory values for Truhlar's set of 24 representative barrier heights. The usual procedure of comparing calculated DFT values with experimental atomization energies is hampered by comparatively large experimental uncertainties in many experimental values and compounds errors due to deficiencies in the DFT functional with those resulting from neglect of relativity and finite nuclear mass. Comparison with accurate, explicitly nonrelativistic, ab initio data avoids these issues. We then proceed to explore the performance of B2x-PLYP-type double hybrid functionals for atomization energies and barrier heights. We find that the optimum hybrids for hydrogen-transfer reactions, heavy-atoms transfers, nucleophilic substitutions, and unimolecular and recombination reactions are quite different from one another: out of these subsets, the heavy-atom transfer reactions are by far the most sensitive to the percentages of Hartree-Fock-type exchange y and MP2-type correlation x in an (x, y) double hybrid. The (42,72) hybrid B2K-PLYP, as reported in a preliminary communication, represents the best compromise between thermochemistry and hydrogen-transfer barriers, while also yielding excellent performance for nucleophilic substitutions. By optimizing for best overall performance on both thermochemistry and the DBH24-W4 data set, however, we find a new (36,65) hybrid which we term B2GP-PLYP. At a slight expense in performance for hydrogen-transfer barrier heights and nucleophilic substitutions, we obtain substantially better performance for the other reaction types. Although both B2K-PLYP and B2GP-PLYP are capable of 2 kcal/mol quality thermochemistry, B2GP-PLYP appears to be the more robust toward nondynamical correlation and strongly polar character. We additionally find that double-hybrid functionals display excellent performance for such problems as hydrogen bonding, prototype late transition metal reactions, pericyclic reactions, prototype cumulene-polyacetylene system, and weak interactions.     

石墨烯中的单空位缺陷对铂原子催化解离氧分子的影响

采用密度泛函理论中的B3LYP方法研究了石墨烯中的单空位缺陷对铂原子(Pt)催化解离O_2分子的影响.计算发现O_2分子首先通过[2+1]或[2+2]环加成作用吸附在以单空位缺陷石墨烯为载体的Pt上(Pt-SV),并以不同的路径进行解离,吸附能分别为-158.23和-152.45kJ/mol.由于石墨烯片上单空位缺陷的存在,O_2分子更容易吸附在单空位缺陷处的Pt上,并且O_2在Pt-SV上解离的能垒(130.25kJ/mol)也明显比在Pt-pristine上解离的能垒低(76.23kJ/mol).因此石墨烯上单空位缺陷的存在提高增加了Pt的催化能力.