高级从头计算法精确研究基态He和Li原子间的相互作用势

高级从头计算法精确研究基态Ｈｅ和Ｌｉ原子间的相互作用势王一波（贵州大学化学系,贵阳,５５００２５）关键词ＨｅＬｉ,弱相互作用势,ＭＰ４ＳＤＴＱ,ＱＣＩＳＤ（Ｔ）,ＣＣＳＤ（Ｔ）基态Ｈｅ和Ｌｉ原子之间的ｖａｎｄｅｒＷａａｌｓ作用势处于１０－６ａ,ｕ．数...     

价键理论中的组态相互作用

近 2 0年来 ,从头计算水平的价键 (VB)方法得到了人们的重视 ,并广泛应用于化学反应等问题的研究[1～ 5] ,然而目前价键理论的计算方法仍然很不完善 .用 VBSCF方法进行计算虽然比较简单 ,能正确地描述化学反应的形成机理 ,但数值结果不理想 ;而用 BOVB方法[4 ] 进行计算虽然可以得到较好的计算结果 ,但存在收敛困难等问题 .分子轨道理论中的组态相互作用是一种简单直接的电子相关能计算方法 ,显然这一方法可以应用于价键方法中 .然而与分子轨道理论方法不同 ,在价键方法中 ,无法直接得到空轨道 ,此外如何选取激发价键函数使得计算结果…     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

分子间相互作用的Ab Initio计算研究——氢键复合物体系(HF)_2中的非键相互作用

本文以6-311 G~(**)为基组对氢键复合物(HF)_2进行了Ab Initio计算,给出了2个HF在距离及取向平衡点附近的势能面。用势诱导最小二乘拟合法计算了体系中各原子的原子电荷。用(exp-6-1)势函数分析法分解了2个HF间的氢键相互作用能,表明其中主要的相互作用是电荷转移而不是静电性的。2个HF分子间的势能曲线呈Morse函数性质。     

反散射理论导出He-Ne相互作用势

以Ｈｅ原子束与Ｎｅ原子束交叉碰撞实验数据为基础。以反散射方法为手段，理论上导出Ｈｅ－Ｎｅ原子相互作用势，所得的势呈显出ｖａｎ ｄｅｒ Ｗａａｌｓ势的特征，即短程为强排斥，中程为碗形吸收，远程为零，这样的势能正确的，符合实验的弹发型用射微分截面。     

CH3F分子间相互作用的ab initio研究

在MP2/6-311++G(3d,3p)电子相关校正水平上, 对CH3F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证, 发现3种势能面上有极小点的构型. 进一步在高级电子相关校正的MP4SDTQ、 CCSD(T)/6-311++G(3df,3pd)方法水平上, 对其中总能量最小的构型进行精确计算, 得到二聚体的结合能为-9.707 kJ/mol. 研究结果支持了由光谱实验结果推测的构型, 解释了CH3F二聚体的谐振频率的多样性.      

价键结构函数与键函数 Ⅰ.理论处理及甲烷的应用实例

采用键函数形式对键表进行逐级展开,其中键函数的强正交积作为零级波函数,进而构造了各级校正函数以描述多对键之间的相关作用.同时采用Heitler-London函数作为键函数,对甲烷进行了键表计算,结果表明CH_4中一、二、三级校正能近似有4:2:-1的关系,即校正函数的收敛性能较好.     

规则星形高分子之间的相互作用势的自洽场研究

用自洽场理论方法计算了星形高分子之间的相互作用势, 并与标度理论的推测以及基于粒子的分子动力学模拟方法的结果进行了比较. 自洽场理论计算得到的星形高分子之间的二体相互作用势与标度理论推测的结果基本一致: 当两个星形高分子核心之间的距离小于其对应的等效软球的半径R时, 二体相互作用势U与星形高分子核心之间的距离d及官能度f之间的关系为U(d)∝-f 3/2 ln(d/R). 当星形高分子核心间距离较小, 相互作用势与星形高分子核心间距离的对数呈线性关系时, 三体相互作用为吸引作用, 但其强度约等于二体作用的10%. 因此, 以二体作用为基础的基于粒子的模拟方法是基本合理的.     

叠氮乙烷二聚体分子间相互作用的理论研究

Energetic materials are aggregative and mixed systems. The intermolecular interactions play significantroles in the physical,chemical and explosive property. The study on intermolecular interactions of energetic materials has attracted wide attention. The organic azides are an important category of energetic materials and widely used in many fields. Ethyl azide is the simple model having the explosive property for the organic azides energetic compound. Ethyl azide monomer（Ⅰ）and all its possible stable clusters（Ⅱ,Ⅲ and Ⅳ）are fully optimized by ab initio method at the HF/6-311++G** level. Vibrational frequencies calculated to ascertain each structure are characterized to be the stable structure（no imaginary frequencies）. The proportions of correlated interaction energies to their total interaction energies ΔE（MP2）are 65.14%,63.76% and 65.62% for Ⅱ,Ⅲ and Ⅳ respectively. In addition,the basis set superposition error（BSSE）correction energies are 7.82,7.61 and 4.40 kJ/mol for Ⅱ,Ⅲ and Ⅳ respectively. The zero point energy （ZPE） corrections for the interaction energies are much less than those of MP2 electron correlation and BSSE correction energies. After MP2 electron correlation correction,BSSE and ZPE correction,the greatest corrected intermolecular interaction of the dimers is -10.45 kJ/mol. The charge redistribution mainly occurs on the adjacent N?H atoms between submolecules. The charge transfer between two subsystems is very small. Natural bond orbital（NBO）analysis is performed to reveal the origin of the interaction. Based on the statistical thermodynamic method,the standard thermodynamic functions,heat capacities（C0p）,entropies（S0m）and enthalpies（H0m）and the changes of thermodynamic properties from the monomer to dimer with the temperatures ranging from 200. 00 K to 800. 00 K have been obtained.     

C_（60）分子间相互作用势及其晶体压缩系数的计算

根据Ｃ（６０）晶体升华焓的实验值求得了Ｃ（６０）分子间相互作用势的解析表达式。用该表达式计算了Ｃ（６０）分别为面心立方和六方晶格时晶体的结合能，讨论了不同晶格的稳定性问题，解释了Ｃ（６０）晶体一般是面心立方的原因。本文还根据Ｃ（６０）分子间相互作用势计算了Ｃ（６０）晶体的压缩系数，并与文献的实验值进行了比较，所得结果与实验值相符合。     

Theoretical study on the mechanisms of the nucleophilic substitution reactions of hydrosulfide ion and halomethanes

Three species involved in the nucleophilic substitution reaction of hydrosulfide ion and halomethanes are investigated by ab initio calculations. Geometries for stationary structures along the reaction paths are fully performed with the second‐order Møller–Plesset perturbation approximation with the cc‐pVDZ basis set. The monomer geometries determined by the MP2 method match the experimental results very well. Single point energy calculations are carried out at the coupled cluster with perturbative triple excitations CCSD (T) theory with aug‐cc‐pVDZ basis set. Halomethanes have three conformers here, which lead to the three product channels, HSCH₃ + F^?1, HSCH₃ + Cl^?1, and HSCH₃ + Br^?1. The investigation encompasses the six complexes formed among three channels, respectively. By selecting the six complexes as the model, we investigate the binding energy, topological property of the electron charge density and their Laplacian in detail theoretically. Electrostatic density potential maps of halomethanes are generated for the determination of attractive interaction sites. It is proved that the similar misshaped electron clouds of the three halogen atoms result in the similar properties of the carbon‐halogen bonds, and reveals that the product ion‐dipole complexes interactions are predominantly electrostatic in nature. The calculated results predict the binding energy of the most stable complex in six complexes is ?47.06 kcal/mol at the MP2 level of theory. The second channel has the lowest energy barrier, which is ?3.63 kcal/mol at the CCSD (T) levels of theory, is expected to be the most important pathway. It occurs via C? Cl cleavage accompanied by C? S bond formation. The other two channels have higher energy barriers and are expected to have smaller rates. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

He2F-体系的分子间相互作用势能面

使用高水平的从头算CCSD(T)/aug-cc-pVTZ方法, 经过Counterpoise校正, 计算了He2F-体系的分子间相互作用势能面. 在He2F-体系的相互作用势能面的最小值处, 发现了一个等腰三角形的稳定结构. 在这个结构中, He…F- 距离是 0.334 nm, He…He 的距离是 0.295 nm, ∠HeF-He 为 52.5°. 计算了此稳定结构的频率、相互作用能、二体相互作用能和三体相互作用能. 在CCSD(T)/d-aug-cc-pVTZ水平下, 相互作用能为-1.727 kJ/mol.     

Rg-HX分子间势的精确从头计算研究

在用非迭代的三重激发项来校正CCSD的CCSD(T)理论水平下,采用aug cc pVQZ基函数对He HF的分子间势进行了系统的研究.结果表明：He HF以线型结构存在.在极限基的情况下,复合物两种线型极小点结构He H F和He F H势阱深分别为46.614 cm－1和25.026 cm－1,对应He原子到HF分子质心的距离Rm分别为0.3149 nm和0.3012 nm.讨论了不同的基函数和理论方法在研究此类弱束缚态复合物的分子间势时的可靠性及其对结果的影响,并研究了HF分子中H－F键长的改变对势能的影响,同时也给出了势函数的解析形式.     

The MIDI! basis set for quantum mechanical calculations of molecular geometries and partial charges

We present a series of calculations designed to identify an economical basis set for geometry optimizations and partial charge calculations on medium-size molecules, including neutrals, cations, and anions, with special emphasis on functional groups that are important for biomolecules and drug design. A new combination of valence basis functions and polarization functions, called the MIDI! basis set, is identified as a good compromise of speed and accuracy, yielding excellent geometries and charge balances at a cost that is as affordable as possible for large molecules. The basis set is optimized for molecules containing H, C, N, O, F, P, S, and Cl. Although much smaller than the popular 6-31G* basis set, in direct comparisons it yields more accurate geometries and charges as judged by comparison to MP2/cc-pVDZ calculations.     

亚甲基氰自由基(HCCN)与一氧化氮(NO)反应势能面的理论研究

在CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p)+ZPE水平上对反应HCCN+NO的二重态反应势能面进行了计算,得到了4种产物:P1(HCN+NCO),P2(OH+NCCN),P₃[HCN+(CNO)]和P4(HCN+CNO).其中产物P1为主要产物,P2为次要产物,P₃和P4很难得到.在G2(B3LYP/MP2/CC)水平,对产物P1和P2的反应通道的单点能量进行了校正.     

Ab Initio MO Studies on the Reaction Mechanism for Carbonyl Insertion Catalyzed by Carbonyl Cobalt Complex

IntroductionSome of ah l'nito MO studies on elementary reactions catalyzed by organometallic compounds have been published[' 41, which are of momentous significance for better understanding the catalytic action and developing the catalytic agent. At the same time, there still aremany unknown aspects that remain to be studied. Carbonyl insertion is one of the key reac.tlons, and it can go in three possible pathsL'--'j, i. e., carbonyl migration, ethyl migrationand concerted reaction paths. Th…     

Bond functions and the topological properties of the bonds

The effect of bond functions on the topological properties of the bonds is studied by performing Hartree-Fock calculations for O₂, Cl₂, F₂, CN^–, N₂ and LiCl molecules. The results show that bond functions increse the degree of covalency of the bonds and the electronic density at the (3, –1) bond critical point. This effect is quite noticeable for intermediate interactions molecules (F₂ and CN^–) where a is obtained instead of a positive value. The optimal position of the bond functions seems to be at or near the (3, –1) bond critical point.     

第一过渡系金属硅烯配合物的理论研究

以HF/6-311+G^*基组研究了硅烯SiH₂同第一过渡系金属的配合物MSiH₂的分子轨道特征及键解离能.MSiH₂为共平面构型.其中基态的3TiSiH₂和4CoSiH₂带有明显的双键特征.M-Si键具有共价性质.M-Si的键解离能,从Sc到Cu呈现周期性变化,这种变化趋势同M的金属离子激发能之间存在近似的线性关系.     

Ab Initio energetics of Si?O bond cleavage

A multilevel approach that combines high‐level ab initio quantum chemical methods applied to a molecular model of a single, strain‐free Si O Si bridge has been used to derive accurate energetics for Si O bond cleavage. The calculated Si O bond dissociation energy and the activation energy for water‐assisted Si O bond cleavage of 624 and 163 kJ mol⁻¹, respectively, are in excellent agreement with values derived recently from experimental data. In addition, the activation energy for H₂O‐assisted Si O bond cleavage is found virtually independent of the amount of water molecules in the vicinity of the reaction site. The estimated reaction energy for this process including zero‐point vibrational contribution is in the range of −5 to 19 kJ mol⁻¹. © 2017 Wiley Periodicals, Inc.