中科大MaterialsStudio培训教程4

反应物产物过渡态(能量鞍点)反应势垒反应能反应坐标目的:介绍如何使用DMol3和ReactionPreview工具进行过渡态搜索的计算。对简单反应，这种方法是有效的。模块:MaterialsVisualizer,DMol3前提:用局域内坐标对固体进行结构优化。用LST/QST搜索过渡态背景对任何反应的势能面的探索都要求知道反应进程中每一步的结构和能量，或者动力学和热动力学的快照(snapshots)。特别重要的是决定反应速率的那一步，这通常需要找到那些难以捕获的过渡态结构。有一些方法对找到过渡态的结构是很有效果的，其中比较知名的就是线性同步度越(linearsynchronoustransit,LST)和二次同步度越(quadraticsynchronoustransit,QST)。本例中，我们将介绍DMol中的LST/QST工具的使用，将会看到如何使用LST/QST搜索乙烯醇转变为乙醛的H转移反应的过渡态结构。CH2CHOH→CH3CHO本例包括以下内容：􀁺1.建立一个计算模型􀁺2.优化分子结构􀁺3.定义原子对􀁺4.用LST/QST的方法计算过渡态􀁺5.优化过渡态结构1.建立一个计算模型选择creatinganewproject，建立名为vinylOH的project。在本单元中，你要在两个不同的3DAtomistic界面中建立反应物和产物模型。第一步就是打开一个新的3DAtomistic界面，构建反应物乙烯醇(vinylalcohol)。点击工具栏里的Newbutton，选择3DAtomistic。在Sketch工具条上，点击SketchAtom按钮。将鼠标移至3DAtomistic界面，连续点击三次鼠标，画三个连接的碳原子。按一下键盘上的ESC键。改为球棍显示在3DAtomistic上，点击选择第三个C原子。点击ModifyElement按钮上的选择箭头，选择氧元素。刚才被选的原子由碳原子变成了氧原子。点击碳-碳键一次，选中。点击Sketch工具条上的ModifyBondType键，选择双键，从而把单键变成双键。点击别处，取消选择碳-碳键。按下AdjustHydrogen按钮，点击一次Clean按钮，拖动结构模型，使得和下图相似，以球棍模型显示。在Project浏览器内，右击3DAtomistic.xsd，选择Rename，将其重新命名为reactant.xsd。在3DViewer上选择Selection按钮，双击乙烯醇结构中的任何一个原子。这样乙烯醇的每一个原子都被选上，颜色显示为黄色。*在这个新的3D界面文件中，点击O-H键。按下键盘上的DELETE键。点击SketchAtom按钮，然后是孤立的H原子，以及亚甲基团中的C原子。如果画分子有问题，则删除3DAtomistic.xsd上的原子，再进行CTRL+C、CTRL+V。*在键盘上按下CTRL+C。选中的结构文件被复制到了剪贴板。*用File/New…打开一个新的3DAtomistic文件，按下键盘上的CTRL+V。*结构模型被粘贴到刚刚新打开的3D界面上。现在需要改变化学键和对原子重新排布以得到产物结构。点击一次Clean按钮。现在结构就和下面的看上去相似了。点击一次C-O键，由单键改为双键。连续双击C-C键，C-C键就会由双键变为三键，然后又变成单键。现在需要把该结构的文件名改为product.xsd。右击工作浏览器Project内的3DAtomistic.xsd，将其名称改为product.xsd，回车。2.优化分子结构为了优化LST/QST的计算性能，需要对反应物和产物的结构进行优化。这个工作可以通过DMol3的几何优化功能来完成。点击别处，取消选择结构。按下工具条上的DMol3按钮，然后选择下拉条中的Calculation。DMol3的计算对话框显示出来。将Task由Energy改为GeometryOptimization。确认Quality设为Medium。将泛函改为GGA\BP。刚才指定了使用的Hamiltonian和计算的精度。精度决定了使用的基组(basisset)和轨道的截断cutoff。这里基组为DND。可以在Electronic栏里检查这些参数的设置。现在需要应用电子分布热平滑thermalsmearing来加快结构优化的收敛。点击Electronic标签。检查SCF是不是设为Medium。按下More…按钮，显示了DMol3的Electronic选项对话框。在SCF标签栏里，勾选上Usingsmearing选项。关闭DMol3Electronic选项对话框。现在准备开始计算了。让reactant.xsd成为当前工作文件。点击JobControl标签。按下More…按钮，显示了DMol3的工作控制选项对话框。确认Updatestructure，Updategraphs和Updatetextualresults三项被勾选上。关闭JobControl选项对话框，点击Run按钮。当第一个计算结束后，对product.xsd重复刚才的操作。把当前工作文件换为product.xsd，点击DMol3计算对话框上的Run按钮。计算过程中，计算的进程用图表和文本文件的形式展现出来。当两个计算都完成的时候，两个新的文件夹出现在工作浏览器中，分别叫做reactantDMol3GeomOpt和productDMol3GeomOpt。最后的优化结构包含在reactant.xsd和product.xsd文件中，计算的输出结果在reactant.outmol和product.outmol文件中。几何优化文件夹包含了.xtd文件，这是能量最小化过程中的轨迹文件，可以显示几何优化过程。下面演示反应物的结构优化过程。从reactantEnergies.xcd图中可以看出，反应物经过12步才优化结束，我们可以看到每一步结构的变化。在Project中双击reactant.xtd文件，动画显示工具按钮Animation激活。如果动画(Animation)工具条是不可见的，则按右侧的操作，使用观看(View)菜单让它显示。设置显示方式，按播放键。B-O近似，体系的电子能量是核构型的函数。在继续工作之前，需要关闭MaterialsVisualizer中的所有文件。关闭DMol3计算对话框。选择File\SaveProject，然后Window\CloseAll。双击几何优化子文件夹中的reactant.xsd和product.xsd。现在工作区域中只有两个优化了的结构。3.定义原子对用DMol3进行过渡态搜索，反应物和产物的所有原子都必须配对对应。这个可以通过使用工具栏里的反应预览(ReactionPreview)功能来实现。从菜单条中选择Window|TileVertically，使反应物和产物以肩并肩的形式显示。现在准备开始对反应物和产物结构中的原子进行对应配对(物质不灭)从菜单条中选择Tools|ReactionPreview，打开了反应预览(ReactionPreview)对话框。从Reactant的下拉树形图中的几何优化文件夹中选择reactant.xsd。同样从Product的下拉树形图中的几何优化文件夹中选择product.xsd。按下Match…按钮。寻找等价原子(FindEquivalentAtoms)对话框显示出来了，从中可以看到，一个原子(O)匹配了，而仍有六个原子(C、H)没有匹配。双击反应物栏(reactantcolumn)中的2xC。在产物栏里的对应的文件夹同时打开了。反应物栏包含了1:C和2:C，它们应该直接和产物栏里的对应物相匹配，以下的步骤将对此加以确认。分别点击反应物框里的1:C和产物框里的1:C。两个对话框里的碳原子被选上，两个3D文件里的碳原子也同时被选中。认为反应物和产物里的两个1:C原子是等价原子，点击AutoFind。寻找等价原子(FindEquivalentAtoms)算法匹配了2个C和1个H原子。还有3个H原子未匹配。仍有3个H原子没有匹配，重复上面的手动匹配步骤。双击反应物栏(reactantcolumn)中的3xH。在产物栏里的对应的文件夹同时打开了。反应物栏包含了4:H、5:H和7:H，分别点击反应物框里的4:H和产物框里的4:H，两个对话框里的H原子被选上，两个3D文件里的碳原子也同时被选中。点击SetMatch，还剩下5H、7H未匹配。重复这个过程，继续点击SetMatch来匹配剩下的没有配对的原子。反应物、产物的原子已配对。现在可以预览一下反应物和产物之间原子的匹配情况。点击反应物或者产物栏中列表里的任意一个原子，可以看到匹配的另一原子。考察匹配情况，直到满意为止。关闭FindEquivalentAtoms对话框。运用DMol3的LST/QST功能来搜索过渡态，需要在反应物和产物之间创建一条通道，这也是DMol3计算时所要求的输入文件。在反应预览(ReactionPreview)对话框中，把桢数提高到100；勾选上Superimposestructures；单击Preview按钮，关闭反应预览(ReactionPreview)对话框。在几秒钟内，一个名为reactant-product.xtd的新的3DAtomisticTrajectory文件显示出来；可以对这个文件进行DMol3计算；可以使用动画(Animation)工具条来播放轨迹文件。动画用Bounce模式观看效果最佳。把化学键监测(bondmonitoring)打开，这样每变化一步，就会对化学键重新计算。点击Build|Bonds，打开BondCalculation对话框，勾选上化学键计算对话框上的Monitorbonding，关闭对话框。设置动画演示方式，按下动画工具条上的Play按钮，观看反应物到产物的变化。看完后按下Stop按钮。4.使用LST/QST/CG方法计算过渡态Note：reactant_product.xtd包含了DMol3需要的重要信息，第一桢是反应物的，最后一桢是产物的。现在准备设置使用DMol3计算过渡态。从菜单条中选择Modules|DMol3|Calculation，或使用DMol3图标，DMol3计算对话框就显示出来。在设置(Setup)标签栏里，把Task由几何优化改为TSSearch。确定计算精度为Medium，泛函为GGA和BP。点击More…按钮显示DMol3过渡态搜索(DMol3TransitionStateSearch)对话框。确认搜索协议(Searchprotocol)设置为CompleteLST/QST，精度为Medium。关闭DMol3TransitionStateSearch对话框。电子Hamiltonian的设置与几何优化计算的设置一样。这次需要计算频率(Frequency)相关的性质。点击Properties标签栏，勾选上Frequency。最后，需要对工作描述(JobDescription)加以设置。点击JobControl标签，确认Automatic没有被勾选上；在JobDescription一栏里打上TS。按下Run按钮。关闭DMol3Calculation对话框。等待计算完毕。Note：在计算期间，数个不同的文件和一个LST/QST图显示在工作区。它们报告了计算的状态。特别地，通过显示energyvs.LST、QST和CG(conjugategradient，共轭梯度)这些路径坐标(pathcoordinate)的变化曲线，LST/QST图监测了过渡态搜索的进程。可以在TS.outmol文件里看到计算文本结果。如果该文件没有自动显示出来，双击工作浏览器内的TS.outmol，按下CTRL+F，搜索能量势垒。由TS.outmol中的数据可知，反应能量大约为-15kcalmol-1，反应势垒大约为50kcalmol-1。当LST/QST计算完成，打开TS.xsd文件就可以看到过渡态。如果文件没有自动显示出来，双工作浏览器(ProjectExplorer)内的TS.xsd。一个3DA