NH自由基与臭氧反应机理的理论研究

用量子化学从头计算方法,在HF/6-31 ++ G**水平研究了臭氧与NH三线态活性自由基反应的微观机理,优化得到反应物、过渡态、中间体和产物的几何构型.用MP2/6-31++G**//HF/6-31 ++ G**方法计算能量,同时进行零点能校正.研究结果表明:NH三线态活性自由基与O3反应首先生成稳定中间体HNO3,然后中裂解生成HNO和O2.     

Al6(OH)18(H2O)6的结构及成键方位的从头算及密度泛函分析Ⅱ

用RHF法6-31G、6-31G*、6-31G**水平的量子化学从头算及RB3LYP/6-31G、B3LYP/6-31G*,B3LYP/6-31G**水平的密度泛函方法,采用SCRF=Dipole溶剂模型,计算了三水铝石有利生长基元最优结构Al6(OH)18(H2O)6的总能量、布居数、原子静电荷等.计算结果表明,桥联OH基团更易成键,Al6(OH)18(H2O)6较为有利的成键方位是桥联OH基团方位.     

3-羟基丙烯醛单键旋转异构光反应机理的从头算研究

应用从头算分子轨道法分别在RHF/6-31G**和UHF/6-31G**水平上对3腔┤┓肿拥幕?S0)和三重激发态(T1)单键旋转异构反应机理进行了研究,优化出反应物和产物在S0态和T1态的4种平衡态和过渡态的几何构型,通过振动分析得到的虚频和计算的内禀反应坐标对过渡态进行了确认,并得到了零点能,根据基态和激发态反应途径分析了光反应机理。计算结果表明,基态和激发态的3腔┤┓肿拥ゼ旃狗从ξ焕萁系停嘏湃菀追⑸永砺凵涎橹ち斯夥从κ笛橹屑钢忠旃固宓拇嬖凇?     

1,2-C4H6→2-C4H6异构化反应机理的理论研究

采用自洽场分子轨道UHF/6-31G**从头算法,研究了1,2-C4H6→2-C4H6异构化反应机理,优化了基态势能面上反应物、过渡态、中间体和产物的几何构型,并对各驻点能量进行了零点能校准.结果表明该反应经历一个1-甲基环丙烯生成产物比经两步氢迁移反应历程更易发生.     

3-羟基-3-甲基-2-丁酮和苯甲酰甲酸甲酯热分解反应的 理论研究

用从头算的方法在6-31G水平上研究了3-羟基-3-甲基-2-丁酮(1)和苯甲酰甲酸甲酯(2)热分解反应的机理。结果是:前一反应是经历五元环过渡态到达氢键中间体,它接着直接分解成乙醛的异构体和丙酮,最后乙醛的异构体异构化成乙醛;后一反应经历六元环过渡态形成中间体1(INT1),中间体1(INT1)直接分解成中间体2(INT2)和甲醛,中间体2(INT2)经过第二个过渡态分解成苯甲醛的异构体和一氧化碳,最后苯甲醛异构体异构化成苯甲醛。其中氢迁过程是反应的速控步骤。在MP~2/6-31G//HF/6-31G+ZPE水平上,对应于这两个反应速控步骤的活化位垒分别是251.42kJ/moL和247.94kJ/mol。采用传统过渡态理论计算了两反应的热反应速率常数,理论的计算结果与实验值吻合较好。     

单重态CCl2与O3反应机理的理论研究

用量子化学从头算方法,研究了单重态CCl_2与O_3反应的机理.在HF/6-31G(d)水平上用梯度解析技术全参数优化上述反应的反应物、中间体、过渡态和产物构型,MP2/6-31G(d)//HF/6-31G(d)方法计算能量.给出了有关化合物的结构数据.结果表明:CCl_2与O_3首先生成富能中间体CCl_2O_3,然后中间体裂解生成CCl_2O和O_2.该反应为强放热反应,放出的热量为516.88kJ·mol~(-1)[MP2/6-31G(d)//HF/6-31G(d)].通过内禀反应坐标(IRC)计算,获得了沿反应途径的势能剖面.     

硅炔氢加成生成甲硅烷基硅烯的量子化学研究

利用量子化学从头算和密度泛函理论(DFT), 研究了硅炔和氢气分子加成生成甲硅烷基硅烯的反应机理. 在B3LYP/6-311G**水平上, 全参数优化了反应通道上各驻点(反应物、过渡态和产物)的几何构型, 计算出了它们的振动频率和零点振动能(ZPVE), 并对它们进行了振动分析, 以确定过渡态的真实性. 各物质总能量由QCISD(T)/6-311G**// B3LYP/6-311G**给出, 并对能量进行了零点能校正. 计算表明, 硅炔与氢分子加成反应可生成稳定的甲硅烷基硅烯. 热力学与动力学计算表明, 反应过程是一个放热、熵减少而自发趋势和反应程度较大的反应.     

合成环氧乙烷新途径的从头算研究

在MP4/6-31G*//RHF/6-31G*理论水平上，对单态甲撑插入甲醛中碳氧双链生成环氧乙烷的反应(CH2_(~1A)＋CH_2O→C_2H_4O)进行了从头算研究。发现其反应历程由两步组成：1) 反应物沿反应坐标接近，体系能量单调下降，生成平而松散环状的分子复合物(MC)；2) 分子复合物沿反应坐标经由一过渡态(TS)重排为产物环氧乙烷，此步的势垒只有36.99 kJ·mol~(-1)。进而计算了该反应的热力学函数和动力学性质，并进行了讨论。当设法将甲醛引入产生活性中间体CH_2(~1A)的体系中，该反应有可能成为在室温下制备环氧乙烷的非催化途径。     

可乐定分子构象异构和互变异构的理论研究

采用HF/6-31G*方法, 对氨基型和亚胺型可乐定中性分子以及质子化的可乐定分子进行了构象分析. 之后采用B3LYP/6-31G**方法对处于势能面上的能量极小构象异构体进行全自由度几何优化和频率分析. Onsager反应场溶剂模型用于水相的计算. 结果表明, 在气相和水相中可乐定中性分子主要以亚胺型存在, 这同已有实验结果一致. 进一步, 寻找构象异构化过渡态和氨基型亚胺型互变异构化过渡态, 探讨质子化和水溶剂化效应对异构化过程的几何结构和能量的影响. 为了考察氯苯的共轭效应对可乐定互变异构体稳定性的影响, 在B3LYP/6-31G**水平上, 研究了2-氨基-2-咪唑啉的互变异构化反应机理.     

气态烯丙醇光异构化反应的理论研究

本文采用自洽场分子轨道UHF/6-31G从头计算法, 辅以能量梯度法研究气相烯丙醇光异构化反应的机理。全部优化了T1态势能面上反应物、过渡态、中间体和产物的几何构型。基于Fukui提出的内禀反应坐标理论(IRC)计算这一体系的反应途径, 并针对各驻点进行MP2/6-31G的相关能校正, 得到该反应在激发态进行并为一经历双自由基中间体的分步反应的结论。支持了实验工作者提出的机理。     

DFT法研究3-羟基丙烯醛的双键旋转异构反应机理

利用密度泛函理论(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了基态3-羟基丙烯醛分子在双键旋转异构反应过程中的平衡态以及过渡态的几何构型，分析了反应过程中键参数的变化，计算了该反应的内禀反应坐标(IRC)，发现在重排反应途径上存在一个四元环骨架的中间体.通过振动分析对平衡态和过渡态进行了确认，并得到了零点能.计算结果表明，基态3-羟基丙烯醛分子的双键旋转异构反应经过两步完成，第一步反应位垒稍高，第二步反应位垒较低，存在着发生重排反应的可能性.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

DFT法研究离子液中EMIM＋催化丁烯双键异构反应机理

利用密度泛函方法(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了离子液体中1-乙基-3-甲基咪唑阳离子(EMIM＋)催化丁烯双键异构反应过程中的反应物、产物以及过渡态的几何构型,分析了反应过程中键参数的变化.通过振动分析对平衡态和过渡态进行了验证,并得到了零点能.通过计算内禀反应坐标(IRC),确认了对应于过渡态的反应物和产物. 计算结果表明,EMIM＋催化丁烯双键异构可以基元反应的方式一步完成,1-丁烯异构化为2-丁烯的活化能约为192 kJ•mol－1, 逆反应活化能约为208 kJ•mol－1, 可在室温或高于室温条件下进行.     

乙烯醇钠的从头算研究(英文)

在ＲＨＦ／３－２１Ｇ和ＲＨＦ／６－３１＋Ｇ水平上优化了乙烯醇钠的的平衡构型及过渡态，比较了不同构型的稳定性。在ＲＨＦ／６－３１＋Ｇ水平上探讨了气相中由乙醛和氢氧化钠生成乙烯醇钠的机理。并在ＭＰ２（ｆｕｌ）／６－３１＋Ｇ／／ＭＰ２（ｆｕｌ）／６－３１＋Ｇ水平下计算了反应热     

HNNH3的从头计算研究

Ab-Initio studies applying the 3-21G, 6-31G and 6-31G** basis sets and also including the MP2 correction were carried out on H2NNH2, HNNH3 and the transition state molecule of the reaction H2NNH2=HNNH3. First, the geometry of the three molecules was optimized using the theoretical methods mentioned in the Hartree-Fock(HF) scheme. The energies of the molecules corresponding to RHF/6-31 G** geometries were subsequently calculated including electron correlation effects at the level of the second-order Moller-Plesset(MP2) perturbation theory. The vibrational frequencies, net charges and dipole moments were obtained from the theoretical calculations. The results of our ab-initio calculations indicate unambiguously that H2NNH2 is thermodynamically more stable than HNNH3. On the other hand, an isolated HNNH3 molecule once created would be rather stable since barriers for its unimolecular isomerization and decomposition are relatively high. Nevertheless, HNNH3 can be considered as an intermediate in chemical processes involving N2H4. This assumption has been supported by further calculation at higher level of theory.     

Mechanistic study of the deamination reaction of guanine: a computational study

The mechanism for the deamination of guanine with H(2)O, OH(-), H(2)O/OH(-) and for GuaH(+) with H(2)O has been investigated using ab initio calculations. Optimized geometries of the reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), B3LYP/6-31G(d), and B3LYP/6-31+G(d) levels of theory. Energies were also determined at G3MP2, G3MP2B3, G4MP2, and CBS-QB3 levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and an intermediate in the last step that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with OH(-) was found to have an activation energy barrier of 155 kJ mol(-1) compared to 187 kJ mol(-1) for the reaction with H(2)O and 243 kJ mol(-1) for GuaH(+) with H(2)O. The lowest overall activation energy, 144 kJ mol(-1) at the G3MP2 level, was obtained for the deamination of guanine with H(2)O/OH(-). Due to a lack of experimental results for guanine deamination, a comparison is made with those of cytosine, whose deamination reaction parallels that of guanine.     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

异硫氰酸与甲亚胺环加成反应机理的理论研究

对异硫氰酸与甲亚胺形成四元或六元环产物的环加成反应进行了理论研究,结果表明,一分子异硫氰酸与一分子甲亚胺形成四元环(1:1)产物P1的反应(1)为经过一个两性离子中间体的分步反应,其中第二步为速控步骤,其活化热垒为107.86kJ/mol.此外,反应(1)的中间体还可与另一甲亚胺或异硫氰酸分子继续反应形成两个不同的六元环(1:2或2:1)产物P2或P3;这两个反应均为协同反应,其活化势垒分别为15.88和21.82kJ/mol.这些结果与当异硫氰酸酯与亚胺发生环加成反应时,只有类似于P2和P3的两种六元环产物生成的实验事实一致。     

锗烯与异硫氰酸反应机理的量子化学研究

采用密度泛函理论B3LYP方法研究了GeH2自由基与HNCS的反应机理,并在B3LYP/6-311++G**水平上对反应物,中间体,过渡态进行了全几何参数优化,通过频率分析和IRC确定中间体和过渡态。为了得到更精确的能量值,用QCISD(T)/6-311++G**方法计算了各个驻点的单点能,计算结果表明单重态的锗烯与异硫氰酸的反应有抽提硫、插入N-H键、抽提亚氨基的路径,而经由三元环中间体的抽提硫反应GeH2+HNCS→IM3→TS2→IM4→TS3→IM5→GeH2S+HNC(P1),反应能垒最低,为主反应通道,甲锗硫醛和异氰氢酸为主产物。锗烯经由四元环中间体抽提硫的反应为竞争反应通道。     

乙烯醇锂的从头算研究

用限制的HF/3-21G和HF/6-31G^*优化乙烯醇锂的几种可能构型,比较了它们的稳定性.用限制的HF/3-1G,从乙醛开始,探讨了气相反应生成乙烯醇锂的机理,并在MP2水平上用6-31G^*基组计算了反应热.