分子动力学在材料科学中的应用.docx

分子动力学在材料科学中的应用 摘要：分子动力学方法是进行物质原子或分子层次计算机模拟时所采用的一种基本方法。通过分子动力学模拟，可以给出原子尺度上材料及其演化过程细节的可能性，具有史无先例的准确性，使材料设计和性能的预测成为可能。本文介绍了分子动力学模拟的基本原理和算法。综述了用分子动力学模拟方法研究固休的相结构、表面问题、界面问题以及薄膜形成过程的方法。 关键词：材料;分子动力学;计算机模拟； Abstract： Molecular Dynamics method is one of the basic methods used to simulate the atomic or molecular level computer simulation. By molecular dynamics simulation，we can give the possibility of the details of the evolution process of the material on the atomic scale，and it has the possibility of no precedent，so that the material design and performance prediction can be made. This paper introduces the basic principle and algorithm of molecular dynamics simulation，and summarizes the methods of studying the phase structure，surface problem，interface problem and the formation process of the thin film by molecular dynamics simulation. Key words： materials; molecular dynamics; computer simulation; 分子动力学(Molecular Dynamics ，简称MD)用于计算以固体、液体、气体为模型的单个分子运动，它是探索各种现象本质和某些新规律的一种强有力的计算机模拟方法，具有沟通宏观特性与微观结构的作用，对于许多在理论分析和实验观察上难以理解的现象可以做出一定的解释。MD方法不要求模型过分简化,可以基于分子(原子、离子)的排列和运动的模拟结果直接计算求和以实现宏观现象中的数值估算。可以直接模拟许多宏观现象，取得和实验相符合或可以比较的结果，还可以提供微观结构、运动以及它们和体系宏观性质之间关系的极其明确的图象。MD以其不带近似、跟踪粒子轨迹、模拟结果准确而倍受研究者的关注，在物理、化学、材料、摩擦学等学科及纳米机械加工中得到广泛而成功的应用。本文主要讲述MD方法在材料科学中的应用及MD的必威体育精装版研究进展。 1、分子动力学基本原理 分子动力学将连续介质看成由N个原子或分子组成的粒子系统，各粒子之间的作用力可以通过量子力学势能函数求导得出，忽略量子效应后，运用经典牛顿力学建立系统粒子运动数学模型，通过数值求解得到粒子在相空间的运动轨迹，然后由统计物理学原理得出该系统相应的宏观动态、静态特性，图1所示是MD模拟过程。 MD具体的做法是计算机上求运动方程的数值解。通过适当的格式对方程进行近似，使之适于在计算机求数值解。从使用连续变数和微分算符的描述过渡到使用离散变数和有限差分算符的描述，显然会有误差，误差的阶数取决于具体的近似机制，即所用的算法。模拟首先是规定初始条件。为了确定起见可令初始位置在格子的格点上，而初始速度则从波尔兹曼分布得出。一个按上述办法建立的系统不会具有所要的能量，而且很可能这个状态并不对应于一个平衡态。为了推动系统到达平衡，需要一个趋衡阶段。可以通过增加或从系统中移走能量，对运动方程向前积分若干时间步，使系统弛豫到平衡态。接着是物理量的计算阶段，沿着系统在相空间中的轨道计算一切令人感兴趣的量。 模拟中,MD采用周期边界条件和最小镜像原理,可以大幅度减少计算工作量。周期边界条件是将一定数量的粒子N集中在一定的容积V中，这个容积V称为原胞，原胞周围的部分可以看作是原胞的复制，它们称作镜像细胞。这些镜像细胞的尺寸和形状与原胞完全相同，并且每个镜像细胞所包含的N个粒子是原胞中粒子的镜像，原胞在各个方向周期复制便形成了宏观物质样本。这样只需根据原胞周围的边界条件计算原胞内粒子的运动幅度,减少了工作量。原子间作用势能模型的构造对于MD 法的应用至关重要。最简单的偶势模型只考虑两体作用，而与其它原子无关，在模拟中运算量小。20世