Si60团簇与Si（111）面碰撞的紧束缚分子动力学模拟

利用普适的紧束缚分子动力学模拟研究在不同入射条件下Si60团簇与Si(111)面的碰撞机理，结果表明Si60分子是“非弹性分子”，入射角度对碰撞结果有重要影响，Si60分子易和Si(111)面发生化学反应，其保持结构特性吸附于Si(111)表面的能量域值在50～80eV之间．     

Si10团簇与Si(111)面碰撞的紧束缚分子动力学模拟

利用紧束缚分子动力学方法模拟研究了Si10团簇与Si(111)表面碰撞的微观过程．研究了入射动能在10—60eV之间的碰撞过程，给出了不同时刻体系的原子位置的微观状态图．结果表明：Si10与Si(111)面碰撞时，吸附在晶体表面的能量域值为20eV，替位域能为30eV，损伤域能为40eV；入射能量大于50eV时，晶体表面严重损伤并且在60eV时，发现在外延层有Si7团簇生成；Si10团簇保持其结构特性吸附在Si(111)表面，需要的最佳能量在10～20ev之间，进而通过对碰撞结果的讨论得到了改变轰击能量可以控制外延生长的结构的结论。     

Si在Si（111）表面吸附过程的分子动力学模拟

利用Tersoff半经验多体相互作用势和分子动力学模拟方法研究了荷能的硅原子在Si(111)表面的吸附过程.对100K时,初始入射动能分别为5,15,25,30eV的硅原子从6个不同位置轰击Si(111)表面进行了模拟,观察到了硅原子在Si(111)表面形成的吸附结构,并讨论了不同入射位置和入射能量对吸附结构的影响.结果表明,硅原子相对于表面不同的局部构型发生不同的碰撞过程,而硅原子入射能量的提高有利于成键过程的发生,从原子尺度模拟了沉积机制.     

紧束缚分子动力学研究硅团簇(n=20)的稳定结构

利用紧束缚分子动力学模拟了硅团簇(n=20)在不同条件下的稳定结构,分析并比较了它们的构型和结合能的大小,结果发现模拟退火温度在1700K 时得到的结构是最稳定的.     

硅团簇拉伸/压缩性能的分子动力学模拟

采用基于Stillinger-Weber(SW)势的分子动力学方法模拟了Si13团簇温度T分别为300K,50 K,5 K下的拉伸/压缩过程,得出了Si13团簇拉伸/压缩的力与应变的曲线图,讨论了不同温度下Si13拉伸/压缩的力学特性以及拉伸与压缩特性的差异。研究表明,拉伸/压缩初始阶段力与应变呈线性关系,拉伸时温度越高受力越大,压缩时正好相反,而且同等应变下拉伸时的力明显高于压缩时的力;当拉伸应变达到承载极限,继续拉伸,力突然变得很小,随后力与应变又呈现出与初始阶段相同的线性关系,这是团簇所特有的现象。     

紧束缚分子动力学模拟:C_(60)与石墨(0001)面的碰撞

运用紧束缚势的分子动力学模拟方法，研究了Ｃ６０分子对石墨表面的碰撞，观察在入射动能为２４０ｅＶ时碰撞的微观过程以及动能和势能的变化趋势．结果发现在弹回过程中Ｃ６０分子质心的运动可被看做是在准谐势下的运动，通过碰撞Ｃ６０分子在石墨表面形成了薄膜．     

团簇在金属表面沉积过程的分子动力学模拟

在冷喷涂过程中,喷涂粒子被超音速气流加速到较高的速度,在低于喷涂材料熔点的温度下撞击基体,发生剧烈的塑性变形而沉积形成涂层.但是由于高速粒子碰撞变形的瞬时特点,不能对粒子变形沉积过程进行直接观察,通过对Ni团簇在Cu基体上的沉积过程的分子动力学模拟,可以观察到团簇撞击基体并在基体上沉积的过程,以及团簇和基体的形貌变化;另外,通过计算团簇原子进入基体表面层的数量探讨了影响团簇沉积过程的主要因素.     

Si团簇的结构特性研究

采用由M．Menon等人的非正交紧束缚模型，结合遗传算法和分子动力学方法，对n=2-20的中小尺寸Si团簇的结构进行了优化，计算了团簇的束缚能Eb、平均配位数CN，裂化能Efrag，总能的二阶差分△2E，HOMO和LUMO的能隙Eg，得出了随着团簇尺寸的增加，构型发生了近球形-长椭球形-近球形的转变、n=4，7，10等团簇具有较强稳定性，团簇中原子之间共价键性结合占支配地位等主要结论。所得的主要结论与已有的从头计算和实验结果相一致。     

Fe团簇在Fe（110）表面上扩散和结构稳定性的分子动力学研究

采用嵌入原子势和分子动力学方法研究了Fe的小团簇（小于10个原子）在Fe（110）表面上的结构稳定性和扩散特性.计算结果表明体心立方结构的Fe表面,四原子和七原子团簇有较高的转变能和扩散势垒,说明这两种团簇比其他原子团簇更稳定和容易观测到.这两种团簇可能是薄膜生长过程中的临界晶核.研究了三聚体的旋转扩散和团簇扩散过程中二聚体剪切机制,发现二聚体扩散可能是Fe的小团簇在Fe（110）表面扩散的一种重要扩散方式.     

Si（111）表面离子束辅助沉积对类金刚石薄膜结构影响的计算机模拟

本文选C2分子和Ar离子作为沉积源和辅助沉积粒子,采用分子动力学（MD）方法在Si（111）面上模拟研究了离子束辅助沉积（IBAD）类金刚石（DLC）膜的物理过程.重点讨论了C2分子和Ar离子的入射能量及到达比（Ar/C）对平均密度和sp^3键含量的影响,并与Si（001）-（2×1）表面生长类金刚石膜的结果进行比较.结果表明,到达比和入射能的改变,对薄膜结构的影响不同;Si（111）面上生长类金刚石膜,薄膜在衬底的附着力更强.     

碳团簇结构的紧束缚理论方法

利用紧束缚分子动力学模拟退火方法研究了碳团簇Cn(n=2-8)的结构和能量,通过与前人工作结果的比较,发现本理论方法的结果与ab inito方法计算的结论相符。因此,用紧束缚分子动力学方法可对较大碳原子团簇进行计算。     

C60-C60分子对撞的计算机模拟

采用紧束缚方法模拟了两个 Ｃ６０分子对撞的过程,画出了对撞时的分子图形,并给出了不同能量下对撞得到的结果     

紧束缚分子动力学模拟纳米硅管的稳定结构

利用紧束缚分子动力学退火方法模拟研究了纳米硅管(SiNT)的稳定结构和基态能量,结果表明:几何结构特征对纳米硅管的结合能有重要影响,平均键长为0.236 4 nm,表面为双层原子面,具有很高的亚稳性,锯齿型管比扶手椅型管的原子结合能大0.051 eV/atom,SiNT(m,m)(m=2-4)不具有管状结构,SiNT(2,2)可以作为单元形成硅纳米线.     

C_(60)分子的旋转速度对对碰结果的影响

运用紧束缚分子动力学模拟的方法研究了两个C60 分子的旋转对碰 ,与无旋转情况下对碰的结果比较表明 ,旋转这一因素对碰撞结果起着重要的作用 .     

液态碲碘化合物结构的分子动力学模拟

利用分子动力学方法模拟液态碲碘化合物Te1-xIx的结构随着原子浓度的变化而改变．研究结果表明随着碲原子的浓度的增加，液态碲碘化合物经历了由非金属性到金属性的转变，其转变发生在Te0.6I0.4到Te0.8I0.2；偏配位数NTeTe，NTeI，NITe和NII在富碘、富碲的一侧有不同的变化规律．配位数、角分布函数和原子快照的分析显示随着碘原子浓度的增加，液态碲碘化合物Te1-xIx经历了由类似于液态碲的链状结构到类似于液态TeI4结构的转变．     

分子动力学模拟锐钛矿的压缩行为

作者利用壳层模型分子动力学对锐钛矿在室温下的压缩行为进行了研究.基于常压下的状态方程作者得到了锐钛矿的晶格常数、弹性常量和体变模量.与此同时作者还研究了晶格常数a和c,以及原胞体积和压强的关系.其计算结果和已有的实验结果符合得较好.