




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、.第二讲第二讲 分子动力学模拟分子动力学模拟.A. 基础知识 晶体材料 组成材料的原子排列成三维周期性重复的花样。金属和许多非金属固体都是晶体材料。空间点阵单胞每个点有相同的空间环境.A. 基础知识 晶面OA,OB,OC为点阵参数任何一晶面,可用其它与三个主轴的截距OA,OB,OC来表征,通常取这些截距与对应单胞尺寸的倒数。表示为 (OA/OA, OB/OB, OC/OC)然后简化成具有相同比例的最小整数。左图中,ABC, ABE, CEA, CEB的晶面指数(miller指数)为(1,1,1), (1,1,1), (1,1,1), (1,1,1)。这些晶面具有相同晶体学类型形成晶面族1,1,
2、1.对于晶向LM,可用通过原点且平行于LM的直线来表示。对于左图,为OE。表示为OE在三个主轴上的分量与单胞尺寸的最小整数比。表示为:OA/OA, OB/OB, OC/OC=1 1 1A. 基础知识 晶向左图中,CG,AF,DB和EO的方向分别为111, 111, 111, 111。它们属于同一类晶体学类型的方向。形成方向族,表示如下:晶面也可以通过晶面的法线晶向来表示。?晶面()圆滑;晶向棱角 .A. 基础知识 晶体结构a=2r,r为原子半径简单立方晶格.A. 基础知识 晶体结构1,1,2面上的堆垛顺序其面间距为体心立方晶格(BCC)6/a钼(Mo)、钨(W)、钒(V)、-铁(912)等 配
3、位数:配位数为晶格中与任一个原子相距最近且距离相等的原子数目。配位数越大, 原子排列紧密程度就越大。体心立方晶格的配位数为8。 .A. 基础知识 晶体结构面心立方晶格(FCC)具有这种晶格的金属有铝(Al)、铜(Cu)、镍(Ni)、金(Au)银(Ag)、- 铁( -Fe, 9121394)等 配位数:12100与110面原子的堆垛顺序是ABABAB而111面原子排列最紧最密,是塑性滑移最易启动面。其堆垛顺序为ABCABC.A. 基础知识 晶体结构密排六方晶格(HCP)具有这种晶格的金属有镁(Mg)、镉(Cd)、锌(Zn)、铍(Be)等 ra2配位数:12.B.分子动力学模拟简介 分子动力学模拟
4、是一种用来计算一个经典多体系统的平衡和传递性质的方法。 经典这个词意味着组成粒子的核心运动遵守经典力学定律。(忽略量子效应). 不仅可以得到原子的运动轨迹,还可以观察到原子运动过程中各种微观细节。它是对理论计算和实验的有力补充。广泛应用于力学、材料科学和生物物理等。 在分子动力学模拟中,我们一般采用经验势来代替原子间的相互作用势,如Lennard-Jones势 、Mores势、EAM原子嵌入势、F-S多体势。然而采用经验势必然丢失了局域电子结构之间存在的强相关作用信息,即不能得到原子动力学过程中的电子性质。 .C.分子动力学模拟控制方程U为系统的总势能求解方程基本思想是将时间的连续计量离散化为
5、t的时间步长(实际上是微小的时间间隔)。每个粒子在时刻t所受力等于其余所有粒子对它的作用力的矢量和,根据粒子的受力可以计算出时刻t的加速度a(t)。假设在t时间间隔内粒子受力保持恒定,从而得到t+t时刻的位置r(t+t)和速度v(t+t),根据各粒子的新位置又可以计算新的受力进而得到t+2t时刻的位置和速度,如此往复可以得到一段时间内各粒子的运动轨迹,以及整个微结构的演化,通过统计平均的方法进而得到诸多物理力学性质随时间的变化。关键:原子的作用力、运动方程积分算法、边界条件处理.D.原子作用势 分子动力学计算结果的精度关键取决于原子作用力的精确性。 而相互作用力的精确性又取决于原子作用势的好坏
6、与选择。 原子间势的发展经历了从对势到多体势的发展过程。.对势 对势认为原子间相互作用是两两间相互作用,与其他原子的位置无关。 常见的对势包括Alder和Wainwright使用的硬球模型、Lennard-Jones势、Morse势等。L-J势L-J势是针对惰性气体之间相互作用而建立的。它表达的作用力较弱,描述材料的行为比较柔韧。也有人用它来描述铬、钼、钨等体心立方过渡金属。.对势Morse势Morse势可以描述金属,如Cu。与之类似的对势还有Johnson势,常用于描述afa-Fe。. 实际上,在多原子体系中,一个原子位置不同,将影响空间一定范围内的电子云的分布,从而影响其它原子间的有效相互
7、作用。尤其在固体中,这种影响是非常强的。 因此,研究纳米固体力学时多采用多体势。其基本形式一般如下:多体势第i个原子的空间矢量N=3体势能函数系统处于重力场、静电场中外力场的势能。常忽略。N=2体势能函数,两原子之间相互作用。常忽略. 各向同性势函数常被用来描述金属晶体的性质。 这种势函数一般由内聚能和减聚能两项组成 内聚能主要是由该原子所在格点处的局域电子密度决定,局域电子密度来自于目标原子格点的近邻原子的贡献。 减聚能则用来反映原子间的静电排斥作用的贡献,一般以对势的形式表示多体势各向同性多体势原子i处的局域电子密度相互吸引的内聚能相互排斥的减聚能,为各向同性对势。一般通过拟合实验数据得到
8、. 为了解决各种材料的具体问题,目前提出了多种方法来构造具体的各向同性多体势函数。 由 D a w 和 B a s k e s 提 出 的 嵌 入 原 子 方 法(Embedded Atom Method,EAM)便是其中影响较大,使用较广的一种方法多体势各向同性多体势EAM嵌入原子势这种方法将各原子埋入局域电子密度为i的电子云中,并由嵌入原子的能量导出嵌入函数F(i),其中,i由近邻原子的球对称(电子)电荷密度决定。如果进一步考虑电子云的非球形对称分布(共价键),得到修正的原子嵌入法(MEAM)。. 铜晶体 Masao Doyama构造的EAM势函数多体势各向同性多体势EAM嵌入原子势原子J
9、对I的电子密度的贡献函数的截断半径. 键级势函数模型最早由Tersoff 提出,这种势函数根据量子力学中键级(bond order)的概念来描述共价键系统(金刚石,Si,Si-C,石墨,无定形碳,碳纳米管)。多体势键级多体势原子对i和j之间的键级近邻原子间的排斥力近邻原子间的吸引力键级b大小是化学键强度的度量,它与配位数Z的平方根成反比。也就是临近的原子越多,键越弱。截断函数表面上看,Tersoff势是两体势,实际上键级函数bij计及了与i原子存在共价键作用的其它所有原子的信息,因此Tersoff势是一种多体势。.多体势键级多体势Tersoff给出的键级函数bij如下截断函数如下. 它描述分子体系中原子的拓扑结构和运动行为。 包括所有原子的类型表,和一套势能函数表达式及相应的力常数表。 Mayo提出的分子力场DREIDING,能够描述大量有机物、生物分子和所有主族无机分子的普适力场。 DREIDING力场将分子体系的总势能表达为键合能和非键合能之和分子力场.DREIDING力场键合能非键合能Mayo et al. The Journal of Physical Chemistr
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 视听制作服务协议
- 澄清剂采购合同
- 粮油产品品牌故事与情感营销考核试卷
- 矿山带电工作与电气安全管理考核试卷
- 纸品抗菌涂层研究考核试卷
- 红外热成像仪在公共安全中的应用考核试卷
- 毛皮服装设计与时尚界大师访谈考核试卷
- 信用合作社法人治理优化考核试卷
- 肥料制造工程与自动化生产技术研究考核试卷
- 五金行业人才培养与产业转型升级研究考核试卷
- 《化工工艺概论》解析
- 山泉水水厂可行性方案
- 医疗器械经营安全培训必备知识
- 化工公司bluesign认证资料准备清单20201201
- 建设用地报批服务投标方案(技术方案)
- 短引线保护引出线保护以及T区保护
- 完美公司玛丽艳美的观念
- 浙摄影版(2020)信息技术三年级上册第一课认识计算机(课件)
- 校园安全常识测试题卷
- 建筑用玻璃ccc标准
- 第一课中国人民站起来了
评论
0/150
提交评论