毕业论文：《原子的距离和势能之间的关系研究》.docVIP

本 科 毕 业 论 文 第 PAGE 3 页 共 NUMPAGES 24 页 1 引言 分子动力学它是一门结合物理，数学和化学的综合技术。分子动力学是一套分子模拟方法，该方法主要是依靠牛顿力学来模拟分子体系的运动，以在由分子体系的不同状态构成的系综中抽取样本，从而计算体系的构型积分，并以构型积分的结果为基础进一步计算体系的热力学量和其他宏观性质。假设N为模拟系统的原子数,第i个原子的质量为mi,位置坐标向量为ri,速度为： 加速度为： 受到的作用力Fi,原子i与原子j之间距离为rij=|ri-rj|,原子j对原子i的作用力为fij,原子i和原子j相互作用势能为φ(rij),系统总的势能为： 所有的物理量都是随时间变化的,即A=A(t),控制方程如下: 为了在原子级别上研究体系的特性，对于大型体系，人们通常使用经典的分子动力学方法。在研究过程中，人们往往都要实现估计运行结果、沿着一定的方向试探计算体系的特征。这就需要对于微观原子间的作用力有一定的认识，才能比较准确地把握探索的方向。但目前为止，人们对于原子的作用力性质并无专门研究论文。因此，本次课题选择在分子动力学中应用较多的GEAM作用势，研究不同状况下，原子之间作用力的特性，试图揭示原子间作用力的一些特征。本次课题的意义如下：由于分子动力学模拟是在真实的的自然过程中，需要使用作用势。然而不同的作用势有不同的特征，不同的原子也有不同的特征。正是这种不同的特征导致了物理体系有不同的特征。目前人们对作用势的研究还不是很多，只是把作用势应用到不同的物理体系。如果事先有一点估计，且对分子动力学模拟有一定了解，则在计算前就可能更好地使模拟过程变得更简洁，更准确，为以后的计算提供一定的基础。 1. 1分子动力学方法 随着计算机模拟技术的发展,使实验上尚无法获得或很难获得的大量重要信息的获取成为可能,虽不能完全代替实验但为科研工作者们提供了重要的参考和指导实验、验证某些理论假设,降低试验的盲目性、成本低廉广等特点。分子动力学模拟(molecular dynamics simulation , MD)是在评估和预测材料结构和性质方面模拟原子和分子的一种物质微观领域的重要模拟方法,通过计算机对原子核和电子所构成的多体体系中的微观粒子之间相互作用和运动进行模拟,在此期间把每一原子核视为在全部其他的原子核和电子所构成的经验势场的作用下按照牛顿定律进行运动,进而得到体系中粒子的运动轨迹,再按照统计物理的方法计算得出物质的结构和性质等宏观性能. 简而言之即是应用力场及根据牛顿运动力学原理所发展的一种计算机模拟方法.分子动力学模拟是一种非常有效的计算机技术已成为重要的科学研究的方法之一.要想理解什么是分子动力学模拟,就必须首先清楚地理解力场等基本概念. 同时,在分子动力学模拟领域中,系综、周期性边界条件、积分步长等也是经常提及的术语名词,对它们的正确理解也影响着对分子动力学的深入理解. 1.1.1 GEAM作用势的引入 本次课题选用GEAM作用势作为研究对象，即普适嵌入原子势。之所以选择GEAM作用势是由于这种势目前的研究较少，很多的领域尚未涉足。本次研究能给这方面的进一步研究奠定一定的实验和理论基础。通常来说分子势能是很复杂的，可以有很多的表示方法，可以用势函数来表示，不同的材料它们的势函数也不一样。常用的有LJ势，金属多用EAM势，多原子分子如高分子材料的势函数还与键长、键角、双面角等因子有关。钨、铅、铜等原子之间的作用势则是典型的GEAM势。因此对GEAM作用势的充分理解是这次研究成败的关键。在EAM模型中，由N个原子组成的系统总能量可表示为: 其中： 其中φij是原子之间的两体作用对势，Fi是把原了i埋入电了密度为门的位置处所需要的能量.其中为系统中所有其他原了在i原了处产生的局域背景电子密度，ρ(rij)是第j个原子的核外电子在第i个原子处贡献的电荷密度。势中的第一项代表不同原子芯之间的排斥作用,第二项描述由于原子嵌埋于电子密度中(电子键合)引起的相互作用，可见这一项描述了势能的多体性质。 Daw和Baskes的EAM势的模型参数不是分析形式的，不能直接和具体的物理参数相联系，只适应处理特定的系统，对合金的计算则需用合金特性参数进行拟合，这就使EAM理论的应用受到很人限制，特别不适合用来系统地处理有关合金方面的问题。于是X.W.ZHOU等人提出了一种普适的EAM势（generalized EAM potential，简称GEAM）。GEAM是一种普适的作用势，适用于16(Cu,Ag,Au,Ni,Pd,Pt,Al,Pb,Fe,Mo,Ta,W,Mg,Co,Ti,和Zr)组成的合金。[21] 他们导出的对势采用如下形式： rc是平衡间距，A,B,α,β是四个参数，κ,λ是截断半径