分子动力学与原子多体势.doc

一．分子动力学简介 随着纳米科技的到来，许多新的学科产生了，例如纳米电子学、纳米生物学、纳米材料学、纳米机械学等。人们的注意力逐渐从宏观物体转向小尺度及相应的器件，其中微机械系统 (mieromachine)或称微型机电系统(mieroe一eetro一meeh耐ealsystem，MEMs)尤其取得了成功，并正被拓展应用于各种工业过程。由图 1可知，分子动力学正是处于尺度下的研究方法。 1957年Alder和ainwright[1]开创了分子动力学(Moleeularnynamies，MD)方法，之后经过多位科学家的努力，拓展了分子动力学方法的理论、技术及应用领域，尤其是在20世纪80年代由Andersen等先后完成的恒温、定压分子动力学方法，标志着分子动力学方法的科研应用进入了一个新阶段。 分子动力学方法是研究纳米尺度物理现象的重要手段。随着越来越多的材料原子间作用势函数被精确描述并经过实验验证、计算机硬件水平的快速更新以及高效率新算法的提出，分子动力学模拟被广泛应用于纳米尺度力学行为和纳米材料力学性能的研究。 在纳米尺度下，材料由离散的原子排列而成，由于比表面积大、表面效应明显，材料的力学性能和力学行为将与宏观材料迥异。基于连续性假设的宏观连续介质理论在研究材料的损伤演化、失效过程时，往往在时间和空间上将原子尺度的缺陷进行平均化处理，但这种处理仅适用于大量缺陷分布在材料中计算区域的情形，而对许多细微观材料和力学实验观测到的现象都无法解释，如疲劳与蠕变过程中的位错模式、塑性变形的不均匀性、脆性断裂的统计本质、尺寸效应等。因此，连续介质理论显然难以准确求解纳米尺度的力学问题。同时，如果直接从第一原理出发进行计算，除了类氢原子以外其他材料的薛定愕方程求解难度都太大，而且局域密度泛函近似理论并不是总能满足实际问题的需要。另一方面，材料本身在空间、时间和能量等方面存在藕合和脱祸现象，直接从头开始的量子力学计算难以很好地应用到几百个原子以下的计算规模中，无法达到一般纳米材料和器件的模拟要求。此外，由于实验条件控制的困难和合成、制备方式不同，各种纳米材料力学性能的有关实验结果分散性较大甚至相反，以至于目前难以通过纳米力学实验得到普遍适用的定量力学规律。 鉴于理论和实验上的困难，通过分子动力学方法模拟纳米尺度的力学性能和行为来探索纳米尺度的一般规律，是进行纳米力学研究的有效方法。分子动力学最早用于热动力学和物理化学，计算不同物理系统如固体、气体、液体的整体或平均热化学性能。1957年Alder首次提出并采用分子动力学方法分析刚性球系统的固液相变问题取得成功，此后，分子动力学开始逐渐应用于材料领域。随着上世纪80年代计算机硬件水平的提高和各种描述原子间作用的势函数的提出，分子动力学模拟日益活跃。通过分子动力学模拟不仅可以深入了解复杂的机制，发现本质上崭新的现象，而且可以定量模拟真实固体中所发生的过程，是诸如表面结构和扩散中的动力学和稳定性等许多结果的唯一来源。在EAM理论逐渐成熟和Baskes实验室、Ackland实验室等精确测出大量常用材料的EAM参数以后，分子动力学方法在模拟材料的物理性能和现象方面逐渐显示出强大的计算能力和较高的精度，大量的模拟尤其是固体结构、位错运动、表面界面现象、力学性能、变形机制和流体中的电泳、电渗透流等都得到了理想的结果。只要能将基于物理的模型建立起来，通过分子动力学计算就可以揭示出物理现象的本质，逐渐被广泛应用于凝聚态物理、材料学、力学、生物学、微电子学和微纳米加工等领域。 晶体是由大量的原子有序排列而成，材料的强度来源于原子间的相互作用，塑性来源于原子间的相互运动。因此，直接从原子尺度上对材料的微观力学行为进行研究显得非常有必要。分子动力学模拟技术既能得到原子的运动轨迹，还能像做实验一样观察。对于平衡系统，可以在一个分子动力学观察时间 (ObservationTime)内做时间平均来计算一个物理量的统计平均值，对于一个非平衡系统过程，只要发生在一个分子动力学观察时间内(一般为1一10ps)的物理现象也可以用分子动力学计算进行直接模拟。特别是许多与原子有关的微观细节，在实际实验中无法获得，而在计算机模拟中可以方便的得到。这种优点使得分子动力学方法广泛运用于材料科学与工程中，如材料设计、断裂分析等。 近年来，利用计算机模拟技术研究材料的力学性能日益成为人们感兴趣的课题。由于计算机处理速度的迅速提高，计算机模拟已经和实验观察、理论分析并列成为本世纪科学研究的三种方法。计算机模拟的数据(从模型中得来)可以用来比较、验证各种近似理论;同时，计算机模拟方法还可用来对实验和模型进行比较，从而提供了评估一个模型正确与否的手段。计算机模拟方法还有一个优点就是可以沟通理论和实验。某些量或行为可能是无法或难以在实验中测