《计算材料学》第8章势函数.pptx

分子动力学——势函数摘要一．内容回顾二．本章介绍1. 势函数的引入2. 四种势函数模型简介3. 实例: 模拟Au晶体的熔化三．本章小结计算材料学的层次和方法纳观微观介观宏观时间 / s 长度 / m10-1510-1010-610-710-310-4100100缺陷动力学 相变动力学 连续动力学有限元 有限差分分子动力学 蒙特卡罗方法量子力学 密度泛函理论电子原子分子动力学(MD)多粒子体系按照其内部的动力学规律跟踪系统中每个粒子的个体运动，然后根据统计物理规律，得到微观量与宏观量之间关系MD模拟五个要素：? 原子间相互作用势? 初始条件? 边界条件? 积分算法? 原子运动轨迹势函数在MD中的重要性初始条件、边界条件 势函数V(R)演化的 驱动力F = -?V(??)??` =??0+a???R` =R+??`???计算原子受力计算原子加速度和速度a= F /m更新原子坐标势函数的合理性决定了 模拟结果的可靠性统计热力学信息摘要一．内容回顾二．本章介绍1. 势函数的引入2. 三种势函数模型简介3. 实例: 模拟Au晶体的熔化三．本章小结物质的形态与原子间结合力? 世界的精彩离不开原子间结合力……? 结合强度：固体液体气体? 结合力是物质的重要属性材料性能与原子间结合力强度韧性缺陷表面能原子结合力熔点热传导……原子结合的本质-复杂的电子相互作用V(R,r)复杂电子相互作用≠1V= Kx22原子弹簧模型量子力学经典势场：??????增长计 算 量“指数墙”自洽求解精确，但计算量大极限：~上千个原子原子数MD中势函数V(R)? 复杂相互作用≈简洁的函数？? 理论上可行！绝热近似电子和原子核的运动分开考虑想象：电子高速运动形成电子云，原子核近似看作在电子的等效势场V(r)中运动10LAMMPS中的势函数势函数的特点? 函数形式不确定性经验势? 电荷密度相关，猜测合理的形式? 拟合实验值或量子力学结果，获得势参数摘要一．内容回顾二．本章介绍1. 势函数的引入2. 四种势函数模型简介3. 实例: 模拟Au晶体的熔化三．本章小结1. Lennard-Jones PotentialThe original formula:V(r) = 4? [(?/r)12 - [(?/r)6]units of energy units of length?Analytic pair interaction potential?Extensively used for many other materials modeling? useful for qualitative atomistic modeling分子间作用：Lennard-Jones势126?????????? ??= 4???Ar 电荷密度描述分子间弱相互作用，特点：? 二体势? 分子（原子）电子满壳层? 只和原子间距相关，无方向性适用于：中性原子/分子间相互作用: 如稀有气体、液体分子15Renormalized LJ PotentialShortcoming of LJ Potential对于很多真实材料，事实并非如此！Limitations of pair potentialsLimitations of pair potentials Three classes of potentials ionic bonding:Shell modelcovalent bonding: Bond-order potentialsmetallic bonding:Embedded-Atom method21离子键-Shell model?带电粒子之间库伦相互 作用?极化离子之间偶极相互 作用适用于：纯的离子键体系共价键：Tersoff势简写为：三体作用109.5?键长对势 能的贡献键角对势 能的贡献金刚石电荷密度描述共价键相互作用特点：? 三体势? 有显著方向性适用于：非金属半导体及其化 学物: 如C、Si、BN、GaN等金属键：Embedded Atom Method(EAM)势+=idia????12??? ij ijF (? ) ?? (r )E ?Cu电荷密度ii?j (?i )i描述电子密度背景下原子间 的相互作用，特点：? 多体势? 电子密度适用于：金属及合金原子间相互作用: 如Fe、Cu、PtAu等摘要一．内容回顾二．本章介绍1. 势函数的引入2. 四种势函数模型简介3. 实例: 模拟Au晶体的熔化三．本章小结晶体的熔化熔化：晶体从固态转变为液态的过程，对应的温度为熔点计算熔点：通过某一物理量的变化标识固液相变过程 标识方法：1. 体积变化2. 均方位移(mean squared distance, MSD)其中 是对组内的所有原子进行平均 标识原子偏离其平衡位置的程度手心