Material Studio建模知识讲解.docx

《Material Studio建模知识讲解.docx》由会员分享，可在线阅读，更多相关《Material Studio建模知识讲解.docx（25页珍藏版）》请在三一办公上搜索。

1、Mate rial Studio 建模铁基块体非晶合金-纳米晶转变的动力学模拟过程Discover 模块1原子力场的分配在使用Discover模块建立基于力场的计算中，涉及几个步骤。主要有：选 择力场、指定原子类型、计算或指定电荷、选择non-bond cutoffs。在这些步骤中，指定原子类型和计算电荷一般是自动执行的。然而，在某 些情形下需要手动指定原子类型。原子定型使用预定义的规则对结构中的每个 原子指定原子类型。在为特定的系统确定能量和力时，定型原子使工作者能使 用正确的力场参数。通常，原子定型由Discover使用定型引擎的基本规则来自 动执行，所以不需要手动原子定型。然而，在特殊情

2、形下，人们不得不手动的定型原子，以确保它们被正确地设置。图3-11）计算并显示原子类型：点击EditAtom Selection，如图所示Build 工sF 旦品 Zoom Inttrl +Z令十丫Efe 3py 危也m 牛咖# Intert From.W+XCtrl+CSfilKtCtrl4 AEdjt SetsFind PMterris弹出对话框，如图所示nil Atom SelectionSelection modeCreate a rew selectiori from mil visible atoms厂Select from the ewisbng selectionL Add t

3、o the existing selectionS elect I Help从右边的的元素周期表中选择Fe,再点Select，此时所建晶胞中所有Fe原子都将被选中，原子被红色线圈住即表示原子被选中。再编辑集合，点击Edit-Edit Sets，如图所示Edfel Miwy跚Id 工p# Undo Zoom In Ctrl+Z电 LopyCtrli-C龟由盘Ctr I -h VX Delete己| Jpsfrrt Frann_s.S-elact 1:1Cirl+A轴lEimEdh SetsFind PatiE!5弹出对话框见图，点击New.,给原子集合设定一个名字。这里设置为Fe,则3D视图中会

4、显示“ Fe ”字样，再分配力场： 在工具栏上点击Discover按钮从下拉列表中选择Setup,显示Discover Setup对话框，选择Typing选项卡。旬 Edit Set|g| | Discover GetupEnergyNon-BondJob ControlAutomationForcefield typing and chargingCalculate forcefield types:Calculate partial charges:Forcefield typesNameElementDescriptionaOLmAgsilver, metalagLOAgSilver,

5、+1 oxideal_mAl日luoninium, metalal3oAlalurminium 3+ cwidesT* 1rrr _JkA List all forcefield typesAssign图 3-2 Discover Setup 对话框 Typing 选项卡在Forcefield types里选择相应原子力场，再点Assign(分配)按钮进 行原子力场分配。注意原子力场中的价态要与Properties Project里的原子价 态(Formalcharge)一致。2力场的选择1) Energy力场的选择：力场是经典模拟计算的核心，因为它代表着结构中每种类型的原子与围绕着它 的原子

6、是如何相互作用的。对系统中的每个原子，力场类型都被指定了，它描 述了原子的局部环境。力场包括描述属性的不同的信息，如平衡键长度和力场 类型对之间的电子相互作用。常见力场有COMPASS、CVFF和PCFF。Select下拉菜单中有三个选项： COMPASS力场：COMPASS力场是第一个把以往分别处理的有机分子体系 的力场与无机分子体系的力场统一的分子力场。COMPASS力场能够模拟小分 子与高分子，一些金属离子、金属氧化物与金属。在处理有机与无机体系时， 采用分类别处理的方式，不同的体系采用不同的模型，即使对于两类体系的混 合，仍然能够采用合理的模型描述。 CVFF力场：CVFF力场全名为一

7、致性价力场(consistant valence force field)，最 初以生化分子为主，适应于计算氨基酸、水及含各种官能团的分子体系。其 后，经过不断的强化，CVFF力场可适用于计算多肽、蛋白质与大量的有机分 子。此力场以计算系统的结构与结合能最为准确，亦可提供合理的构型能与振 动频率。PCFF力场：PCFF为一致性力场，增加一些金属元素的力参数，可以模拟含有相应原子的分子体系，其参数的确定除大量的实验数据外，还需要大量的量 子力学计算结果。3非键的设置：|,| | Discover SetupHNp I非键作用力包括范 德华力和库伦力。这里将两者都选上，为的是后 期做 minimi

8、zer优化原子位置时精确度更高，因为考虑了作用力因素多，即两者都考 虑了。Summation method （模拟方法）： Atom Basedatom based基于原子的总量，包括一个原子的截断距离，一个原子的缓冲宽度距离；为直接计算法，即直接计算原子对之间的非键相互作用，当原子对 超出一定距离（截断半径cutoff distance）时，即认为原子对之间相互作用为零（注：cutoff distance指范德瓦尔斯作用力和库仑力的范围，比如：设定截断半 径为5,则表示已分子或原子中心为圆心，以5为半径作圆，半径以外的作用 力都不考虑）。此方法计算量较小，但是可能导致能量和其导数的不连续性。

9、 当原子对间距离在Cut Off半径附近变化时，由于前一步考虑了原子对之间的相 互作用，而后一步不考虑，由此会导致能量发生跳跃。当然，对于较小的体 系，则可以设置足够大的Cutoff半径来保证所有的相互作用都被考虑进来。 Group Basedgroup based基于电子群的，总量中包括一个原子的截断距离，一个原子的 缓冲宽度距离；大多数的分子力场都包括了每个原子之间点电荷的库仑相互作 用。甚至在电中性的物种中也存在点电荷，例如水分子。点电荷实际上反映了 分子中不同原子的电负性。在模拟中，点电荷一般是通过电荷平衡法(charge equilibrium)评价或者力场定义的电荷来分配的。当评价

10、点电荷时，一定要小 心不要在使用Cutoff技术时引入错误的单极项。要了解到这一点，可以参看如 下事实：两个单极，当只有1e.u.电荷时，在10A的位置上其相互作用大约为 33Kcal ；而对于由单位单极分离1A所形成的两个偶极，相同距离其相互作用能 不超过 0.3Kcal/mol。很明显，忽略单极-单极相互作用会导致错误的结果，而忽略偶极-偶极相 互作用则是适度的近似。然而，如果单极相互作用处理不清的话，仍然会出问 题。当non-bond Cutoff使用基于原子-原子基组时，就可能发生，会人为将偶极 劈裂为两个“假”的单极(当一个偶极原子在Cutoff内，另一个在其外)。这就 不是忽略了相

11、对较小的偶极-偶极相互作用，而是人为引入了作用较大的单极- 单极相互作用。为了避免这种人为现象，Materials Studio引入了在Charge Groups 之上的 Cutoff。一个“Charge Group”是一个小的原子基团，其原子彼此接近，净电荷为0或 者接近于0。在实际应用中，Charge Group一般是常见的化学官能团，例如羰 基、甲基或者羧酸基团的净电荷接近于中性Charge Group。Charge Group之间 的距离为一个官能团中心到另一个官能团中心的距离R，Cutoff设置与Atom Based相类似。 Ewald SummationEwald是在周期性系统内计

12、算Non-bond的一种技术。Ewald是计算长程静电 相互作用能的一种算法。Ewald加和方法比较合适于结晶固体。原因在于无限 的晶格内，Cutoff方法会产生较大的误差。然而，此方法放也可以用于无定形 固体和溶液体系。Ewald计算量较大，为o(NM/2)，体系较大时，会占用较多 的内存并花费较长的时间【分子模拟一从算法到应用】。 cell multipole cell based只能用于基于指定数量层。一般情况下，基于Atom适合于孤立体系，对于周期性体系计算量较小， 但是准确性较差；基于Group适合于周期性和非周期性体系，计算的准确性好 一些，计算量最小；Ewald适合于周期性能体系

13、，计算最为准确，但计算量最 大。hM Atom Bas-ed CutoffCutoff dislance:冒虱ASpline width;p-00ABuffer width:f0?50&kOOI- Applv long-rarige correctionOK I Cancel I KelpCutoff distance （截断距离）：指的是范德瓦尔斯作用力和库仑力的范围Spline width :Buffer width :缓冲宽度距离。Setup其他选项保留默认设置即可。4结构优化_在工具栏上点击Discover按钮I ，然后选择Minimizer。或者从菜单栏选 择 Modules | D

14、iscover | Minimizer。显示Discover Minimizer对话框，可以进行几何结构优化计算。注：优化前 （Min），先查看所有原子是否都已分配力场，如果没有，可以手动添加，在 Properties Explorer中双击Forcefield type，然后修改力场类型即可。其次 在Min之前，需要把晶体结构所有原子重新固定。minimizer只是对结构进行优 化，以达到能量最小化。在作动力学（dynamics）之前最好minimizer一下，因 为如果不minimizer计算收敛时间会比较长，能量波动会比较大，而且计算有可 能出错。|51| Discover Minimi

15、zation丞More.Help优化方法Mathod：最陡下降法（Steepest Descent）、共轭梯度法（ConjugateGradient）、牛顿方法（Newton）和综合法（Smart Minimizer）。Convergence level：收敛精度水平。Maximum iteration: 最大迭代数。Optimize cell选中的话表示优化晶胞参数和原子位置。MS Discover结构优化原理分子的势能一般为键合(键长、键角、二面角、扭转角等)和非键合相互 作用(静电作用、范德华作用等)能量项的加和，总势能是各类势能之和，如 下式：总势能二范德华非键结势能+键伸缩势能+键角

16、弯曲势能+双面角扭曲势能+离平面振动势能+库伦静电势能+除了一些简单的分子以外，大多数的势能是分子中一些复杂形势的势能的 组合。势能为分子中原子坐标的函数，由原子不同的坐标所得到的势能构成势 能面(Potential Energy Surface， PES)。势能越低，构象越稳定，在系统中 出现的机率越大；反之，势能越高，构象越不稳定，在系统中出现的机率越 小。通常势能面可得到许多极小值的位置，其中对应于最低能量的点称为全局 最小值(Global Energy Minimum)，相当于分子最稳定的构象。由势能面求最 低极小值的过程称为能量最小化(Energy Minimum)，其所对应的结构为

17、最优 化结构(Optimized Structure)，能量最小化过程，亦是结构优化的过程。通过最小化算法进行结构优化时，应避免陷入局部最小值(local minimum)，也就是避免仅得到某一构象附近的相对稳定的构象，而力求得到全 局最小值，即实现全局优化。分子力学的最小化算法能较快进行能量优化，但 它的局限性在于易陷入局部势阱，求得的往往是局部最小值，而要寻求全局最 小值只能采用系统搜寻法或分子动力学法。在Materials Studio的Discover 模块中，能量最小化算法有以下四种：1) 最陡下降法(Steepest Descent)，为一经典的方法，通过迭代求导， 对多变量的非线

18、性目标函数极小化，按能量梯度相反的方向对坐标添加一位 移，即能量函数的负梯度方向是目标函数最陡下降的方向，所以称为最陡下降 法。此法计算简单，速度快，但在极小值附近收敛性不够好，造成移动方向正 交。最陡下降法适用于优化的最初阶段。2) 共轭梯度法(Conjugate Gradient)，在求导时，目标函数下降方向不 是仅选取最陡下降法所采用的能量函数的负梯度方向，而是选取两个共轭梯度 方向，即前次迭代时的能量函数负梯度方向与当前迭代时的能量函数负梯度方 向的线性组合。此法收敛性较好，但对分子起始结构要求较高，因此常与最陡 下降法联合使用，先用最陡下降法优化，再用共轭梯度法优化至收敛。3）牛顿方

19、法（Newton），以二阶导数方法求得极小值。此法的收敛很迅 速，也常与最陡下降法联合使用。4）综合法（Smart Minimizer），该方法可以混合最陡下降法，共轭梯度 法和牛顿法进行结构优化，在MS中是可选择的。Smart Minimizer中，牛顿法可以设定最大的原子数，如果体系的原子数大于 所设定的值，则计算是会自动地转为前面设定的收敛法（共轭梯度法或最陡下 降法），收敛精度会改为共轭梯度法的默认收敛精度值。点开各种方法后面的More，可设定收敛精度（Convergence），算法（Algorithm）和一维搜索（Line search，指每一次迭代中的精度）等。r_旬 Smart

20、MinimizerConvergerice:|l 000Line search: |0.5kcsl/rftol/?OK. I匚吕 ncel Helpal Smart Minimizer - Conjugate Gradi.当Job结束后，结果被返回到Disco Min目录，最小化的结构被命名为3DAtomistic.xsd，并被保存在 “3D Atomistic Disco Min”目录。还生成一个名为 “3DAtomistic.out ”的文本文档，它包含了有关计算的所有能量信息。同时还生成 “SimulationEnergies.xcd”，它显示了能量随迭代次数的变化情况。如图所示。本次模

21、拟得到如图所示的结构，5高温弛豫打开discover下的Dynamics，如图所示|5il Discove-r Mole-cular DynamicsDynamics | Stress)Ensemble:T emperature:Pressure:Number of steps:Time step:DynamicsTrajectorjJFrame output ev叫RestartMore.rhernnostat.Ensemble （系综）：NVE、NVT、NPT、NPH。Temperture ：目标温度。Pressure :给系统所施加的压力。Number of steps :整个动力学所运

22、行的总步数。Timp step :每一动力学步骤所花费的时间。Dynamics time : Number of stepsxTimp step。Trajectory Save : Coordinates表示保存坐标；Final Structure表示只保存最终结构；Full表示保存所有。Frame output every :若输入5000，则表示每5000步输出一次体系构型文件。此操作表示结构在2000K的温度下进行弛豫，此过程原子通过迁移、运动或者扩散，逐步降低原来的高内能态，向稳定的低内能态转变。运行结束后，可以通过调用Animation观看三维动画，见图动画工具条可以控制三维窗口中动

23、画文件的显示。它包含以下命令:Play Backwards :倒映动画文件。Step Backwards :每次向后放一帧Stop：停止放映。Step Forwards :每次一帧加速放映。Play：放映动画。Pause：暂停放映，再按一次后继续放映。Animation Mode :显示动画模式下拉菜单，6.1系综简介系综（ensemble）是指具有相同条件系统（system）的集合。平衡态的分子动力学 模拟，总是在一定的系综下进行。系综是统计力学中非常重要的概念，系统的 一切统计特性基本都是以系综为起点推导得到的。实际应用时，要注意选择适 当的系综，如（N，T，P）常用于研究材质的相变化等。

24、1）在微正则系综（micrononical ensemble）中，模型体系的粒子数N、体积V 及内能（热力学能）E （在热力学通常用U表示内能）。孤立、保守的系统。值 得注意的是：体系总能量，即势能和动能的总和，是保持守恒的，常被用来判 断积分的精度固定不变。它对应于绝热过程，即体系与环境没有热交换，不存 在温度T和压力p的控制因素。由于体系的能量E是守恒的，体系的动能和势 能之间互转化。一般说，给定能量的精确初始条件是无法得到的。能量的调整 通过对速度的标度进行，这种标度可能使系统失去平衡，迭代弛豫达到平衡。2）NVT系综（正则系综）正则系综（canonical ensemble）中，体系的

25、粒子数N、体积V及温度T保持 不变，且总动量保持不变。因此正则系综动力学有时也被称为恒温动力学。为 了控制体系的温度，就需要设置一个“虚拟”的热浴环境，与体系进行能量交 换。常用的热浴（bath）包括：Nose-Hoover， Berendsen， Andersen 以及 “velocity scaling （速度标定）”方法等。3）NPT系综（恒温恒压系综）恒温恒压系综中，体系的粒子数N、压力p、温度T都是恒定不变的。恒温恒压系 综允许体系的“体积”发生变化。这里的体积的变化有两种方式，一种是只变化 尺寸而保持形状（比如对于晶体来说，晶格类型维持不变，但是晶胞参数中的 a，b，c可以变化），

26、另一种是同时变化形状和尺寸（即晶格类型和晶胞参数都可以变化）。压强P与体积共轭控压可以通过标度系统的体积来实现。目前有 许多调压的方法都是采用的这个原理。4）NPH系综（恒焓恒压系综）NPH系综中体系的粒子数N、压力p及体系的焓H（H=E+pV）是守恒的，例如节流 膨胀就是一恒焓过程。在模拟中较少见。6.2系综控温机制系综的控温：温度调控机制可以使系统的温度维持在给定值，也可以根据外界 环境的温度使系统温度发生涨落。一个合理的温控机制能够产生正确的统计系 综，即调温后各粒子位形发生的概率可以满足统计力学法则。系综控温机制主要有 :Velocity Scale、 Nose、 Berendsen。

27、T hermostat:Il Andersen3T emp. difJ&rence|10.0K.Q ratioICollision ratioI1-0Dec即 constant|Cancel |HelpThermostatThermostat下拉菜单有四个： 1）Velocity Scale （直接速度标定法）：系统温度和粒子的速度直接相关，可 以通过调整粒子的速度使系统温度维持在目标值。实际分子动力学模拟中，并 不需要对每一步的速度都进行标定，而是每隔一定的积分步，对速度进行周期 性的标定，从而使系统温度在目标值附近小幅波动。直接速度标定法的优点是 原理简单，易于程序编制。缺点是模拟系统无法

28、和任何一个统计力学的系综对 应起来；突然的速度标定引起体系能量的突然改变，致使模拟系统和真实结构 的平衡态相差较远。2) Nose :该方法可以把任何数量的原子与一个热浴耦合起来，可以消除局域的 相关运动，而且可以模拟宏观系统的温度涨落现象。Andersen:3) Berendsen控温机制：又称Berendsen外部热浴法。其基本思想是假设系统和 一个恒温的外部热浴耦合在一起，通过热浴吸收和释放能量来调节系统的温 度，使之与恒温热浴保持一致。对速度每一步进行标定，以保持温度的变化率 与热浴和系统的温差(Tbath-T(t)成比例。6.3系综空压机制1OTBarostatB arostdl:B

29、erendsen ParrinellaCell ndss:And&rseriarnuBerendsenLL 1 K.CancelHelr下拉菜单有3项：Andersen :假定系统与外界“活塞”耦合，当外部压强不能补偿系统内部压强时，“活塞”运动引起系统均匀地膨胀或收缩，最终使得系统压强等于外部压强。Andersen方法具有重要的意义，后来的各种压力控制方法基本都是基于Andersen思想发展起来的。Berendsen :这种方法是假想把系统与一“压浴”相耦合。Parrinello :这种方法允许原胞的形状与体积同时发生变化，以达到与外压平 衡。这种方法是对Anderson调压方法的一种扩展，

30、可以实现对原胞施加拉伸剪 切以及混合加载情况的模拟，因此在对材料的力学性质的分子动力学模拟中， 得到了广泛地应用。运行结束后可以2 Forcite 模块2.1 Quench （快冷）在工具栏上点击按钮，选择calculation，弹出对话框，如图所示1.选择Quench （快冷，淬火），再点击More出现如图所示对话框:旬 Fore it e Quench DynamicsQuench ever |5000stepsDriarriits options:More.GeomBtrAJ optimization options:阴陷.|Help再点击Dynamics options的more出现如

31、图所示：Inilial velocities:|RandomRandomTemperature:扇航检Initial velocities :第一次由于设置速度，所以只能选择Random （随即速度）第二次以及以后运行则可选择Current （当前速度）了，此时速度为上一次结束 的速度。Fo-rcitff DynamicsDiinamics | T吨rm瑚BarostatIntegratior toleranceEnergy deviation: 150000.0 kcd/molT raiectaryW Include velocities际 Includ已 forcesHalj注意：模拟退火

32、的时候要加力。运行结束后会得到一些文件，有1) 3D Atomistic.xtd，这是快冷后得到的结 构，从这里可以得知得到的是非晶合金，见图。2) Status.txt以及3D Atomistic.txt包含了快冷过程的相关参数设置以及结果数据。3)3D AtomisticTemperature.xcd描述了温度与时间的关系，见图。4) 3D Atomistic Energies.xcd 描述了几种能量(势能、动能、非键能以及总能量)最时间的变化关系(见 图)等。2.选择退火（anneal）如图所示点击more出现下图：E3151| Forcite Anneal DynamicsAnneal

33、ing cycles :运行一次退火所作的退火循环次数。Initial temperature :一次退火循环的起始温度也是退火循环的终止温度。Mid-cycle temperature : 一次退火循环包括升温过程和降温过程中的最高温度。Heating ramps per cycle :一次循环中加热过程的温度梯度步数，冷却过程的温 度下降梯度(cooling ramps per cycle)步数与加热过程的温度梯度步数相等。Dynamics steps per ramp :每一温度梯度的动力学步数。Total number of steps : Annealing cycles x(Hea

34、ting ramps per cycle+coolingHeating ramps per cycle)xDynamics steps per ramp (即上图中的总步数=5x10x500)园 For cite DynarriicsIHH Fare its DynamicsFrame output ever:stepsDynamics Thermostat Barostat L.3ncedIntegratior talerance|50000.0Energy deviation:T 旧 tciryInclude velocitiesaIrclude forces：kcsl/nriolHelp

35、 Help目标温度根据快冷得到700K的结构而设定为700K，而两组模拟中，中间最高温度(Mid-cycle temperature)分别设为835K和830K，因其晶化温度大致在825K左右【】。以及得到其他信息，可在Project中点击查看。再对上面的结构作X衍射。3 Reflex 模块如下图所示，调出Powder Diffraction 工具从工具栏选择Reflex匚工具，然后选择Powder Diffraction,或者从菜 单栏选择 Modules | Reflex。里y甲姑a 山2丑5画向皿HI B end s CASTERCCDCCn farmersDRB4Mol3DPPDl5t

36、OTTffrEquilibria-Forcite fiau sranGULP ttnetix MfiwDynMeso citeWorp_holcg/NETEPPDlmorphQMERA也四膻e Pra-ssing 史DtflracFu?ner RefiriementU% Powder QPA l Gwdrrr Cryslcillinlty .L Aowd er olve司 ,口McrEpt.-Jeb IdSflyw Setting j.LA.出虾1吁Indexing I J1 Powd |Hi| Reflex Powder DiffractionTemperature Factors | Af

37、ynmetry | EMpenmental Data | Display DiffractometerRadiation Profiles SampleH angeMinMawStep size2-T held：5.000，45.000,0.0500*Vdhkk0.056628 V?0.4968112?Line shift G eometrjj：| Bregg-Brunton。_jZero point:0.0000 Shihttl:0.0000ShiftflZ |0.0000n Monitor diffractionResetHelp显示 Reflex Powder Diffraction 对

38、话框，如图 91 所示。Powder Diffraction对话框由8个不同的选项卡组成，包括你需要的所有设置。Diffractometer -设置基本的扫描设置，例如2-theta范围和线性变化；Radiation-设置不同的衍射线类型，可以选择X射线、电子和中子射线；Profiles -设置粉末衍射图显示的峰形函数并加宽显示衍射图；Sample -设置样品尺寸；Temperature Factors -包括控制修正原子热振动对衍射图的影响；Asymmetry -控制用于修改峰形的任何不对称性修正；Experimental Data -允许你添加实验数据进行对比；Display -设置常规的显示属性，这对控制图形数据是很重要的。第一步是计算衍射图。