转动取向影响因素的多元非线性分析

在不同拓扑特征的London-Eying-Polanyi-Sato(LEPS)型势能面上，应用经典轨线法计算了不同动态条件下的58个样本，应用多元非线性分析方法研究了势能面及动态条件对转动取向的影响.结果表明, 对于A＋BC→AB＋C类反应, 反应物分子的质量比、势能面出口处的宽度和相对平动能对产物分子的转动取向有重要作用，且发现不同影响因素之间的乘积对产物的转动取向具有更重要的作用.     

长寿命轨线形成的模式识别研究

设计了6个参数来定量表征共线势能面的几何特征,并寻48个具有不同几何特征的LEPS型势能面,在以7个参数表征的不同的动态条件下应用Venus94程序进行准经典轨线(QCT)计算。随后,应用模式识别研究了势能面的特征及动态因素对长寿命轨线形成的影响,初步证实了若干知的并发现了若干新的影响长寿命轨线形成的主要影响因素并据此建了长寿命轨线形成的定性预测模型。     

碰撞能量对反应Sr＋HF转动取向的影响

在推广LEPS势能面上,用经典轨线方法,研究了反应碰撞能量对反应Sr＋HF的转动取向的影响.计算结果与产物轨道角动量模型进行比较.计算结果表明,随着碰撞能量的增加,产物转动取向越强烈.     

电场中B2分子特性研究

利用密度泛函理论(DFT)B3LYP/6-3l1＋G(2d)方法研究在不同方向电场(0～＋0.02 a.u.)作用下的B2分子的基态键长、总能量、偶极矩、最高占据轨道(HOMO)能量、最低空轨道(LUMO)能量、能隙及势能曲线的变化规律. 结果表明: 在一定外加电场范围内, 随电场强度的增大, 分子键长变大; 总能量降低; 偶极矩增大; HOMO能级、LUMO能级均降低; 能隙依赖外电场方向, 平行分子轴(Z)方向电场使能隙递减, 垂直分子轴(X)方向电场使能隙递增; 分子势能降低, 平行分子轴线(Z)方向电场对分子势能的影响随着核间距的增大而增大, 原有的“势能平台”遭到破坏.     

氧在C－Ni（100）表面反应动力学的研究──Ⅱ．动力学计算结果

根据Ｏ２＋Ｃ－Ｎｉ（１００）表面反应的推广ＬＥＰＳ势能面，并用ＱＣＴ方法研究了该体系的反应动力学行为．结果给出了分子化学吸附和解离原子化学吸附及ＣＯ脱附反应的分子动态特征，以及态－态过程的吸附几率．考察了各种能量形式对吸附、脱附几率的影响和各类吸附、脱附几率的变化规律．     

基态Li2O＋分子的结构和势能面研究

以6-311＋＋G(d)为基函数, 采用CASSCF方法优化出Li2O＋分子的稳定构型为线形Li-O-Li (C∞V), 电子组态为2∏, 并对平衡核间距、离解能和基态简正频率进行了计算. 根据原子分子反应静力学原理, 导出了Li2O＋分子的合理的离解极限. 并运用多体展式理论方法首次导出了基态Li2O＋分子的分析势能函数, 绘出了势能等值图, 其势能等值图准确地再现了Li2O＋分子的平衡结构特征.     

PuC和PuC2的分子结构与势能函数

采用密度泛函B3LYP方法和相对论有效原子实理论模型优化出PuC和PuC2分子稳定构型,其电子状态分别为X5Σ－和X5A2.PuC2分子为C2v构型,其∠CPuC=147.67°,平衡核间距Re=0.22819 nm, 离解能De=5.543 eV, 并计算出谐振动频率:ν1=61.736 cm－1、ν2=229.894 cm－1、ν3=305.582 cm－1.在此基础上,运用多体项展式理论方法,导出了基态PuC2分子的分析势能函数,该势能面准确地再现了PuC2分子的稳定结构,并根据势能面等值图讨论了PuC＋C反应和Pu＋C2反应的势能面静态特征.     

~1CH_2 N_2O反应的势能面

利用密度泛函理论（B3LYP）计算了~1CH_2＋N_2O反应的反应物、中间体、过渡态及产物的几何构型.进而用从头算方法（QCISD（T））计算了单点能量.由此描绘了反应的势能面,确定了反应的最终产物通道为N_2＋H_2CO和NO＋HCN＋H.后者比前者有更大的分支比.N_2、H_2CO、NO、HCN的存在有待于实验检测.作者认为,反应在室温下是加成-消除机理,而在高温下可以通过直接取代的机理获得N_2＋H_2CO.     

1CH2＋N2O反应的势能面

利用密度泛函理论（B3LYP）计算了1CH2＋N2O反应的反应物、中间体、过渡态及产物 的几何构型.进而用从头算方法（QCISD（T））计算了单点能量.由此描绘了反应的势能面, 确定了反应的最终产物通道为N2＋H2CO和NO＋HCN＋H.后者比前者有更大的分支比.N2、H2CO 、NO、HCN的存在有待于实验检测.作者认为,反应在室温下是加成-消除机理,而在高温下 可以通过直接取代的机理获得N2＋H2CO.     

离子速度成像方法研究C8H17Br分子的光解动力学

用离子速度成像方法, 研究了长链C8H17Br分子在234 nm激光下的光解过程. 通过2+1共振增强多光子电离探测了两种光解产物Br*(2P1/2)和Br(2P3/2), 得到了它们的相对量子产率. 从光解产物Br*(2P1/2)和Br(2P3/2)的速度图像得到了能量和角度分布. 并根据相对量子产率和角度分布, 计算了不同解离通道的比例. 实验发现C8H17Br分子解离过程中大部分能量都转化为内能, 该能量分配可以较好地用软反冲模型来解释, 并分析了这种能量分配跟烷基大小的关系.     

Application of the modified Shepard interpolation method to the determination of the potential energy surface for a molecule-surface reaction: H2 + Pt(111)

We have used a modified Shepard (MS) interpolation method, initially developed for gas phase reactions, to build a potential energy surface (PES) for studying the dissociative chemisorption of H2 on Pt(111). The aim was to study the efficiency and the accuracy of this interpolation method for an activated multidimensional molecule-surface reactive problem. The strategy used is based on previous applications of the MS method to gas phase reactions, but modified to take into account special features of molecule-surface reactions, like the presence of many similar reaction pathways which vary only slightly with surface site. The efficiency of the interpolation method was tested by using an already existing PES to provide the input data required for the construction of the new PES. The construction of the new PES required half as many ab initio data points as the construction of the old PES, and the comparison of the two PESs shows that the method is able to reproduce with good accuracy the most important features of the H2 + Pt(111) interaction potential. Finally, accuracy tests were done by comparing the results of dynamics simulations using the two different PESs. The good agreement obtained for reaction probabilities and probabilities for rotationally and diffractionally inelastic scattering shows clearly that the MS interpolation method can be used efficiently to yield accurate PESs for activated molecule-surface reactions.     

A global potential energy surface and time‐dependent quantum wave packet calculation of Au + H2 reaction

A global potential energy surface (PES) corresponding to the ground state of AuH₂ system has been constructed based on 22 853 ab initio energies calculated by the multireference configuration interaction method with a Davidson correction. The neural network method is used to fit the PES, and the root mean square error is only 1.87 meV. The topographical features of the novel global PES are compared with previous PES which is constructed by Zanchet et al. (Zanchet PES). The global minimum energy reaction paths on the two PESs both have a well and a barrier. Relative to the Au + H₂ reactants, the energy of well is 0.316 eV on the new PES, which is 0.421 eV deeper than Zanchet PES. The calculation of Au(²S) + H₂(X¹Σ_g⁺) → AuH(X¹Σ⁺) + H(²S) dynamical reaction is carried out on new PES, by the time‐dependent quantum wave packet method (TDWP) with second order split operator. The reaction probabilities, integral cross‐sections (ICSs) and differential cross‐sections are obtained from the dynamics calculation. The threshold in the reaction is about 1.46 eV, which is 0.07 eV smaller than Zanchet PES due to the different endothermic energies on the two PESs. At low collision energy (<2.3 eV), the total ICS is larger than the result obtained on Zanchet PES, which can be attributed to the difference of the wells and endothermic energies.     

Quantum dynamics of the dissociation of H2 on Cu(100): dependence of the site-reactivity on initial rovibrational state

We perform six-dimensional (6D) quantum wavepacket calculations for H2 dissociatively adsorbing on Cu(100) from a variety of rovibrational initial states. The calculations are performed on a new potential energy surface (PES), the construction of which is also detailed. Reaction probabilities are in good agreement with experimental findings. Using a new flux analysis method, we calculate the reaction probability density as a function of surface site and collision energy, for a variety of initial states. This approach is used to study the effects of rotation and vibration on reaction at specific surface sites. The results are explained in terms of characteristics of the PES and intrinsically dynamic effects. An important observation is that, even at low collision energies, reaction does not necessarily proceed predominantly in the region of the minimum potential barrier, but can occur almost exclusively at a site with a higher barrier. This suggests that experimental control of initial conditions could be used to selectively induce reaction at particular surface sites. Our predictions for site-reactivity could be tested using contemporary experimental methods: The calculations predict that, for reacting molecules, there will be a dependence of the quadrupole alignment of j on the incident vibrational state, v. This is a direct result of PES topography in the vicinity of the preferred reaction sites of v = 0 and v = 1 molecules. Invoking detailed balance, evidence for this difference in preferred reaction site of v = 0 and 1 molecules could be obtained through associative desorption experiments.     

Greedy 3-Point Search (G3PS)—A Novel Algorithm for Pharmacophore Alignment

Chemical features of small molecules can be abstracted to 3D pharmacophore models, which are easy to generate, interpret, and adapt by medicinal chemists. Three-dimensional pharmacophores can be used to efficiently match and align molecules according to their chemical feature pattern, which facilitates the virtual screening of even large compound databases. Existing alignment methods, used in computational drug discovery and bio-activity prediction, are often not suitable for finding matches between pharmacophores accurately as they purely aim to minimize RMSD or maximize volume overlap, when the actual goal is to match as many features as possible within the positional tolerances of the pharmacophore features. As a consequence, the obtained alignment results are often suboptimal in terms of the number of geometrically matched feature pairs, which increases the false-negative rate, thus negatively affecting the outcome of virtual screening experiments. We addressed this issue by introducing a new alignment algorithm, Greedy 3-Point Search (G3PS), which aims at finding optimal alignments by using a matching-feature-pair maximizing search strategy while at the same time being faster than competing methods.     

NO3-+Cl2→ClONO2+Cl-反应势能面和势能阱

采用密度泛函理论B3LYP方法和6-311+G(d)基组, 计算构建离子-分子气相反应NO3-+Cl2→ClONO2+Cl-的三维势能面. 三维反应势能面证明该反应没有过渡态和势能垒, 但是存在一个深达-55.0 kJ/mol的势能阱(以氯气分子和硝酸根离子相隔无穷远为参量). 在势能阱底部, 有个化合物(O2NOClCl)- 称为势阱化合物, 依赖于势能阱而稳定存在. 理论红外光谱预测低温红外光谱能检测该势阱化合物. 低温条件下, 该反应由热力学控制, 反应产物是势阱化合物(O2NOClCl)-. 当温度升高, 该反应由动力学控制, 势阱化合物(O2NOClCl)-不稳定, 发生分解反应, 重新生成NO3-和Cl2. 研究结果可用来解释低温时ClONO2与Cl-气相反应不能产生Cl2的原因.     

气相中Ir＋循环催化N2O与CO反应机理的理论研究

采用密度泛函UB3LYP方法和Stuttgart赝势基组, 计算研究了气相中循环催化N2O ＋CO →N2 ＋CO2 反应的微观机理. 通过对相关物种亲氧性的计算, 证明了Ir＋循环催化作用在热力学上是可行的. 不同自旋态反应势能面的计算结果表明, 循环催化的两步反应均为自旋禁阻反应, 各存在不同自旋态势能面的交叉, 并运用Yoshizawa的内禀坐标单点垂直激发计算的方法找出了势能面交叉点; 两步反应均为放热反应, 总放热量为358.9 kJ&#8226;mol－1.     

Electronic structure and chemical bonding in MO(n)- and MO(n) clusters (M = Mo, W; n = 3-5): a photoelectron spectroscopy and ab initio study

Photoelectron spectroscopy (PES) and ab initio calculations are combined to investigate the electronic structure of MO(n)- clusters (M = W, Mo; n = 3-5). Similar PES spectra were observed between the W and Mo species. A large energy gap between the first and second PES bands was observed for MO3- and correlated with a stable closed-shell MO3 neutral cluster. The electron binding energies of MO4- increase significantly relative to those of MO3-, and there is also an abrupt spectral pattern change between MO3- and MO4-. Both MO4- and MO5- give PES features with extremely high electron binding energies (>5.0 eV) due to oxygen-2p-based orbitals. The experimental results are compared with extensive density functional and ab initio [CCSD(T)] calculations, which were performed to elucidate the electronic and structural evolution for the tungsten oxide clusters. WO3 is found to be a closed-shell, nonplanar molecule with C3v symmetry. WO4 is shown to have a triplet ground state (3A2) with D2d symmetry, whereas WO5 is found to be an unusual charge-transfer complex, (O2-)WO3+. WO4 and WO5 are shown to possess W-O* and O2-* radical characters, respectively.     

A barrier-free molecular radical-molecule reaction: {^{3}C_{2} (a^{3}\Pi) {+} O_{2} (X^{3} \Sigma)}

The reaction of ³C₂ (a³Π) radical with O₂ (X³Σ) molecule has been studied theoretically using ab initio Quantum Chemistry method. Both singlet and triplet potential energy surfaces (PES) are calculated at the CCSD(T)/aug-cc-pVDZ//B3LYP/6-311+G(d) + ZPE and G3B3 levels of theory. On the singlet PES of the title reaction, it is shown that the most feasible pathway should be the O-atom of O₂ attacking the C-atom of the  ³C₂ molecule first to form the adduct 1 CCOO, followed by the O-shift to give intermediate 2 CC(OO), and then to the major products P1 (2CO). Alternatively, 1 can be directly dissociated to P1 via transition state TS1-P1. The other reaction pathways are less competitive due to thermodynamical or kinetic factors. On the other hand, the pathways on the triplet PES are less competitive than those on the singlet PES in low temperature range, whereas it is not the case in high temperature ranges. On the basis of the analysis of the kinetics of all pathways through which the reactions proceed, we expect that the competitive power of reaction pathways may vary with experimental conditions for the title reaction. The reaction heats of formation calculated are in good agreement with that obtained experimentally.