丙烯在MoO3和WO3缺氧表面吸附的原子簇模型的从头算研究

本文选取原子簇模型, 用从头算UHF方法优化几何构型, 并进行二级MΦller-Plesset微扰(MP2)相关能计算, 研究了丙烯在MoO3和WO3表面的吸附态结构。优化得到的吸附前后的几何构型及计算的原子净电荷和吸附能都一致地表明,WO3对丙烯的吸附作用强于MoO3; 丙烯均用π成键电子向中心原子Mo或W配位而被吸附在MoO3的缺氧表面, 使丙烯中的碳-碳双键减弱致使容易断裂, 进而起到促进歧化反应的作用。从MP2计算的WO3对丙烯的吸附能0.934eV大于MoO3对丙烯的吸附能0.506eV及吸附过程均不必翻越能垒, 还可以定性理解用WO3作催化剂时的丙烯歧化反应的最佳温度(675K)高于用MoO3作催化剂(均以Al2O3为载体)的最佳反应温度(480K), 且WO3在675K时的催化活性高于MoO3在480K时的活性。     

含S-P-N键的磷烯正离子及其相关物的结构和性质的从头计算

本文对5-甲硫基-3-N-苯基-2,3-2(H)-1,3,4,2-噻二唑磷茂啉的磷烯正离子等三个相关分子的平衡几何构型进行了从头算解析能量梯度方法的全优化计算. 优化结果表明, 二配位磷烯正离子与相应的三配位磷母体分子的结构有本质的差异.形成二配位磷烯正离子后,磷所在的五员环形成共轭体系而使原来的单键键长变短, 原来的双键键长变长, 且使原来不共面的五员环共面. 在优化的平衡几何处进行单点CI计算的结果表明, 基态分子中磷原子上的正电荷的相对多少次序与实验测得的摩尔电导率及^3^1P NMR谱的化学位移的相对大小次序完全一致.     

二苯二硫和二苄二硫润滑性能的量子化学研究

The molecular geometries optimization and electronic structures of diphenyl disulfide (DPDS) and dibenzyl disulfide (DBDS) compounds were investigated by density functional theory (DFT) and ab initio method at the 6-31G basis set level. The active atoms and bonds of reaction were provided by frontier molecular orbital theory. The molecular orbital parameters of DPDS and DBDS compounds and iron atom cluster were calculated by using density functional theory. The interaction pattern between the organic disulfide compounds and iron atom cluster was discussed based on the approximate rule of orbital energy. Some parameters characterizing the action strength between the organic disulfide compounds and iron atom cluster, including the bonding strength, reactive strength and static action strength, were analyzed by using frontier electron density, super de-localizability, net atomic charge and the interaction energy of chemical adsorption as criteria. The results indicate that S-S chemical bond and C-S chemical bond of the compounds are inclined to be broken when DPDS and DBDS interact with the metal. The anti wear ability order of DPDS and DBDS compounds is DPDS>DBDS, and the extreme pressure ability order of DPDS and DBDS compounds is DBDS>DPDS, and the prediction results based on quantum chemistry calculations are in good accordance with the friction and wear test results.     

水在石墨(0001)面簇模型桥位上吸附的量子化学研究

用从头计算方法对水在石墨(0001)面桥位上的吸附进行了研究.用C₆H₈原子簇模拟石墨表面,在6-31G^*水平上计算了水在不同方向和位置上的吸附能量.研究表明:水在石墨面上的吸附很弱,属于物理吸附;在中性或带负电荷的石墨表面,当水分子中的氢原子靠近石墨面时,体系存在能量最小值,而在带正电荷的表面,当氧原子靠近石墨面时存在稳定的吸附点;不论表面带正电荷还是带负电荷,均对水分子的吸附起增强作用.     

H2O在Cu(100)表面吸附的从头算研究

用量子化学从头计算方法,以原子簇Cu₅为模拟表面,研究了水在Cu(100)面上不同吸附位的吸附情况,结果表明:水分子通过氧原子与表面成键,顶位是其最佳吸附位,吸附能约为70kJ/mol,平衡距离为0.213nm,氢原子远离表面.在氧原子不加极化函数时,水分子的二次轴垂直于表面时能量最低,但倾斜至50°所需能量仅在10kJ/mol以内.当考虑O原子d轨道的影响时,水分子倾斜时能量较低,得到了与实验相符的吸附构型.另外还研究了表面电荷对吸附体系的影响,结果表明:表面带正电荷时,水与表面间的相互作用增强,水上所转移电荷增多,Cu-O间平衡距离减小;表面带负电荷时,情况与之相反,且氢原子靠近表面时,势能曲线有最低点.     

Adsorption of atomic and molecular oxygen on the LaMnO3(001) surface: ab initio supercell calculations and thermodynamics

We present and discuss the results of ab initio DFT plane-wave supercell calculations of the atomic and molecular oxygen adsorption and diffusion on the LaMnO(3) (001) surface which serves as a model material for a cathode of solid oxide fuel cells. The dissociative adsorption of O(2) molecules from the gas phase is energetically favorable on surface Mn ions even on a defect-free surface. The surface migration energy for adsorbed O ions is found to be quite high, 2.0 eV. We predict that the adsorbed O atoms could penetrate the electrode first plane when much more mobile surface oxygen vacancies (migration energy of 0.69 eV) approach the O ions strongly bound to the surface Mn ions. The formation of the O vacancy near the O atom adsorbed atop surface Mn ion leads to an increase of the O-Mn binding energy by 0.74 eV whereas the drop of this adsorbed O atom into a vacancy possesses no energy barrier. Ab initio thermodynamics predicts that at typical SOFC operation temperatures (approximately 1200 K) the MnO(2) (001) surface with adsorbed O atoms is the most stable in a very wide range of oxygen gas pressures (above 10(-2) atm).     

Extended electron distributions applied to the molecular mechanics of some intermolecular interactions

Summary Extended electron distributions (XEDs) have been added to the molecular mechanics Coulombic term and applied to a selection of intermolecular interactions. The results from this approach have been compared with the commonly used atom-centred charges and more rigorous AM1-derived natural atom orbital point densities. The use of XEDs generally improves the simulation of experimental and ab initio results over the other two charge allocations and corrects geometries in those cases for which the others yield wrong results.     

铁原子与氮分子间的相互作用——单端位构型

用杂化密度泛函B3LYP方法在6-311+G(d)基组水平上研究了Fe 原子与N2分子相互作用的单端位构型的直线形和弯曲形两种结构的平衡几何结构、电子结构、轨道布局及红外光谱等性质. 计算结果表明, 由于强的σ-σ电子对互斥作用, 基组态4s²3d⁶的Fe原子不能与N2分子发生化学作用; 当Fe 原子呈现可与N2之间发生σ-π授予反馈作用的激发组态时, Fe 与N2分子之间可形成稳定的结构; 在得到的多个电子态中, 能量最低的是直线形的1³∑^-, 比Fe(a⁵D)和N2(¹∑⁺g )能量高21.6 kJ·mol-1, 同时存在几个能量相近的电子态, 如1³∏、1³Φ; 弯曲形都是不稳定态, 可能是连接直线形和单侧双配位构型的过渡态; 单端位构型产物相对于基态的反应物均是热力学不稳定的; 单端位构型中Fe对N2的活化作用很小, N—N 键长增加不超过7 pm.     

A CNDO/INDO molecular orbital formalism for the elements H to Br. applications

The CNDO/INDO molecular orbital formalism introduced in the preceding paper has been applied to a large number of atom combinations up to bromine under the inclusion of the first transition metal series. The results are compared with experimental data (geometries, ionization potentials, dipole moments) or with the results of sophisticatedab initio calculations (one electron energies, net charges, atomic populations). The semiempirical model reproduces for a wide range of molecules the experimental andab initio data with remarkable success.     

Spin Polarization at Organic-Ferromagnetic Interface: E ect of Contact Con guration

Based on ab initio theory, the interfacial spin polarization of a benzene-dithiolate molecule vertically adsorbed on a nickel surface is investigated by adopting di erent microscopic con-tact con gurations. The results demonstrate a strong dependence of the interfacial spin polarization on the contact con guration, where the sign of spin polarization may vary from positive to negative with the change of contact con guration. By analyzing the projected density of states, an interfacial orbital hybridization between the 3d orbital of the nickel atom and the sp3 hybridized orbital of the sulfur atom is observed. We also simulated the interfacial adsorption in mechanically controllable break junction experiments. The magne-toresistance obtained from Julliere model is about 27% based on the calculated interfacial spin polarization, which is consistent with experimental measurement.     

Computational investigation of TiSin (n=2-15) clusters by the density-functional theory

The geometries, stabilities, and electronic properties of TiSin (n=2-15) clusters with different spin configurations have been systematically investigated by using density-functional theory approach at B3LYP/LanL2DZ level. According to the optimum TiSin clusters, the equilibrium site of Ti atom gradually moves from convex to surface, and to a concave site as the number of Si atom increases from 2 to 15. When n=12, the Ti atom in TiSi12 completely falls into the center of the Si outer frame, forming metal-encapsulated Si cages, which can be explained by using 16-electron rule. On the basis of the optimized geometries, various energetic properties are calculated for the most stable isomers of TiSin clusters, including the average binding energy, the highest occupied molecular orbital and lowest unoccupied molecular orbital (HOMO-LUMO) gap, fragmentation energy, and the second-order difference of energy. It is found that at size n=6,8,12 the clusters are more stable than neighboring ones. According to the Mulliken charge population analysis, charges always transfer from Si atoms to Ti atom. Furthermore, the HOMO-LUMO gaps of the most stable TiSin clusters are usually smaller than those of Sin clusters.     

共轭烯烃的分子几何构型、电子结构和生成热的DMM和AM1、PM3的比较研究

运用Ｄｅｌｆｔ分子力学（ＤＭＭ）力场和程序以及半经验分子轨道ＡＭ１和ＰＭ３方法计算研究了丁二烯、苯、甲苯、联苯、苯乙烯、富烯、、环辛四烯、［２，２］对环烷和菲等１０个共轭烯烃分子的几何构型、电子结构和生成热．ＤＭＭ计算的几何构型和生成热与实验结果相吻合，电荷分布结果与从头计算结果较接近．ＡＭ１和ＰＭ３计算的几何构型较好，但计算的生成热与实验结果偏差较大．ＰＭ３计算值比ＡＭ１的稍好．     

Multilayer adsorption of water at a rutile TiO2)(110) surface: towards a realistic modeling by molecular dynamics

We present a model combining ab initio concepts and molecular dynamics simulations for a more realistic treatment of complex adsorption processes. The energy, distance, and orientation of water molecules adsorbed on stoichiometric and reduced rutile TiO(2)(110) surfaces at 140 K are studied via constant temperature molecular dynamics simulations. From ab initio calculations relaxed atomic geometries for the surface and the most probable adsorption sites were derived. The study comprises (i) large two-dimensional surface supercells, providing a realistically low concentration of surface oxygen defects, and (ii) a water coverage sufficiently large to model the onset of the growth of a bulk phase of water on the surface. By our combined approach the influence of both, the metal oxide surface, below, and the bulk water phase, above, on the water molecules forming the interface between the TiO(2) surface and the water bulk layer is taken into account. The good agreement of calculated adsorption energies with experimental temperature programmed desorption spectra demonstrates the validity of our model.     

An ab initio study on thermal rearrangement reactions of 1-silylprop-2-en-1-ol H3SiCH(OH)CH=CH2

The thermal rearrangement reactions of 1-silylprop-2-en-1-ol H3SiCH(OH)CH=CH2 were studied by ab initio calculations at the G2(MP2) and G3 levels. The reaction mechanisms were revealed through ab initio molecular orbital theory. On the basis of the MP2(full)/6-31G(d) optimized geometries, harmonic vibrational frequencies of various stationary points were calculated. The reaction paths were investigated and confirmed by intrinsic reaction coordinate (IRC) calculations. The results show that the thermal rearrangements of H3SiCH(OH)CH=CH2 happen in two ways. One is via the Brook rearrangement reactions (reaction A), and the silyl group migrates from carbon atom to oxygen atom passing through a double three-membered ring transition state, forming allyloxysilane CH2=CHCH2OSiH3. In the other, the reactant undergoes a dyotropic rearrangement; the hydroxyl group migrates from carbon atom to silicon atom coupled with a simultaneous migration of a hydrogen atom from silicon atom to carbon atom, forming allylsilanol CH2=CHCH2SiH2OH (reaction B). The barriers for reactions A and B were computed to be 343.5 and 203.7 kJ/mol, respectively, at the G3 level. The changes of the thermodynamic functions, entropy (DeltaS), entropy (DeltaS(doubledagger)) for the transition state, enthalpy (DeltaH), and free energy (DeltaG) were calculated by using the MP2(full)/6-31G(d) optimized geometries, and harmonic vibrational frequencies of reactants, transition states, and products with statistical mechanical methods, and equilibrium constant K(T) and reaction rate constant k(T) in canonical variational transition-state theory (CVT) with centrifugal-dominant small-curvature tunneling approximation (SCT) were calculated over a temperature range 400-1300 K. The conventional transition-state theory (TST) rate constants were also calculated for the purposes of comparison. The influences of the vinyl group attached to the center carbon of the alpha-silyl alcohols on reactions were discussed.     

某些含氧酸中原子净电荷与pK_a的关系

在量子化学从头算 (分别采用HF/STO 3G ,HF/STO 3G ,HF/ 3 2 1G和MP2 /STO 3G)水平下 ,对 11种含氧酸分子 (H2 O ,H2 O2 ,HNO2 ,HNO3 ,HClO ,HClO2 ,HClO3 ,HCO2 H ,CH3 CO2 H ,HBrO和HIO)进行几何构型全优化和电子结构计算 ,将获得的电子结构数据与 pKa 实验值进行逐步回归分析 ,结果发现含氧酸中羟基氢原子净电荷、含氧酸中与羟基氧相连非氢原子净电荷两项与 pKa 呈良好的二元线性关系 ,表明含氧酸中羟基氧所连原子的净电荷在决定含氧酸的酸性上起重要作用     

Comparative study of imidazole hydration: Ab initio and electrostatic calculations vs. Cambridge structural database analysis

We studied geometries and energies of complexes between water and neutral or protonated imidazole by ab initio molecular orbital calculations using the 4-31G basis set with and without the counterpoise correction. Positions of hydration sites and relative binding energies could be also estimated by using the electrostatic field map of imidazole as calculated by our bond increment method. The reliability of the calculations is confirmed by comparing the geometries of the imidazole-water complex to the experimental ones from the Cambridge Structural Database. These were obtained by X-ray diffraction studies on crystals with water bound to a molecule containing the imidazole fragment.     

银表面甲醇氧化反应机理的ab initio计算研究 1: 银表面静态吸附物种的几何构型和吸附性质

本文用原子簇模型(CM)的从头计算方法, 计算了银表面甲醇氧化反应中的静态吸附物种的优化几何构型及吸附性质。计算表明在清洁银表面甲醇、甲醛只存在物理吸附; 当表面存在吸附氧原子时, 甲醇可在银表面形成两种分子态吸附;甲醛与表面羟基OH(a)或氢原子H(a)共存时在银表面能够形成化学吸附, 且CH2O(a)极易与O(a)反应生成深度氧化中间体η^2-甲二氧基; 中间产物甲氧基在无氧的银表面能够形成稳定吸附, 在富氧银表面极易进一步氧化脱氢生成产物甲醛。通过计算与实验结果的对照, 我们对反应机理作了初步讨论。     

甲醛在CeO2(111)表面吸附的密度泛函理论研究

采用基于第一性原理的密度泛函理论和周期平板模型, 研究了甲醛在以桥氧为端面的CeO2(111)稳定表面上的吸附行为. 通过对不同覆盖度, 不同吸附位的甲醛吸附构型、吸附能及电子态密度的分析发现, 甲醛在CeO2(111)表面存在化学吸附与物理吸附两种情况. 化学吸附结构中甲醛的碳、氧原子分别与表面的氧、铈原子发生相互作用, 形成CH2O2物种; 吸附能随着覆盖度的增加而减小. 与自由甲醛分子相比, 物理吸附的甲醛构型变化不大, 其吸附能较小. 利用CNEB(climbing nudged elastic band)方法计算了甲醛在CeO2(111)表面的初步解离反应活化能(约1.71 eV), 远高于甲醛脱附能垒, 这与甲醛在清洁CeO2(111)表面程序升温脱附实验中产物主要为甲醛的结果相一致.     

Modeling the global potential energy surface of the N + N2 reaction from ab initio data

A new global potential energy surface for the N + N2 exchange reaction has been built from ab initio data. To overcome the difficulty of carrying out ab initio calculations for a large set of geometries the alternative strategy of fitting the minimum energy paths of the surface and their angular dependence using a modified LAGROBO functional form has been adopted. In this way we have been able to reproduce all the main features of the potential using a fairly small set of ab initio values.     

L-缬氨酸旋光异构的两种光反应可能途径

利用密度泛函理论(DFT)和从头算分子轨道理论对L-缬氨酸的旋光异构光反应机理进行了研究. 分别用B3LYP和MP2方法在6-311++G(d, p)基组级别上全优化得到了S₀和T₁态反应路径上的反应物、产物、中间体以及过渡态结构的几何构型, 给出了反应能垒, 利用含时密度泛函理论(TD-DFT)中的B3LYP/6-311++G(d, p)方法优化得到了S₁态反应路径上的平衡态几何构型. 通过分析反应途径上各个驻点的几何构型特征, 确定了L-缬氨酸在激发态可能通过手性碳上的氢原子以羰基氧或氨基氮为中转媒介发生质子迁移来完成旋光异构反应. 进一步用自洽反应场理论中的极化连续模型(PCM)方法探讨了溶剂化效应对旋光异构反应机理的影响.