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