不饱和硼烯HCIN=B：重排反应的DFT研究

采用量子化学中的密度泛函方法，在B3LYP/6-311G**水平上研究了不饱和硼烯HCIN=B：的重排反应机理。结果表明，无论是氢原子H迁移还是氯原子Cl迁移都经过1个三元环过渡态，生成直线型产物，但氢原子H迁移是在面内进行，而氮原子Cl迁移是在面外进行，根据计算结果，详细研究了不饱和硼烯HClN=B:重排反应的热力学及动力学函数，在此基础上讨论了不饱和硼烯HClN=B:的存在寿命问题。     

分子CH_3NH=B:的结构及重排反应的理论研究

采用量子化学中的从头计算方法, 在MP2/6-31G(d,p)水平上研究了不饱和硼烯CH3NH=B:的结构及重排反应机理。结果表明, CH3NH=B:的单线态结构比三线态结构稳定, 该分子的基态是单线态。分子CH3NH=B:可以发生3种不同的重排反应。本文找到了这3种重排反应的过渡态, 并详细计算了不饱和硼烯CH3NH=B:重排反应的动力学函数, 据此讨论了不饱和硼烯CH3NH=B:的稳定性问题。     

分子CH3NH=B：结构及重排反应的理论研究

采用量子化学中的从头计算方法，在MP2／6—31G(d，p)水平上研究了不饱和硼烯CH3NH=B：的结构及重排反应机理。结果表明，CH3NH=B：的单线态结构比三线态结构稳定，该分子的基态是单线态。分子CH3NH=B:可以发生3种不同的重排反应。本文找到了这3种重排反应的过渡态，并详细计算了不饱和硼烯CH3NH=B：重排反应的动力学函数，据此讨论了不饱和硼烯CH3NH=B：的稳定性问题。     

不饱和类硅烯H2C=SiNaCl的密度泛函理论研究

采用密度泛函理论方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类硅烯H2C=SiNaCl的结构.结果表明,不饱和类硅烯H2C=SiNaCl共有4种平衡构型,其中非平面的p-配合物型构型能量最低,是不饱和类硅烯H2C=SiNaCl存在的主要构型.对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒.计算预言了最稳定构型的振动频率和红外强度.     

不饱和类硅烯H2C=SiNaF的DFT研究

用密度泛函理论方法, 在B3LYP/6-31+G(d, p)水平上研究了不饱和类硅烯H2C=SiNaF的结构. 结果表明, 不饱和类硅烯H2C=SiNaF共有四种平衡构型, 其中非平面的p-配合物型构型能量最低, 是不饱和类硅烯H2C=SiNaF存在的主要构型. 对平衡构型间异构化反应的过渡态进行了计算, 求得了转化势垒. 计算预言了最稳定构型的振动频率和红外强度.     

LiBH4-X（X=O,F和Cl）体系解氢性能的第一原理计算

采用基于密度泛函理论的第一原理方法,计算了LiBH4-X(X=O,F和Cl)体系的晶体与电子结构及解氢性能.生成热和H原子解离能的计算结果表明:O原子掺杂优先占据LiBH4间隙位,F置换氢原子位,而Cl则取代BH4单元;O,F和Cl掺杂的LiBH4体系结构稳定性发生变化,其中O提高体系解氢效果明显,而F和Cl掺杂受H原子区域环境的影响.态密度、Mulliken电子占据数和电子密度的分析结果表明:B—H之间较强的共价键是LiBH4结构稳定、解氢困难的电子结构根源,O,F和Cl对LiBH4解氢能力影响主要是掺杂改变了H的s态与B的sp态的杂化特性、以及BH4单元与Li的成键作用.     

不饱和类锗烯H2C=GeLiCl的DFT研究

采用密度泛函理论方法, 在B3LYP/6-311G(d, p)水平上研究了不饱和类锗烯H2C=GeLiCl的结构及异构化反应. 结果表明, 不饱和类锗烯H2C=GeLiCl有三种平衡构型, 其中非平面的p-配合物型构型能量最低, 是其存在的主要构型. 对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒. 计算预言了最稳定构型的振动频率和红外吸收强度.     

不饱和类锗烯H2C=GeNaBr的密度泛函理论研究

应用密度泛函理论DFT方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类锗烯H2C=GeNaBr的结构及异构化反应.结果表明,不饱和类锗烯H2C=GeNaBr有3种平衡构型,其中非平面的p-配合物型构型能量最低,是其存在的主要构型.并对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒,预言了最稳定构型的振动频率和红外强度.     

不饱和类硅烯H_2CSiNaCl的密度泛函理论研究

采用密度泛函理论方法,在B3LYP/6-311G(d,p)水平上研究了不饱和类硅烯H2CSiNaCl的结构.结果表明,不饱和类硅烯H2CSiNaCl共有4种平衡构型,其中非平面的p-配合物型构型能量最低,是不饱和类硅烯H2CSiNaCl存在的主要构型.对平衡构型间异构化反应的过渡态进行了计算,求得了转化势垒.计算预言了最稳定构型的振动频率和红外强度.     

HCHO+X(X=F、Cl、Br)的反应机理

采用CCSD/6-311++G(d,p)//B3LYP/6-311++G(d,p)方法研究了HCHO与卤素原子X(X=F、Cl、Br)的反应机理. 计算结果表明, 卤素原子X(X=F、Cl、Br)主要通过直接提取HCHO中的H原子生成HCO+HX(X=F、Cl、Br). 另外还可以生成稳定的中间体, 中间体再通过卤原子夺氢和氢原子直接解离两个反应通道分别生成HCO+HX(X=F、Cl、Br)和H+XCHO(X=F、Cl、Br). 其中卤原子夺氢通道为主反应通道, HCO和HX(X=F、Cl、Br)为主要的反应产物; 且三个反应的活化能均较低, 说明此类反应很容易进行, 计算结果与实验结果符合很好. 电子密度拓扑分析显示, 在HCHO+X反应通道(b)中出现了T型结构过渡态, 结构过渡态(STS)位于能量过渡态(ETS)之后. 并且按F、Cl、Br的顺序, 结构过渡态出现得越来越晚.     

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

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

采用密度泛函理论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),反应能垒最低,为主反应通道,甲锗硫醛和异氰氢酸为主产物。锗烯经由四元环中间体抽提硫的反应为竞争反应通道。     

SiH2自由基与HNCO反应机理的理论研究

采用密度泛函理论B3LYP方法研究了SiH₂自由基与HNCO的反应机理, 并在B3LYP/6-311＋＋G**水平上对反应物、中间体、过渡态进行了全几何参数优化, 通过频率分析和内禀反应坐标(IRC)确定了中间体和过渡态. 为了得到更精确的能量值, 又用QCISD(T)/6-311＋＋G**方法计算了在B3LYP/6-311＋＋G**水平优化后的各个驻点的相对能量. 计算结果表明SiH₂自由基与HNCO的反应有五条反应通道, 其中顺式反应通道SiH₂＋HNCO→IM3→ TS4→IM5→TS5→IM6→SiH₂NH＋CO反应能垒最低, 为主反应通道.     

Theoretical determinations of ionization potentials of dopamine

Adiabatic and vertical ionization potentials (IPs) of nine conformers of dopamine in the gas phase are determined using density functional theory (DFT) B3LYP, B3P86, B3PW91 methods and high level ab initio HF method with 6-311++G** basis set, respectively. And the nine stable cationic states have been found in the ionization process of dopamine. Vertical ionization potentials of nine conformers of dopamine are calculated using the older outer-valence Green’s function (OVGF) calculations at 6-311++G** basis set. Vibrational frequencies and infrared spectrum intensities of G1b and G1b⁺ at B3LYP/6-311++G** level are discussed.     

Ab initio studies of aspartic acid conformers in gas phase and in solution

Systematic and extensive conformational searches of aspartic acid in gas phase and in solution have been performed. For the gaseous aspartic acid, a total of 1296 trial canonical structures and 216 trial zwitterionic structures were generated by allowing for all combinations of internal single-bond rotamers. All the trial structures were optimized at the B3LYP/6-311G* level and then subjected to further optimization at the B3LYP/6-311++G** level. A total of 139 canonical conformers were found, but no stable zwitterionic structure was found. The rotational constants, dipole moments, zero-point vibrational energies, harmonic frequencies, and vertical ionization energies of the canonical conformers were determined. Single-point energies were also calculated at the MP2/6-311++G** and CCSD/6-311++G** levels. The equilibrium distributions of the gaseous conformers at various temperatures were calculated. The proton affinity and gas phase basicity were calculated and the results are in excellent agreement with the experiments. The conformations in the solution were studied with different solvation models. The 216 trial zwitterionic structures were first optimized at the B3LYP/6-311G* level using the Onsager self-consistent reaction field model (SCRF) and then optimized at the B3LYP/6-311++G** level using the conductorlike polarized continuum model (CPCM) SCRF theory. A total of 22 zwitterions conformers were found. The gaseous canonical conformers were combined with the CPCM model and optimized at the B3LYP/6-311++G** level. The solvated zwitterionic and canonical structures were further examined by the discrete/SCRF model with one and two water molecules. The incremental solvation of the canonical and zwitterionic structures with up to six water molecules in gas phase was systematically examined. The studies show that combining aspartic acid with at least six water molecules in the gas phase or two water molecules and a SCRF solution model is required to provide qualitatively correct results in the solution.     

Hydroxylation Reaction Mechanism for Nitrosodimethylamine by Oxygen Atom

The hydroxylation reaction mechanism of nitrosodimethylamine(NDMA) by oxygen atom was theoretically investigated at the B3LYP/6-31G** level. It has been found that the path of the oxydation of the C―H bond is easier than the path involving a Singlet/Triplet crossing. The study of the potential surface shows that both solvent effect at B3LYP/6-31G** level and different method at more credible MP2/6-311G** level in the gas phase have no effect on the hydroxylation reaction mechanism. The oxidation hydroxylation process of NDMA by O is exothermic reaction and easy to occur.     

OXO (X=Cl,Br)与X (2P3/2)反应机理的理论研究

用密度泛函B3LYP/6-311+G**和高级电子相关的组态相互作用QCISD(T)/6-311+G**方法研究了OXO与X (2P3/2)双自由基反应的微观机理.研究结果表明该反应存在两个反应通道,产物分别为XO和X2+O2.由于形成产物XO的活化势垒较低,因而是主要反应通道,这与实验观察到的结果是一致的.而形成X2+O2的通道从动力学上看是不利的.     

Mechanism and Kinetics for the Reaction of NCS and OH Radicals

The mechanism and dynamical properties for the reaction of NCS and OH radicals have been investigated theoretically. The minimum energy paths (MEP) of the reaction were calculated using the density functional theory(DFT) at the B3LYP/6-311 G^** level, and the energies along the MEP were further refined at the QCISD(T)/6-311 G^** level. As a result, the reaction mechanism of the title reaction involves three channels, producing HCS NO and HNC SO products, respectively. Path Ⅰ and path Ⅱ are competitive, with some advantages for path Ⅰ in kinetics. As for path Ⅲ, it looks difficult to react for its high energy barrier. Moreover, the rate constant have been calculated over the temperature range of 8190-2500K using canonical variational transition-state theory (CVT). It was found that the rate constants for both path Ⅰ and path Ⅱ are negatively dependent on temperature, which is similarwith the experimental results for reactions of NCS with NO and NO2, and the variational effect for the rate constant calculation olavs an important role in whole temperature range.     

The influence of protecting the hydroxyl group of β-oxy-α-diazo carbonyl compounds in the competition between Wolff rearrangement and [1,2]-hydrogen shift. Density functional theory study and topological analysis of the charge density

The influence of protecting the hydroxyl group of a β-oxy-α-diazo carbonyl compound on the competition between the Wolff rearrangement (WR) and the [1,2]-hydrogen shift (HS) was investigated theoretically. Stationary points on the potential-energy surface were located with the B3LYP density functional and the 6-31G** basis set. For the basic system geometry optimisations at B3LYP/6-311+G** were performed to validate the reliability of the B3LYP/6-31G** calculations. Single-point energy calculations were carried out at the B3LYP/6-311+G** level on the B3LYP/6-31G**-optimised geometries. Further insight into the processes was achieved with the aid of the theory of “atoms in molecules” of Bader. The calculated energy barriers qualitatively predicted the yields of HS and WR obtained experimentally. In order to rationalise the calculated energy barriers, it was necessary to take into account not only the electronegativity of the protective groups but also the alignment of the migrating groups with the depletion sites at the carbene centre. Further, when the hydroxyl group was not protected the existence of an intramolecular hydrogen bond played an important role in both HS and WR. Received: 30 December 1998 / Accepted: 7 May 1999 / Published online: 4 October 1999