材料设计、计算及模拟.PPTVIP

Schr?dinger 方程 原子核和 电子动能 原子核-电子 相互作用 电子-电子 相互作用 原子核-原子核 相互作用 * Kohn-Sham（沈吕九） 方法 第二项：静电作用能 第三项：无相互作用体系的动能 Exx[ρ]: 含有交换-相关能的项（难点） * 材料电子结构模拟——密度泛函理论（DFT） Bismuth-induced embrittlement of copper grain boundaries Calculated charge density from copper grain-boundary region. G. Duscher, et.al. Nature Materials, 22 August, 2004 * Atomic structures of the SrTiO3/Si(001) interfaces C.J. Forst, Nature 427 53 (2004) * * MD是经典力学方法，针对的最小结构单元不再是电子而是原子 因原子的质量比电子大很多，量子效应不明显，可近似用经典力学方法处理 20 世纪 30 年代, Andrews 最早提出分子力学(MM)的基本思想；40 年代以后得到发展, 并用于有机小分子研究。90年代以来得到迅猛发展和广泛应用 材料原子层次模拟——分子动力学 * 构造出简单体系的势能函数, 简称 势函数 或 力场(force field)。 利用势函数，建立并求解与温度和时间有关的牛顿运动方程，得到一定条件下体系的结构随时间的演化关系。 再将粒子的位置和动量组成的微观状态对时间平均，即可求出体系的压力、能量、粘度等宏观性质以及组成粒子的空间分布等微观结构。 势函数：势能与原子位置的关系。且往往是不知道的需要通过其他方法，如量子化学方法及实验数据获得。 分子动力学——基本原理 * 分子动力学模拟中，忽略了量子效应后，系统中粒子将遵循牛顿运动定律。为了得到原子的运动，可以采用各种有限差分法来求解运动方程。常用的有以下几种算法：Verlet算法、Velocity-verlet算法、Leap—frog算法、Beeman算法、Gear算法、Rahman算法。其中Verlet算法虽然精度比Gear算法稍差，但使用方便，占用存储量少，稳定性好，因此使用较为广泛。 * 粒子系综的控制理论 调压技术 ①Andersen压浴法: Andersen压浴法假想系统与一活塞相接触。其作法为对模拟系统的体积乘以压力标度因子,对原子质心坐标乘以Cp. ②Parrinello—Rahman(P—R)方法: P—R法用于处理晶格的形状和体积都发生变化的情况，可实现对原胞施加拉伸剪切以及混合加载情况的模拟，在材料科学中得到了广泛的应用. * 粒子系综的控制理论 调温技术 ①速度标度法: 速度标度法是保持系统温度恒定最简单的方法。其具体做法是每隔一定的模拟步数，将原子运动的速度乘以修正系数.使体系的动量始终保持不变． ②Nose-Hoover热浴法:Nose-Hoover热浴法假想系统与一个温度为期望值的虚拟热浴相接触。热浴的温度足够大，使所研究的体系的温度随时在热浴中获取和释放． * 分子动力学系综 即表示具有确定的粒子数（N）、体积（V）、温度（T）。总能量（E）和系统压强（P）可能在某一平均值附近起伏变化。正则系综的特征函数是亥姆霍兹自由能 F（N,V,T）。 表示具有确定的粒子数（N）、体积（V）、总能量（E）。系综的温度（T）和系统压强（P）可能在某一平均值附近起伏变化。微正则系综的特征函数是熵 S（N,V,E）。 表示具有确定的粒子数（N）、压强（P）、温度（T）。其总能量(E)和系统体积(V)可能存在起伏。特征函数是吉布斯自由能G(N,P,T)。 表示具有确定的粒子数（N）、压强（P）、焓（H）。这种系综在实际的分子动力学模拟中已经很少遇到了。 微正则系综NVE （micro-canonical ensemble） 正则系综NVT （canonical ensemble） 等温等压NPT (constant-pressure, constant-temperature) 等压等焓NPH (contant-pressure, constant-enthalpy) * 材料原子层次模拟——分子动力学 Surface roughness with adatom incident energy at 300?K after one and fiveML deposition of Al or Ni films on Ni(1?1?1). The black spheres represent Ni adatoms