Material Studio建模学习资料.docx

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

2、情形下，人们不得不手动的定型原子，以确保它们被正确地设置。 3-1 图Atom SelectionEdit 点击1）计算并显示原子类型：，如图所示 弹出对话框，如图所示 原子Fe，再点Select，此时所建晶胞中所有从右边的的元素周期表中选择FeEdit 点击Edit都将被选中，原子被红色线圈住即表示原子被选中。再编辑集合， Sets，如图所示 ，则，给原子集合设定一个名字。这里设置为弹出对话框见图，点击New.Fe Fe”字样，再分配力场：视图中会显示“3D Discover ，从下拉列表中选择Setup，显示按钮在工具栏上点击Discover 对话框，选择SetupTyping选项卡。 T

3、yping选项卡图3-2 Discover Setup对话框（分配）按钮进行AssignForcefield 在types里选择相应原子力场，再点里的原子价态Properties Project原子力场分配。注意原子力场中的价态要与 Formalcharge）一致。（ 力场的选择2Energy ）1 力场的选择：因为它代表着结构中每种类型的原子与围绕着它的力场是经典模拟计算的核心，原子是如何相互作用的。对系统中的每个原子，力场类型都被指定了，它描述了原子的局部环境。力场包括描述属性的不同的信息，如平衡键长度和力场类型对之间的电子相互作用。常见力场有 。PCFF、CVFF和COMPASSSele

4、ct下拉菜单中有三个选项： COMPASS 力场：COMPASS 力场是第一个把以往分别处理的有机分子体系的力场与无机分子体系的力场统一的分子力场。COMPASS 力场能够模拟小分子与高分子，一些金属离子、金属氧化物与金属。在处理有机与无机体系时，采用分类别处理的方式，不同的体系采用不同的模型，即使对于两类体系的混合，仍然能够采用合理的模型描述。 CVFF力场：CVFF 力场全名为一致性价力场(consistant valence force field)，最初以生化分子为主，适应于计算氨基酸、水及含各种官能团的分子体系。其后，经过不断的强化，CVFF 力场可适用于计算多肽、蛋白质与大量的有机

5、分子。此力场以计算系统的结构与结合能最为准确，亦可提供合理的构型能与振动频率。 PCFF力场：PCFF为一致性力场，增加一些金属元素的力参数，可以模拟含有相应原子的分子体系，其参数的确定除大量的实验数据外，还需要大量的量子力学计算结果。 3 非键的设置： minimizer为的是后期做非键作用力包括范德华力和库伦力。这里将两者都选上， 优化原子位置时精确度更高，因为考虑了作用力因素多，即两者都考虑了。 （模拟方法）：Summation method Atom Based 一个原子的缓冲宽基于原子的总量，包括一个原子的截断距离，atom based当原子对超出为直接计算法，即直接计算原子对之间的

6、非键相互作用，度距离；（注：即认为原子对之间相互作用为零cutoff distance）时，一定距离（截断半径，指范德瓦尔斯作用力和库仑力的范围，比如：设定截断半径为5cutoff distance。半径以外的作用力都不考虑）以5为半径作圆，则表示已分子或原子中心为圆心，当原子对间距离在但是可能导致能量和其导数的不连续性。此方法计算量较小，半径附近变化时，由于前一步考虑了原子对之间的相互作用，而后一步OffCut 不考虑，由此会导致能量发生跳跃。当然，对于较小的体系，则可以设置足够大的Cutoff半径来保证所有的相互作用都被考虑进来。 Group Based基于电子群的，总量中包括一个原子的截

7、断距离，一个原子的group based大多数的分子力场都包括了每个原子之间点电荷的库仑相互作用。缓冲宽度距离；点电荷实际上反映了分子中例如水分子。甚至在电中性的物种中也存在点电荷，）charge equilibrium在模拟中，点电荷一般是通过电荷平衡法（不同原子的电负性。评价或者力场定义的电荷来分配的。当评价点电荷时，一定要小心不要在使用技术时引入错误的单极项。要了解到这一点，可以参看如下事实：两个单Cutoff；而对于由33Kcal电荷时，在10A的位置上其相互作用大约为极，当只有1e.u. 。相同距离其相互作用能不超过0.3Kcal/mol单位单极分离1A所形成的两个偶极，偶极相互-单

8、极相互作用会导致错误的结果，而忽略偶极很明显，忽略单极作用则是适度的近似。然而，如果单极相互作用处理不清的话，仍然会出问题。原子基组时，就可能发生，会人为将偶极劈裂使用基于原子-当non-bond Cutoff。这就不是Cutoff内，另一个在其外）为两个“假”的单极（当一个偶极原子在单极相偶极相互作用，而是人为引入了作用较大的单极-忽略了相对较小的偶极-之上Charge Groups互作用。为了避免这种人为现象，Materials Studio引入了在 。的Cutoff0”是一个小的原子基团，其原子彼此接近，净电荷为Charge Group一个“一般是常见的化学官能团，例如羰。在实际应用中，

9、Charge Group0或者接近于之间的。Charge Group基、甲基或者羧酸基团的净电荷接近于中性Charge GroupAtom Based设置与R，Cutoff距离为一个官能团中心到另一个官能团中心的距离 相类似。 Ewald Summation是计算长程静电Ewald的一种技术。bond-Non是在周期性系统内计算Ewald 相互作用能的一种算法。Ewald加和方法比较合适于结晶固体。原因在于无限的晶格内，Cutoff方法会产生较大的误差。然而，此方法放也可以用于无定形固体和溶液体系。Ewald计算量较大，为o(N3/2)，体系较大时，会占用较多的内存并花费较长的时间【分子模拟从

10、算法到应用】。 cell multipole cell based只能用于基于指定数量层。 一般情况下，基于Atom适合于孤立体系，对于周期性体系计算量较小，但是准确性较差；基于Group适合于周期性和非周期性体系，计算的准确性好一些，计算量最小；Ewald适合于周期性能体系，计算最为准确，但计算量最大。 ：指的是范德瓦尔斯作用力和库仑力的范围。Cutoff distance（截断距离） ：Spline width ：缓冲宽度距离。Buffer width 其他选项保留默认设置即可。Setup 结构优化4 。或者从菜单栏选Minimizer在工具栏上点击Discover按钮，然后选择 。择Mo

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

12、（优化方法MathodSteepest Descent）、共轭梯度法（ Smart MinimizerNewtonGradient）、牛顿方法（）和综合法（）。 ：收敛精度水平。Convergence levelMaximum iteration：最大迭代数。 Optimize cell选中的话表示优化晶胞参数和原子位置。 MS Discover 结构优化原理 分子的势能一般为键合（键长、键角、二面角、扭转角等）和非键合相互作用（静电作用、范德华作用等）能量项的加和，总势能是各类势能之和，如下式： 总势能 = 范德华非键结势能 + 键伸缩势能 + 键角弯曲势能 + 双面角扭曲势能 + 离平面振

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

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

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

16、迅速，也常与最陡下降法联合使用。 4）综合法（Smart Minimizer），该方法可以混合最陡下降法，共轭梯度法和牛顿法进行结构优化，在MS中是可选择的。 Smart Minimizer中，牛顿法可以设定最大的原子数，如果体系的原子数大于所设定的值，则计算是会自动地转为前面设定的收敛法（共轭梯度法或最陡下降法），收敛精度会改为共轭梯度法的默认收敛精度值。 点开各种方法后面的More，可设定收敛精度（Convergence），算法（Algorithm）和一维搜索（Line search，指每一次迭代中的精度）等。 3D Min目录，最小化的结构被命名为结束后，结果被返回到当JobDisco

17、3D “”目录。还生成一个名为“Atomistic.xsd，并被保存在3D Atomistic Disco Min”的文本文档，它包含了有关计算的所有能量信息。同时还生成Atomistic.out ”，它显示了能量随迭代次数的变化情况。如图所示。“SimulationEnergies.xcd 本次模拟得到如图所示的结构， 高温弛豫5 打开discover下的Dynamics，如图所示 NPH。、NPT、Ensemble（系综）： NVENVT Temperture：目标温度。 ：给系统所施加的压力。Pressure ：整个动力学所运行的总步数。Number of steps ：每一动力学步骤所

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

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

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

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

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

23、度直接相关，可以Velocity Scale（直接速度标定法）1）并不需通过调整粒子的速度使系统温度维持在目标值。实际分子动力学模拟中，对速度进行周期性的标要对每一步的速度都进行标定，而是每隔一定的积分步，直接速度标定法的优点是原理简单，从而使系统温度在目标值附近小幅波动。定，突然缺点是模拟系统无法和任何一个统计力学的系综对应起来；易于程序编制。致使模拟系统和真实结构的平衡态相差较的速度标定引起体系能量的突然改变， 远。：该方法可以把任何数量的原子与一个热浴耦合起来，可以消除局域的2）Nose 相关运动，而且可以模拟宏观系统的温度涨落现象。Andersen: 外部热浴法。其基本思想是假设系统和

24、Berendsen又称:Berendsen3）控温机制 一个恒温的外部热浴耦合在一起，通过热浴吸收和释放能量来调节系统的温度，对速度每一步进行标定，以保持温度的变化率与热 使之与恒温热浴保持一致。浴和系统的温差(Tbath-T(t)成比例。 6.3 系综空压机制 下拉菜单有3项：耦合，当外部压强不能补偿系统内部压强时，“活塞”Andersen：假定系统与外界运动引起系统均匀地膨胀或收缩，最终使得系统压强等于外部压强。活塞”“Andersen方法具有重要的意义，后来的各种压力控制方法基本都是基于Andersen 思想发展起来的。 ：这种方法是假想把系统与一“压浴”相耦合。Berendsen：这种

25、方法允许原胞的形状与体积同时发生变化，以达到与外压平衡。Parrinello调压方法的一种扩展，可以实现对原胞施加拉伸剪切以这种方法是对Anderson得到了及混合加载情况的模拟，因此在对材料的力学性质的分子动力学模拟中， 广泛地应用。 运行结束后可以 模块2 Forcite 2.1 Quench（快冷） ，弹出对话框，如图所示calculation按钮，选择在工具栏上点击 出现如图所示对话框：（快冷，淬火），再点击More选择1.Quench 出现如图所示：再点击Dynamics options的more ，第二次以及以（随即速度）Random：第一次由于设置速度，所以只能选择Initial

26、 velocities后运行则可选择Current（当前速度）了，此时速度为上一次结束的速度。 注意：模拟退火的时候要加力。，这是快冷后得到的结构，从这里可以3D Atomistic.xtd）运行结束后会得到一些文件，有1包含了快冷过程的相关参3D Atomistic.txtStatus.txt以及2得知得到的是非晶合金，见图。）43D Atomistic Temperature.xcd描述了温度与时间的关系，见图。）数设置以及结果数据。3描述了几种能量（势能、动能、非键能以及总能量）最时间的变化3D Atomistic Energies.xcd 关系（见图）等。 选择退火（2.anneal）

27、如图所示 点击more出现下图： Annealing cycles：运行一次退火所作的退火循环次数。 Initial temperature：一次退火循环的起始温度也是退火循环的终止温度。 cycle temperature：一次退火循环包括升温过程和降温过程中的最高温度。Mid-：一次循环中加热过程的温度梯度步数，冷却过程的温度下降梯度per cycleHeating ramps ）步数与加热过程的温度梯度步数相等。（cooling ramps per cycle ：每一温度梯度的动力学步数。Dynamics steps per rampHeating ramps per cycle+coo

28、ling Heating ramps per （Total number of steps：Annealing cycles ） 10500=5cycle）Dynamics steps per ramp（即上图中的总步数 cycle （，而两组模拟中，中间最高温度Mid-700K700K目标温度根据快冷得到的结构而设定为 。左右【】825K，因其晶化温度大致在830K和835K）分别设为temperature 中点击查看。以及得到其他信息，可在Project 衍射。再对上面的结构作X 3 Reflex 模块 Powder Diffraction工具调出如下图所示，或者从菜单工具，然后选择Ref

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