G3(QCI)模型化学方法

建立了非微扰外推模式下的几种G3(QCI)方法：G3(QCI/fu1)、G3(QCI/fu2)、G3(QCI/ful)//B3PW91和G3(QCI/fu2)//B3PW91.其中,电子能量用QCISD(T,FC)/G3large直接计算,芯电子相关能分别在MP2/6-31G(d)和MP2/6-31G(d,p)级别上计算,对125个G2-1 test set 的计算结果表明,总体精度与G3和G3 //B3LYP相当;平均绝对偏差分别为4.370、4.389、4.061和4.022kJ mol-1,相应G3和G3//B3LYP分别为4.27和4.05kJ mol-1.文章提出的方法排除了G3中外推办法的不确定因素,且更适用于非平衡几何构型体系能量的定量计算。     

含芯电子相关能修正的G2理论——G2(ful)

A general method in considering the core electronic correlation energies has been proposed and introduced into the standard Gaussian-2 (G2)[7] theory by small post-Hartree-Fock calculations. In this paper an additional MP2(FC)/6-31G(d) calculation over the G2 procedures is employed and examined in modification in modification to the flaw of Frozen-Core (FC) approximation of G2 vai eq.:
ΔE(full)= E[MP2(full)/6-31G(d)]-E[MP2(FC)/6-31G(d)]
where the MP2(full)/6-31G(d) energy has been obtained in the molecular geometry optimizations. This energy, ΔE(full), is directly added into the total G2 energy of a molecule in facilitating the effect of core electronic correlations for each molecule in chemical reactions. It has been shown that the over-all average absolute deviation for the 125 reaction energies of the G2 test set (test set 1) is slightly reduced from 5.09 to 5.01 kJ, mol(-1) while for the 55 D0 values, which have been used for the derivation of the A coefficient of the empirical High-Level...更多-Correction (HLC), it is also reduced from 4.99 [for both G2 and G2(COMPLETE)[8]]to 4.77 kJ• mol(-1). In addition, larger errors (greater than ±8.4 kJ•mol(-1) for the D0 energies are improved, especially for the largest error of the D0 of SO2 This error is reduced from 21.3 to 15.4 kJ. mol(-1), in which the experimental geometry would further reduce it by 7.1kJ.mol(-1)[8]. Another improvement is the absolute value of the A coefficient in HLC being reduced from 4.81 for G2 to 4.34 milli-hartrees which is believed to be useful in isolating the relationship between the HLC and the FC approximation. Modifications to the original G2 from this work is denoted as G2(fu 1) and thus the G2 (fu 1) total energy for a molecule is
E[G2(fu 1)]= E[G2]+Δ E(full)h
with a new ΔE[HLC] =-0.19α- 4.34nβ milli-hartree.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

H CH_2CO反应机理的G2计算

分别在UQCISD/6-311G(d,p)和G2理论计算水平上,对CH2CO和H反应可能存在的四条反应通道进行了研究,详细分析了每个通道的反应机理;通过振动分析的虚频数和内禀反应坐标(IRC)计算,确认了反应涉及的每一个过渡态.通过反应位能剖面的比较,发现经过一个中间体生成CH3 CO的一条途径是主反应通道,该通道是个放热反应,总焓变为-146.07kJ·mol-1,速控步骤的位垒为55.09kJ·mol-1.理论计算结果较好地解释了实验观察到的主要产物和副产物并存的现象.     

过充电保护添加剂1,4-二甲氧基苯的反应机理

采用密度泛函理论(B3LYP/6-311+G(d,p))和MP2/6-311+G(d,p)方法,研究锂离子电池过充电保护添加剂1,4-二甲氧基苯(p-DMOB)的作用机理.计算结果表明,在过充时,p-DMOB优先于溶剂分子(乙基甲基碳酸酯、二甲基碳酸酯、碳酸乙酯)发生氧化反应.用B3LYP和MP2计算所得的p-DMOB理论氧化电位接近,分别为4.12和4.05V(vsLi/Li+).p-DMOB氧化时首先失去一个电子,生成p-DMOB+·正离子自由基,用B3LYP和MP2方法计算所得的相应能量变化分别为701.24和728.27kJ·mol-1.失去电子后苯环的共轭性受到破坏,随后p-DMOB+·苯环上的C―H键发生断裂,失去H+并形成p-DMOB·自由基.用B3LYP和MP2方法计算所得的相应能量变化分别为1349.78和1810.99kJ·mol-1.p-DMOB·自由基很不稳定,会在电极表面发生聚合反应形成聚合物膜,用B3LYP和MP2方法计算所得的相应能量变化分别为-553.37和-1331.20kJ·mol-1.     

HPS2异构体结构与相对稳定性

在MP2/6-311++G(d,p)和QCISD(t)/6-311++G(3df,2p)(单点)水平下计算得到了包括9个异构体和15个过渡态的HPS2势能面.其中,异构体trans-HSPS(E1)的能量最低, 其次是cis-HSPS(E2)和HPS(S)(C2_V,E3),能量分别比trans-HSPS(E1)高3.43和14.17 kJ/mol.计算结果表明,异构体E1,E2,E3和立体的三元环结构的异构体HP(S)S(C_S,E4), 及与E4共存的trans-HPSS(E5)和cis-HPSS(E6)具有一定的稳定性.在QCISD(t)/6-311++G(3df,2p)//MP2/6-311++G(d,p)并包含零点能水平下,PH和S₂反应生成的E6和E5分别越过65.75和71.73 kJ/mol的势垒就可以异构化为具有较高动力学稳定性的产物E4, 计算结果对实验认定的产物是cis-HPSS(E6)的结论进行了修正.     

CHF_3与O(~3P)反应机理的理论研究

采用从头算MP2/6-311G(d,p)方法研究了CHF3与O(3 P)的反应机理,优化了所有反应物、产物和过渡态的几何构型,并通过振动频率分析和内禀反应坐标(IRC)方法确证了过渡态的真实性.在QCISD(T)/6-311++G(d,p)水平上精确计算了各反应物种的单点能.结果表明,标题反应共存在4类反应,分别为抽提氢反应(R1)、抽提氟反应(R2),消氟化氢反应(R3)和消氢反应(R4),在QCISD(T)/6-311++G(d,p)//MP2/6-311G(d,p)水平上,R1,R2,R3和R4反应的能垒分别为70.7,378.7,294.7和307.2kJ·mol-1,抽提氢反应为主反应通道.     

CH3F与C2H3反应机理的理论研究

采用密度泛函B3LYP/6-311G(d,p)方法对CH3F与C2H3的反应体系进行了理论研究,获得了反应的势能面信息及可能的微观机理.在QCISD(T)/6-311++G(d,p)水平上精确计算了各反应物种的单点能.结果表明,除抽提氢反应外,标题反应还存在抽提氟(R1)、消氟化氢(R2)、消氢(R3)和自由基形成(R4)四类反应.在QCISD(T)/6-311++G(d,p)//B3LYP/6-311G(d,p)水平上,R1,R2,R3和R4反应的能垒分别是163.9,152.2,209.8和224.2kJ·mol-1,相应反应能为-56.6,-164.3,-2.7和-156.0kJ·mol-1,所有反应均放热,为热力学允许的反应.     

GeH4与HX(X＝F,Cl,Br)间二氢键复合物的理论研究

对GeH4与HX形成的二氢键复合物的结构特征及本质进行了探讨.在MP2/6-311 ++G(3 d,3p)水平优化、频率验证得到复合物的分子结构,通过分子的几何参数及电子密度拓扑分析,确认GeH4与卤化氢已形成了二氢键复合物.MP2/6-311 ++ G(3d,3p)水平下进行BSSE校正后的结合能为3.281到4.5...     

类锗烯H2C=GeLiCl与RH(R=F,OH, NH2)的插入反应

采用DFT B3LYP和QCISD方法研究了不饱和类锗烯H2C=GeLiCl与RH(R=F, OH, NH2)的插入反应. 在B3LYP/6-311+G(d,p)水平上优化了反应势能面上的驻点构型. 结果表明, H2C=GeLiCl与HF、H2O 或NH3发生插入反应的机理相同. QCISD/6-311++G(d,p)//B3LYP/6-311+G(d,p)计算的三个反应的势垒分别为173.53、194.48和209.05 kJ·mol-1, 反应热分别为60.18、72.93和75.34 kJ·mol-1. 相同条件下发生插入反应时, 反应活性顺序都是H—F>H—OH>H—NH2.     

Criegee中间体气相反应热力学的G2理论计算

Criegee intermediate is believed to play an important role in the atmospheric chemistry. Because of its short life and the difficulty in experimental study, we carried out ah initio calculations on the thermochemistry of the Criegee involving reactions in this study. Thermochemistry data of reaction enthalpies and Gibbs free energies for four different stable structures of the Criegee intermediates (singlet CH2OO ①1 A1 in C2v, triplet CH2OO ②3B1 in C2v, singlet CH2OO ③1A' in Cs and triplet CH2OO ④ in C1 symmetry) involved in some of the gas-phase reactions were calculated at the standard Gaussian-2 [G2(MP2) and G2] and a modified G2, G2(fu1)[10],levels of theory. Relative energies among those Criegees and formic acid were compared. Chemical reactions include the formation of Criegees, re-arrangement from Criegee to formic acid, dissociations (producing CH2(3B1)+O2, CH2(1A1)+O2, CO2+H2, CO2+2H, CO+H2O, OH+HCO) and the reactions between Criegee and NO/H2O. Standard equilibrium constants for some reactions were investigated and may be obtained for all of the rest reactions involved in this study by the standard Gibbs free energies. It is shown that the formation of Criegee ①-④ by ethylene and ozone, the re-arrangement from any Criegee to formic acid, the dissociation in producing CO2+O2and CO+H2O and the reactions between any Criegee and NO/H2O are all favourable thermodynamicaly. The dissociation in forming CO2+2H and OH+HCO is less favourable. While the dissociation in forming carbene (either in 3B1 or 1A1 state) is not allowed by ΔrGm? values. Standard enthalpies of formation at 298 K for the four Criegees were predicted at the G2(ful) level of theory. Each value is the average value from ten of the above reactions and they are -4.3, 74.8,98.9 and 244.6 kJ mol-1 at the G2(ful) level for Criegee ① to Criegee ④, respectively. In addition, tile standard enthalpy of formation at 298 K for HOCH2OOH is further predicted to be -315.6 kJ mol-1 at the G2(MP2) level.     

等电子-等自旋与非等旋反应的G2(MP2)和G2研究

我们在前文［1－2］中分别用MP2－4／6－31G**／／MP2／6－31G**及MP4／6－311G（2df，Zpd）／／MP2／6－31G**研究了一些双原子氢化物、卤化物、硫化物和氧化物的化学反应的烂变·这些化学反应按如下类型分为四组，即（1）反应物与生成物之间为等电子一等自旋关系，（2）价层等电一等旅，（3）等施和（4）非等旅·结果表明，MP4／6－3fiG（2才，ZPd）对于（1），（2）和（3）类反应，基本上与实验误差小于士15kJ·mo－‘而对非等旋反应仍有较大误差；MPZ－4／6－31G””只对（1）类反应较好．由于PoPle等人近几年来创立的Gaussi…     

分子几何构型优化方法的系统性比较

对《ＣＲＣ物理化学手册》第７７版中收集的第三周期以前的所有已知实验构型的无机分子的构型,以ＭＰ２、Ｂ３ＬＹＰ、Ｂ３ＰＷ９１级别上进行了构型优化的系统性比较,优化采用基组为６－３１Ｇ（ｄ,ｐ）,６－３１１Ｇ（ｄ,ｐ）,６－３１１Ｇ（２ｄ,ｐ）。对大多数分子另比较了ＱＣＩＳＤ（Ｔ）方法,对多原子分子比较了ＢＰＷ９１方法,对含氢双原子分子也使用ＱＣＩＳＤ方法。结果表明,键长的平均绝对偏差（单位：ｐｐｍ）     

Molecular Dipole Moment Computed with Ab Initio MKS Charges

1 INTRODUCTION Dipole moment (DM) can be measured by a variety of experimental methods[1] or computed with atomic charge distribution directly derived from molecular orbital calculations[2~10]. The average errors in DM computed with AM1, PM3 and MNDO methods for 125 organic compounds are 0.35 D, 0.45 D and 0.45 D, respectively[2]. Feller et al. reported that the average errors in the computed DM at the HF level are 0.65 D, 0.49 D and 0.30 D with the basis set STO-3G, 3-21G and…     

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

F-(H2O)n(n=1～4)体系的结构、偶极矩、极化率和一阶超极化率的从头计算

用MP2方法得到F^-(H₂O)_n(n=1～4)体系的优化几何结构.依据优化结构计算偶极矩μ₀、极化率α₀、各向异性极化率Δα和一阶超极化率β₀,研究了基组效应.讨论了溶剂水分子数、溶剂水分子排布及溶剂层数变化诸因素对体系超极化率等性质的影响.     

HBO2异构体的结构和相对稳定性

用ab initio方法在MP2/6 311＋＋G(d,p)水平下优化得到了HBO2体系的若干异构体和过渡态，并在QCISD(t)/6 311＋＋G(3df,2p)//MP2/6 311＋＋G(d,p)水平下进行了单点能量校正.对计算结果的分析表明，无论是在热力学还是在动力学上，具有链状结构的HOBO异构体(E1)是势能面上最稳定的结构，并对E1的电子结构进行了分析；另一具有C2v对称性的HBO(O)结构的异构体(E2)的能量比E1高381.72 kJ•mol－1，由于E2处于一个较深的势垒中，因此是比较稳定的，可以推断，在适合的实验中应该可以观察到异构体E2.     

Theoretical study on the hydrogen abstraction reaction of CH3CH2F (HFC-161) with Cl atom

A direct dynamics method is employed to study the hydrogen abstraction reaction of CH₃CH₂F+Cl. Three distinct transition states are located, one for -H abstraction and two for β-H abstraction. The potential energy surface (PES) information is obtained at the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p), CCSD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) and G2//MP2/6-311G(d,p) level. Based on the QCISD(T)/6-311+G(3df,2p)//MP2/6-311G(d,p) results, the rate constants of the three reaction channels are evaluated by using the canonical variational transition state theory (CVT) with small-curvature tunneling (SCT) contributions over the temperature range of 220–2800 K. The calculated results indicate that -H abstraction dominates the total reaction almost over the whole temperature range.     

Assessment of ab initio MP2 and density functionals for characterizing the potential energy profiles of the SN2 reactions at N center

The potential energy profiles of five selected bimolecular nucleophilic substitution (S_N2) reactions at nitrogen (N) center have been reinvestigated with the CCSD(T), G3[MP2,CCSD(T)], MP2, and some density functional methods. The basis sets of 6‐31+G(d,p) and 6‐311+G(3d,2p) are used for the MP2 and density functional calculations. Taking the relative energies at the CCSD(T)/CBS level of theory as benchmarks, we recommend the MP2, B97‐K, B2K‐PLYP, BMK, ωB97X‐D, M06‐2X, M05‐2X, CAM‐B3LYP, M08‐SO, and ωB97X methods to generally characterize the potential energy profiles for the S_N2 reactions at N center. Furthermore, these recommended methods with the relatively small 6‐31+G(d,p) basis set may also be used to perform direct classical trajectory simulations to uncover the dynamic behaviors of the S_N2 reactions at N center. © 2012 Wiley Periodicals, Inc.