分子动力学课件.ppt

1、分子动力学分子动力学模拟分子动力学模拟Molecular dynamics simulation单学坤 彭瑶 王志强 刘志敏 张孝松 施娅NanJing Tect Unversity背景介绍 举例：n微电子电路尺度己达微米量级中的导热对流冷却n微槽道微米量级中的对流换热n薄膜材料的热传导n微电子机械系统中的流动与传热共同的特点 换热是在微尺度条件下进行的传统换热微尺度换热比较微尺度传热具有新的特点和规律：p微槽道中对流换热强度比常规尺度对流换热高出两个数量级p超薄薄膜材料导热系数比常规尺度同种材料导热系数低一到两个数量级p微尺寸物体自然对流换热比大空间自然对流换热明显增强这些微尺度条件下的传热

2、现象,其中有些目前还无法由经典的传热理论得到完整的解释。在这样的背景下，分子动力学模拟法应运而生分子动力学方法特征：可预测纳米尺度上的动力学特性。可用于模拟与原子运动路径相关的基本过程。粒子的运动行为是通过经典的Newton运动方程所描述。分子动力学方法是确定性方法，一旦初始构型和速度确定了，分子随时间所产生的运动轨迹也就确定了。原理&模拟方法模拟软件：LAMMPS开源的MD模拟软件可单机运行，也可并行计算可以模拟气体、液体、固体l分子动力学从分子和原子的水平着手l通过研究众多的单个分子运动规律来揭示系统整体的宏观行为l能够在认识热现象的本质的基础上求得不少热现象的宏观参数基本原理：分子原子视

3、为经典粒子统计力学计算机模拟分子动力学模拟的主要步骤有：l 研究系统及其边界以及系统内分子间作用势模型的选取l 模拟算法的建立，实质是求解一组运动方程l 分子的初始位置和动量的设定l 当体系达到平衡后，依据统计公式，获得各宏观参数和输运性质系统达到平衡状态宏观参数不随时间变化相关统计公式所需物理量，例如热力学性质或参数能解决的问题 模拟分子系统得出微观系统的热力学参数。研究热扰动的传播的机理与特点,能够从微观上阐述其物理本质。研究汽液界面之过渡区内热力学参数的变化规律,并用实验来验证和解释微尺度相变的一些特殊规律。应用举例 1 MD模拟纳米尺度结构表面过冷沸腾传热机理 1.控制方程和选择的模型

4、 基本控制方程根据牛顿第二定律：U函数：是i原子的坐标本例子中，水和铜板传热过程将使用MD模拟，因此需要铜，氧和氢原子的参与。CuO,CuH,and CuCu,LennardJones(LJ)的相互作用被应用于：2.计算模型和方法l水分子是基本立方（BC）单元l在铜固体区域中，铜原子设置面心立方单元（FCC）样式每个原子都符合下面的方程：计算域和边界可大致分为四个区域，即水蒸汽，水液体，纳米结构的固体铜表面和恒温器边界区域。如下图：计算域和边界图横截面尺寸图热量条件和不同情况的模型表格在此模拟中使用的参数结果与分析：左图呈现了P-T0模型沸腾过程，时间分别为 0，0.8，1.6和2.4纳秒。右

5、图显示出了水/蒸汽的水温与时间的相互关系。不同模型的过冷沸腾过程Comparison of sub-cooled boiling process in cases P-T2,C-T2 and F-T2 at t=200,300 and 400 ps.(White ellipses indicate where nano-scale gas cavities form).应用举例 2 非平衡分子动力学模拟作为一种计算热力学系数的方法l 假设温度梯度往Z方向形成，给出平衡热力学关系式：NEMD模拟系统比热：热膨胀系数：其中，T是温度，是密度，h表示比焓其中，r表示分子间距离，表示截止距离crL J

6、函数分子通过在偏移力（SF）形式的LJ势相互作用这里执行了两个NEMD模拟，模拟A和模拟B：l 运行A模拟旨在探讨在LE-NEMD的基本特征l 试验B目的是展示在一个高度不均匀系统LE-NEMD的性能计算Cp和p的波动公式：在EMD模拟中，平均值的统计误差通过下面的标准偏差来估计：The specific enthalpy,h,in a local slab of volume V was calculated by比焓和质量密度与温度的关系图模拟A：超临界流体系统左下图示为A模拟得到：剖面ZZ压力分量p ZZ，温度T，比焓h，和密度，圆圈表示由EMD算出的平衡值 温度T，比焓h，和密度的 部分梯度值 黄带表示近似统计误差 圆圈表示由EMD计算出的相应的平衡值MD单元固体结构的快照The green circle at the center in the top view shows the molecule at the origin模拟B：固液共存系统 红线显示了由NEMD模拟得到的局部值 液体和固体状态下的平衡值由圆圈和三角形分别示出n 在这两个模拟中，NEMD运行的局部密度和比焓h等于所述统计误差范围内的平衡值