He-LiH势能面的ab initio计算

采用超分子耦合簇理论CCSD(T)方法和由键函数3s3p2d1f组成的大基组,计算得到了He-LiH体系在冻结LiH键长情况下的全程势能面.该势能面存在2个势阱,较深的势阱对应于线性He-LiH构型,其中Rm为4.25a0,阱深为177.53cm-1 ,较浅的势阱对应于线性He-HLi构型,其中Rm为9.875a0,阱深仅为9.90cm-1,整个势能面表现强的各向异性.     

Ne-HCN体系势能面的ab initio计算

在单双迭代包括非迭代三重激发微扰处理的耦合簇CCSD(T)理论水平下,采用aug-ccPVQZ并加上键函数(3s3p2d2f1g)组成的大基组,计算得到了体系在冻结HCN键长情况下的全程势能面.计算结果表明该势能面存在2个势阱:较深的势阱对应于线性Ne-HCN构型,阱深为61.56cm-1,其中Rm=8.00a0;较浅的势阱对应于T型构型,其中Rm=6.62a0,θ=102°,阱深为52.48cm-1,整个势能面呈现出弱的各向异性.     

铬铁矿的分子动力学势能参数反演

导出一组对铬铁矿及其高压相、分解产物都适用的Buckingham势能参数,用于铬铁矿高压相变、分解以及高压相卸压降温过程的分子动力学模拟.首先,利用密度泛函理论(DFT)计算获得铬铁矿、钙铁相、钙钛相、FeO和Cr2O3不同体积时的结构和总能量.再根据晶格能与总能量之间的关系,由DFT所得结构和总能量,列出势能参数所满足方程组.应用遗传算法这一全局优化算法获得势能参数.为检验所得势能参数的正确性,将其用于分子动力学计算,得出的5种晶体的晶格常数与实验和DFT结果分别差～2%和～1%,得出的5种晶体的V(p)关系与DFT拟合的状态方程曲线相吻合.     

MEHMO法计算F+H_2共线反应势能面

用改进的EHMO(MEHMO)法计算了F+H_2共线反应的势能面,本方法在H矩阵的对角元中引入了一项描述不同原子间排斥力的经验项,基组选用F原子的2p_(?)轨道和H原子的1s轨道,参数值B为0.64,在得到的势能面上,于R_(UH)~*=2.54 a.u.和P_(HH)~*=2.18a.u.处出现的位垒高度为15.72kJ/mo1.反应产物和作用物的能量差为-144.21kJ/mo1.本文所得结果与实验测定、从头计算法的结果一致。     

H原子通过碳纳米管的分子动力学模拟

利用Brenner(#2)半经验多体相互作用势和分子动力学模拟方法研究了氢原子在单壁碳纳米管中的运动.对温度为10K时,初始入射动能分别为1,5,10,15,20,30eV的H原子以不同的入射角度进入碳纳米管进行模拟,观察了氢原子与碳纳米管的吸附情况及在碳纳米管中的运动,并讨论了不同入射能量和入射角度对氢原子在碳纳米管中吸附的影响.结果表明,在能量为20eV时H原子能够很好的吸附在碳纳米管管壁上.     

势能面—过渡状态理论的物理模型

本文较详细地讨论了过渡状态的势能面物理模型。首先从最简单的H—H—H体系的势能面入手,分析了H+H_2反应的详细微观机理及能量变化情况,并且图示了该反应势能面的主要特点,同时阐明了活化能概念的本质,着重讨论了活化络合物的形成与分解方式,然后对较复杂的排斥势能面及吸引势能面进行了探讨,最后指出了该物理模型的优缺点及应用。     

从头计算分子动力学的理论基础

本介绍了基于密度泛函理论的从头计算分子动力学的理论基础,具体推导了HatreeFock方程和Kohn-Sham方程,并给出了运用从头计算分子动力学方法进行现代新材料的设计以及计算晶体材料电子结构的流程图。     

BEBO法计算活化势能的研究

重点研究了Ｆ、Ｃｌ、Ｂｒ、与Ｈ２及ＣＨ４六个反应的活化势能，每个反应都计算出上百个数据，以ｎ对Ｕ势作能面剖面图，该图与文献数据符提出了依据。     

利用岩层三角面图解计算产状的新方法

现有的线性投影图解法求解岩层产状尚有某些局限性。该文经理论推导,在原方法的基础上得出三种新的方法。     

基态的HCN和HNC分子的势能面

运用密度泛函理论(DFT)的B3LYP方法在6-311++G**水平上,对基态HCN和HNC分子的结构进行了优化计算,得到HCN分子的稳定结构为C∞v构型,电子态为X1∑+u,平衡核间距RH-C=0.1066 nm、RC-N=0.1149 nm,离解能De=18.88 eV; HNC分子的稳定结构为C∞v构型,电子态为X1∑+u,平衡核间距RH-N=0.09996?nm、RN-C=0.1169?nm,离解能De=18.256?eV,用多体项展式理论推导了基态HCN和HNC分子的解析势能函数,其等值势能图准确再现了基态HCN和HNC分子的结构特征及其势阱深度与位置.     

过冷态水表面性质的分子动力学研究

采用SPC/E 和TIP4P 两种势函数, 分子动力学模拟了228~293 K 温度范围内水的表面性质. 模拟结果表明, 表面张力随温度的降低而增大. 当温度低于273 K 时, 表面张力的温度系数开始缓慢升高, 但在228 K 附近没有观察到类似某些体热力学参量的奇异发散行为, 这与实验结果相一致. 与实验测量值比较, SPC/E 势函数的模拟结果符合较好, 而TIP4P势低估了表面张力. 造成这种差异的主要原因可能是SPC/E 势较好地描述了过冷态水的表面结构. 在表面水分子的方向分布模拟中, SPC/E 势下的水分子表现出更强的方向有序性, 而且表面势也比TIP4P 势下的相应值高.     

指数平滑法初始值计算与平滑系数选取的新方法

通过迭代方法证明了指数平滑方法经过多次迭代后迅速逼近于全期指数加权平均法;并提出以递推算法最小数据量t*为界,分别采用迭代法(或全期指数加权平均法)(t<=t*时)和递推算法(t>t*时)来分段计算指数平滑值;在强调时序数列平稳性的基础上,提出了以最近k个数据的累计贡献率或以最近前若干期数据之加权系数衰减50%的两种估...     

多晶面银纳米线拉伸形变的分子动力学模拟

本文以低能量的四棱台结构为基本单元,构建了[100]晶向的多晶面单晶Ag纳米线.利用分子动力学方法模拟其在拉伸载荷下的形变行为,考察了纳米线的单元结构、体系温度、拉伸速率等因素对该纳米线机械强度和断裂模式的影响.结果表明:当四棱台结构单元的长边与短边比值小于1.4时{111}表面的强化作用占主导地位,而当比值大于1.4时多晶面的引入对纳米线起到弱化作用.在断裂模式上,多晶面纳米线在超低温和低应变速率下易于表现脆性断裂,而在室温和高应变速率易于表现韧性断裂.后者塑性形变的机理受位错成核、增殖和位错消失的共同作用,而前者的形变机理与位错活动无关,超低温使材料脆性增加,结构单元之间交界处产生裂隙,裂隙的形成和扩张导致了纳米线的迅速断裂.     

基于Lennard-Jones和Soft-Sphere势能模型的流体传递性质分子动力学模拟比较

利用分子动力学模拟计算传递性质方法，采用Green-Kubo关系式基于Lennard-Jones势能模型和Soft-Sphere势能模型模拟了速度相关函数、压力张量相关函数、热流相关函数以及它们对应的传递系数，讨论了Lennard-Jones势能分子间吸引力和Soft-Sphere势能“软度”对流体传递性质的影响．     

田菁胶分子的势能面和几何构型

用ＧＡＭＥＳＳ中ＭＯＰＡＣ的ＰＭ３程序在全局优化中确定了田菁胶大分子的几何构型和稳定构象,计算了该大分子绕两个旋转轴θ,φ旋转的势能面,结果表明大分子基闭之间不能绕θ,φ轴进行自由旋转。     

哈雷慧星运动轨道的分子动力学计算

基于分子动力学理论,对哈雷慧星运动轨道作了数值模拟计算.在此基础上进一步计算得到由于太阳质量损失引起的哈雷慧星运动轨道半长轴随时间的变化率为29.5 cm/a.     

高密度氢的从头计算分子动力学模拟

运用从头计算分子动力学方法模拟研究了两种密度（rs=0.62、1．24）下的高密度氢,确定了切断能与总能量的关系,并给出了离子和电子的互作用随距离的变化关系,最后得出当rs不大于0．62时高密度氢表现为等离子体相。     

力学体系保守力的判据及势能计算

得出了力学体系保守力的判据以及势能的计算方法。     

力学体系保守力的判据及势能计算

得出了力学体系保守力的判据以及势能的计算方法.