中科大MaterialsStudio培训教程4

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,2021/9/27,*,*,反应物,产物,过渡态(能量鞍点),反应势垒,反应能,反应坐标,2021/9/27,1,目的:介绍如何使用 DMol,3,和 Reaction Preview 工具进行过渡态搜索的计算。对简 单反应，这种方法是有效的。模块:Materials Visualizer,DMol,3,前提:用局域内坐标对固体进行结构优化。,用,LST/QST,搜索过渡态,背景,对任何反应的势能面的探索都要求知道反应进程中每一步的结构和能量，或者动力学和热动力学的快照(snapshots)。特别重要的是决定反应速率的那一步，这通常需要找到那些难以捕获的过渡态结构。有一些方法对找到过渡态的结构是很有效果的，其中比较知名的就是线性同步度越(linear synchronous transit,LST)和二次同步度越(quadratic synchronous transit,QST)。,本例中，我们将介绍DMol 中的LST/QST 工具的使用，将会看到如何使用LST/QST 搜索乙烯醇转变为乙醛的H转移反应 的过渡态结构。,CH2CHOH CH3CHO,本例包括以下内容：,1.建立一个计算模型,2.优化分子结构,3.定义原子对,4.用LST/QST 的方法计算过渡态,5.优化过渡态,结构,2021/9/27,2,1.建立一个计算模型,选择 creating a new project，建立名为vinylOH 的project。,在本单元中，你要在两个不同的3D Atomistic 界面中建立反应物和产物模型。第一步就是打开一个新的3D Atomistic界面，构建反应物乙烯醇(vinyl alcohol)。,点击工具栏里的New button，选择3D Atomistic。,在Sketch 工具条上，点击Sketch Atom 按钮 。将鼠标移至3D Atomistic界面，连续点击三次鼠标，画三个连接的碳原子。按一下键盘上的ESC 键。,2021/9/27,3,改为球棍显示,在,3D Atomistic,上，点击选择第三个C 原子。点击Modify Element 按钮上的选择箭头，选择氧元素。刚才被选的原子由碳原子变成了氧原子。,点击碳-碳键一次，选中。点击Sketch工具条上的Modify Bond Type 键 ，选择双键，从而把单键变成双键。点击别处，,取消,选择碳-碳键。,2021/9/27,4,按下Adjust Hydrogen 按钮 ，点击一次Clean 按钮 ，拖动结构模型，使得和下图相似，以球棍模型显示。,在,Project,浏览器内，右击,3D Atomistic.xsd,，选择,Rename,，将其重新命名为,reactant.xsd,。,在3D Viewer上选择Selection按钮 ，双击乙烯醇结构中的任何一个原子。,这样乙烯醇的每一个原子都被选上，颜色显示为黄色。,2021/9/27,5,*在这个新的3D 界面文件中，点击O-H 键。按下键盘上的DELETE 键。点击Sketch Atom 按钮 ，然后是孤立的H 原子，以及亚甲基团中的C 原子。如果画分子有问题，则删除3D Atomistic.xsd上的原子，再进行CTRL+C、CTRL+V。,*在键盘上按下CTRL+C。,选中的结构文件被复制到了剪贴板。,*用File/New打开一个新的3D Atomistic 文件，按下键盘上的CTRL+V。,*结构模型被粘贴到刚刚新打开的3D 界面上。现在需要改变化学键和对原子重新排布以得到产物结构。,2021/9/27,6,点击一次Clean 按钮 。,现在结构就和下面的看上去相似了。,点击一次C-O 键，由单键改为双键。,连续双击C-C 键，C-C 键就会由双键变为三键，然后又变成单键。,2021/9/27,7,现在需要把该结构的文件名改为product.xsd。右击工作浏览器Project内的3D Atomistic.xsd，将其名称改为product.xsd，回车。,2.优化分子结构,为了优化LST/QST的计算性能，需要对反应物和产物的结构进行优化。这个工作可以通过DMol3 的几何优化功能来完成。,点击别处，取消选择结构。按下工具条上的DMol3 按钮 ，然后选择下拉条中的Calculation。DMol3 的计算对话框显示出来。,将Task 由Energy 改为Geometry Optimization。确认Quality 设为Medium。将泛函改为GGA BP。,2021/9/27,8,刚才指定了使用的Hamiltonian 和计算的精度。精度决定了使用的基组(basis set)和轨道的截断cutoff。这里基组为DND。可以在Electronic 栏里检查这些参数的设置。现在需要应用电子分布热平滑thermal smearing来加快结构优化的收敛。,点击Electronic 标签。检查SCF 是不是设为Medium。按下More按钮，显示了DMol3 的Electronic 选项对话框。在SCF 标签栏里，勾选上Using smearing 选项。关闭DMol3 Electronic选项对话框。,2021/9/27,9,现在准备开始计算了。,让,reactant.xsd,成为当前工作文件。点击,Job Control,标签。按下,More,按钮，显示了,DMol3,的工作控制选项对话框。确认,Update structure，Update graphs,和,Update textual results,三项被勾选上。关闭,Job Control,选项对话框，点击,Run,按钮。,2021/9/27,10,当第一个计算结束后，对product.xsd 重复刚才的操作。,把当前工作文件换为product.xsd，点击DMol3 计算对话框上的Run 按钮。,计算过程中，计算的进程用图表和文本文件的形式展现出来。,当两个计算都完成的时候，两个新的文件夹出现在工作浏览器中，分别叫做,reactant DMol3 GeomOpt,和,product DMol3 GeomOpt,。最后的优化结构包含在reactant.xsd 和product.xsd 文件中，计算的输出结果在reactant.outmol 和product.outmol 文件中。,几何优化文件夹包含了.xtd 文件，这是能量最小化过程中的轨迹文件，可以显示几何优化过程。下面演示反应物的结构优化过程。从,reactant Energies.xcd图中可以看出，反应物经过12步才优化结束，我们可以看到每一步结构的变化。,2021/9/27,11,在Project中双击reactant.xtd文件，动画显示工具按钮Animation,激活。,如果动画(Animation)工具条,是不可见的，则按右侧的操作，使用观看(View)菜单让它显示。,2021/9/27,12,设置显示方式，按播放键 。,B-O近似，体系的电子能量是核构型的函数。,2021/9/27,13,在继续工作之前，需要关闭Materials Visualizer 中的所有文件。,关闭DMol3 计算对话框。选择FileSave Project，然后WindowClose All。,双击几何优化子文件夹中的reactant.xsd 和product.xsd。,现在工作区域中只有两个优化了的结构。,2021/9/27,14,3.定义原子对,用DMol3 进行过渡态搜索，反应物和产物的所有原子都必须配对对应。这个可以通过使用工具栏里的反应预览(Reaction Preview)功能来实现。,从菜单条中选择Window|Tile Vertically，使反应物和产物以肩并肩的形式显示。,2021/9/27,15,现在准备开始对反应物和产物结构中的原子进行对应配对(物质不灭),从菜单条中选择Tools|Reaction Preview，打开了反应预览(Reaction Preview)对话框。,从Reactant 的 下拉树形图中的几何优化文件夹中选择reactant.xsd。,2021/9/27,16,同样从Product 的下拉树形图中的几何优化文件夹中选择product.xsd。,按下Match按钮。寻找等价原子(Find Equivalent Atoms)对话框显示出来了，从中可以看到，一个原子(O)匹配了，而仍有六个原子(C、H)没有匹配。,2021/9/27,17,双击反应物栏(reactant column)中的2xC。,在产物栏里的对应的文件夹同时打开了。反应物栏包含了1:C 和2:C，它们应该直接和产物栏里的对应物相匹配，以下的步骤将对此加以确认。,分别点击反应物框里的1:C和产物框里的1:C。,两个对话框里的碳原子被选上，两个3D 文件里的碳原子也同时被选中。,2021/9/27,18,认为反应物和产物里的两个1:C原子是等价原子，点击,Auto Find,。,寻找等价原子(Find Equivalent Atoms)算法匹配了2个C和1个H原子。还有3个H原子未匹配。,2021/9/27,19,仍有3个H原子没有匹配，重复上面的手动匹配步骤。双击反应物栏(reactant column)中的3xH。在产物栏里的对应的文件夹同时打开了。反应物栏包含了4:H、5:H 和7:H，分别点击反应物框里的4:H和产物框里的4:H，两个对话框里的H原子被选上，两个3D 文件里的碳原子也同时被选中。,2021/9/27,20,点击,Set Match,，还剩下5H、7H未匹配。重复这个过程，继续点击,Set Match,来匹配剩下的没有配对的原子。,2021/9/27,21,反应物、产物的原子已配对。现在可以预览一下反应物和产物之间原子的匹配情况。,点击反应物或者产物栏中列表里的任意一个原子，可以看到匹配的另一原子。考察匹配情况，直到满意为止。关闭Find Equivalent Atoms 对话框。,2021/9/27,22,运用DMol3的 LST/QST功能来搜索过渡态，需要在反应物和产物之间创建一条通道，这也是DMol3 计算时所要求的输入文件。,在反应预览(Reaction Preview)对话框中，把桢数提高到100；勾选上Superimpose structures；单击Preview按钮，关闭反应预览(Reaction Preview)对话框。,在几秒钟内，一个名为reactant-product.xtd 的新的3D Atomistic Trajectory 文件显示出来；可以对这个文件进行DMol3 计算；可以使用动画(Animation)工具条来播放轨迹文件。,动画用Bounce模式观看效果最佳。把化学键监测(bond monitoring)打开，这样每变化一步，就会对化学键重新计算。,2021/9/27,23,点击Build|Bonds，打开Bond Calculation对话框，勾选上化学键计算对话框上的Monitor bonding，关闭对话框。,设置动画演示方式，按下动画工具条上的Play 按钮，观看反应物到产物的变化。,看完后按下Stop 按钮。,2021/9/27,24,2021/9/27,25,4.使用LST/QST/CG 方法计算过渡态,Note,：reactant_product.xtd 包含了DMol3 需要的重要信息，第一桢是反应物的，最后一桢是产物的。,现在准备设置使用DMol3 计算过渡态。,从菜单条中选择Modules|DMol3|Calculation，或使用DMol3图标 ，,DMol3 计算对话框就显示出来。,在设置(Setup)标签栏里，把Task 由几何优化改为TS Search。确定计算精度为Medium，泛函为GGA 和BP。,点击More按钮显示DMol3 过渡态搜索(DMol3 Transition Stat