高等计算力学课中-分子动力学概述

高等计算力学分子动力学模拟主要参考资料：1.高等计算力学杨庆生郑代华2.材料变形与破坏的多尺度分析范镜泓3.TheArtofMolecularDynamicsSimulationD.C.RAPAPORT一、分子动力学概述主要内容二、分子动力学基本理论三、分子动力学模拟方法四、时间积分算法和时间步长五、分子动力学模拟过程及应用六、多尺度模拟方法七、在复合材料中的应用一、分子动力学概述1.蒙特卡罗(MonteCarlo)方法基础2.分子动力学(MolecularDynamics)方法基础3.准连续介质法4.原子结构力学法纳米力学：在纳观尺度上研究材料的力学行为分子动力学是在原子、分子水平上求解多体问题的重要的计算机模拟方法，可以预测纳米尺度上的材料动力学特性。获得微观态的方法不同搭一个盒子（一个初始构象），怎么驱动原子分子在设定的外界条件比如温度，压力下访问一些可能的状态？有2种方法：MD,或MC。MD,动力学，按照牛顿老人家的说法，原子分子也就是些球球，它们在力的作用下会运动，那好，我计算原子的每一步的位置，受力，加速度，位移，让它自动去跑，这叫模仿自然演进。MC，蒙特卡洛，是随机动力学方法，不考虑受力，原子的移动是随机的，至于移动的合理不合理，那由下一步的能量决定。能量低，肯定合理，接受；能量高，高多少？越高越不合理，那好，我根据能量高多少来决定接受的概率。由此可见，MD和MC同属于采样方法，但核心思想有较大差别。它们唯一的共同点是，需要计算每一步的能量。（分子动力学MD中计算力，而原子受力是体系能量对原子坐标的一阶导数。一、分子动力学概述分子动力学的时空尺度二、分子动力学基本理论基本思想：通过原子间相互作用势，求出每一个原子所受到力，在选定的时间步长、边界条件、初始位置和初始速度下，对有限数目的分子（原子）建立其牛顿动力学方程组，用数值方法求解，得到这些原子经典运动轨迹和速度，然后对足够长时间的结果求统计平均，从而得到所需要的宏观物理量和力学量。二、分子动力学基本理论分子动力学基本方程两个假设：所有粒子的运动都遵循经典牛顿运动定律；粒子间的相互作用满足叠加原理。1.Lagrange运动方程2.Hamilton运动方程二、分子动力学基本理论原子间作用势函数1对势模型方程右端第一项描述原子间的排斥作用，第二项描述原子间的吸引作用；σ为零势距离，ε表征原子间的吸引/排斥强度。Lennard-Joans势0.811.21.41.61.822.22.42.6-1-0.500.511.52间间间间Lennard-Jones间间间间62ijr二、分子动力学基本理论2多体势模型镶嵌原子势EmbeddedAtomMethod(EAM)在局域背景电子云密度环境下的结合能，将系统总能量分解为二、分子动力学基本理论截断半径cccrrrrrUrU0)()(忽略截断半径之外的所有作用三、分子动力学模拟方法几何模型的建立面心立方结构六方最密堆积•初始条件–随机初始条件给法之一要求大小：条件一：规则给法条件二：随机给法max||Vv)1random2(maxVv三、分子动力学模拟方法假设的粒子总数不变初始条件和边界条件三、分子动力学模拟方法PBC(PeriodicBoundaryConditions)P(x)=P(x+nL)PBC要求元胞的尺寸必须大于两倍的原子截断半径。三、分子动力学模拟方法温度调控1直接速度标定法任给初始条件，模拟到平衡，得到系统平衡态温度T。一般TTeq。令用速度再模拟直到平衡，若所得温度仍不等于Teq，再进行上述过程TTfeq/1/diifvv三、分子动力学模拟方法2外部热浴法引入一个与虚拟粒子碰撞的随机力想象系统浸在热浴当中系统和热浴间的相互作用强度由随机碰撞的频率决定碰撞的几率等于Nu×dt如果一个粒子经历碰撞，它的速度将从约束温度下的Maxwell分布中随机抽取总能量和总动量均不守恒三、分子动力学模拟方法压力调控1直接压力标定法3Parrinello-Rahman等应力控制法2Anderson等压控制法四、时间积分算法和时间步长时间积分算法Verlet算法粒子位置的Taylor展开式：粒子位置：粒子速度：t(t)t)-(t(t)2t)(t2iarrriiit2t)-(tt)(t(t)iiirrv粒子加速度：3i2iit(t)61t(t)21t(t)(t)t)(tbavrrii3i2iit(t)61t(t)21t(t)(t)t)(tbavrrii(t)(t)imiiFa开始运动时需要r(t-Δt)：t(0)(0)t)(ivrr+缺点：Verlet算法处理速度非常笨拙四、时间积分算法和时间步长时间积分算法蛙跳(Leap-frog)算法：半步算法t(t)t)21-(tt)21(tiavviitt)21(t(t)t)(tivrrii2t)21-(tt)21(t(t)iiivvv1.首先利用当前时刻的加速度，计算半个时间步长后的速度：2.计算下一步长时刻的位置：3.计算当前时刻的速度：t-Δt/2tt+Δt/2t+Δtt+3Δt/2t+2Δtvrv开始运动时需要v(-Δt/2)：t/2(0)(0)t/2)(iavv四、时间积分算法和时间步长时间步长室温下，∆t≈1fs(femtosecond10-15s)，温度越高，∆t应该减小太长的时间步长会造成分子间的激烈碰撞，体系数据溢出;太短的时间步长会降低模拟过程搜索相空间的能力五、分子动力学模拟过程及应用Adsorptionofmethaneonzeolite分子筛是具有均匀的微孔、其孔径与一般分子大小相当为一类吸附剂或薄膜类物质。在化学化工中经常被用作催化剂或者催化剂的载体。五、分子动力学模拟过程及应用Layer-CellModelereasilycreatessimulationofwater-benzeneinterface在体系内部物理性质和化学性质完全均一的一部分称为“相”（Phase)。相与相之间在指定的条件下有明显的界面，在界面上，从宏观的角度看，性质的改变是飞跃式的。五、分子动力学模拟过程及应用Ru-Al合金断裂过程动态模拟C.S.Becquart,D.Kim,J.A,Rifkin,andP.C.Clapp,Mat.Sci.Engin.,A170,87(1993)断裂点周围的损坏区域五、分子动力学模拟过程及应用六、多尺度模拟方法多尺度方法：在物理量迅速变化的局部区域采用分子动力学的原子模拟方法，而在其周围区域可以采用连续介质力学方法或较为粗化的数值模拟技术多尺度模拟中的关键技术是不同尺度之间的连接与过渡，特别是分子动力学模拟区域与连续介质区域的连接。六、多尺度模拟方法多尺度分析：1.跨原子/连续介质（第一类）2.跨连续介质微/细/宏观（第二类）3.时间多尺度分析范镜泓等发展了一种基于建立中间尺度表征单元的有效本构方程来定量地将上一大尺度的物理量与下一小尺度的相关物理量和几何参数连接起来。七、在复合材料中的应用利用分子动力学模拟可以在原子级别上对晶体表面和界面进行优化驰豫，并分析计算其界面原子构型、电子密度分布、界面能量状态等。对于复合材料界面，界面的力学性能是另一个更为关注的问题，分子动力学在获得界面平衡构型后，可以对其进行加载模拟，以研究界面应力和载荷传递，计算界面力学性能参数，模拟界面变形和失效。对凝聚相界面进行原子尺度的模拟是分子动力学(或第一性原理分子动力学)方法的一个主要研究内容，这部分研究涵盖了复合材料界面扩散反应、界面结构、界面力学性能以及复合材料界面失效等各个方面。鉴于其对界面平衡构型及其结合能、复合材料界面应力及载荷传递，复合材料界面失效的模拟研究李健，杨延清等分子动力学模拟在复合材料界面研究中的进展，江浩，岳红等，复合材料界面相互作用的分子模拟研究与进展罗旋，钱革非等，界面对复合材料静态及弯曲力学性能影响的分子动力学模拟华南理工大学杨佳等人通过分子动力学模拟，分析了聚氯乙烯(PVC)与玻璃纤维(GF)间的界面黏合行为，并从微观分子角度探讨GF增强PVC界面能量变化及其增强机理西北工业大学江浩等人通过分子动力学模拟，模拟研究了天然橡胶(NR)／顺丁橡胶(BR)混合胶体的力学性能，从微观上解释了NR和BR共混后性能得到改善的原因，并通过界面结合能和径向分布函数分析揭示了混合物组分之间的相互作用的本质。宋海洋等人用分子动力学方法对碳纳米管/金(CNT/Au)复合材料在拉伸下的力学行为进行了模拟研究。计算结果表明：嵌入方式不同，对复合材料的屈服应力和应变有本质影响；Ni涂层的加入均使得复合材料界面力学性能有所提高。MolecularDynamicsStudyofInterfaceStructureinCompositesComprisingSurface-ModifiedSiO2NanoparticlesandaPolyimideMatrix分子动力学软件软件名称操作系统应用范围DL_PolyNNIX做界面体系LAMMPSLINUX做材料体系GROMACS无限制做蛋白体系AMBERNNIX做生物体系MaterialsExploer无限制化学和材料体系表1分子动力学软件对比