class9DMol3离域内坐标法进行几何优化.ppt

使用DMol3中的离域内坐标对固体进行几何优化,DMol3的面向分子的离域内坐标优化机制为大分子系统提供了一套良好的方案。在MSModeling的DMol3中，这个机制被扩展到周期性系统。这个基于离域内坐坐标的新型优化工具还有能力处理以下体系：*高配位体系，比如密堆积固体。*片断系统，比如分子晶体，其中的内坐标并不是遍及整个优化空间。*优化过程中，部分进行了的迪卡尔坐标限制的体系。自身效率工作表明，对周期性体系而言，这个艺术级(state-of-the-art)的离域内坐标优化方案的效率要比迪卡尔坐标方法高出2-5倍，而迪卡尔坐标方法是现在进行固态计算的标准方法。在这部分讲义中，将介绍如何利用DMol3的的离域内坐标方法，对沸石结构进行几何优化。内容1.开始建立名为chabazite的Project，输入沸石结构2.建立DMol3计算任务3.控制工作设置和运行计算任务4.使用服务器控制台控制计算任务5.检验计算结果,目的：介绍volumevisualization工具和离域内坐标在周期体系中的应用模块：MaterialsVisualizer,DMol3前提：简单分子的绘制，用DMol3计算带结构和态密度,1.开始建立名为chabazite的Project，输入沸石结构,a.运行MS程序，选择Createanewproj，单击“确定”按钮。,b.找到Class7文件夹，写入Chabazite，单击OK按钮。这样就在ProjectExplorer中建立了名为chabazite的Project。,注：若MS已打开运行，则用File/NewProject另外建立一个名为chabazite的Project。,在Ddisk中建立Class7文件夹。,下一步是输入需要进行优化的沸石zeolite的结构。MSModeling提供了许多现成的沸石结构。本例中，将优化菱沸石(chabazite)。从工具栏选择FileImport，打开Import文件对话框，找到ExamplesDocuments3DModelCHA.xsd，然后点击输入文件对话栏上的Import按钮。,双击,Aunitcellofchabazite,在3D浏览器中右击鼠标，选择DisplayStyle，把显示方式改为Polyhedron。关闭对话栏。,2.建立DMol3计算任务点击工具栏上的DMol3按钮，选择Calculation，或者从菜单栏里选择Module/DMol3/Calculation，打开计算对话框。,从Task下拖列表中选择GeometryOptimization；把Functional设置从LDA和PWC改为LDA和VWN。LDA/VWN函数适合出来类似zeolite这样的共价晶体。虽然精度上不如GGA函数，但在矛盾的计算成本和精度中找到了择衷点。当把计算任务改为GeometryOptimization的时候，More按钮被激活，从而进行更多与此任务相关的设置。单击More按钮，打开DMol3几何优化对话栏。这里可以通过改变Quality的水平来设置收敛偏差(ConvergenceTolerance)或者单独编辑这些值。默认的设置是Medium，包括以下内容：能量为2.0e-5哈特里的收敛设置，0,004哈特里/埃的最大应力设置，以及0.005埃的最大位移设置。Quality设置选择了基矢、积分网格、收敛阈，这些选择使得在合适的计算成本下，大多数情况都能得到合理的结果。Quality设置为Medium，关闭DMol3GeometryOptimization对话栏。,4,在DMol3Calculation对话栏上，转换到Electronic标签。DMol3上的Electronic标签栏显示出来了，这个栏内包括了与电子哈密顿算符相关的参数。在这个任务中将使用DND基组和LDA/VWN汉密尔顿算符。,点击按钮More，显示DMol3ElectronicOptions对话栏。点击OrbitalCutoff标签，可以看到orbitalcutoff设置的值为3.5埃。关闭对话栏。提示：可以通过点击More按钮对SCF和k点电子汉密尔顿算符参数进行更加精确的操作和控制。关闭DMol3ElectronicOptions对话栏。,在设置计算条件时，我们可以改变一些选项的数值，以降低计算成本。在这个栏内，课改变的设置是k-pointset和Orbitalcuttoffquality。把k-pointset和Orbitalcuttoffquality分别设置为Gamma和Coarse。,点击DMol3Calculation对话栏上的Properties标签，Properties标签显示出来。结构优化后，可以对需要的某些性质进行计算。勾选上ElectronDensity和Electrostatics。,3.控制工作设置和运行计算任务可以使用JobControl标签上的命来来控制DMol3计算。点击DMol3Calculation对话栏上的JobControl标签。这里，我们可以选择把计算任务提交到网内的任何一台机器上，并设置不同的选项，这些选项包括计算任务描述、计算是否使用多个处理器运行和使用的处理器的数目等等。可以点击More按钮来对计算任务进行更多选择，包括实时更新设置和控制计算结束时的任务等等。现在可以运行DMol3计算任务了。点击Run按钮。,产生了一个名为status.txt的文件，里面含有DMol3运行状态。这个文件在计算任务结束以前会隔一段时间自动更新。不久之后，两个名为CHAEnergy.xcd和CHAConvergence.xcd的图表文件显示出来了，它们分别对应于计算的优化和收敛状态。这对于可视化监视计算进程非常有用。,工作浏览器JobExplorer开始工作，包括了计算状态等信息。,4.使用服务器控制台控制计算任务我们可以通过使用服务器控制台ServerConsole来查看实时状态。从菜单栏中选择Tool|ServerConsole。扩展ServerManagement和Jobs节点。计算任务被名为CHADMol3GeomOpt。点击CHADMol3GeomOpt。服务器控制台上结果界面显示了计算任务的实时状态。关闭ServerConsole。当计算任务结束的时候，相关文件就从服务器传回客户端。,ExpandtheServerManagementandJobsnodes,5.检验计算结果当计算结束的时候。计算结果返回到ProjectExplorer的CHADMol3GeomOpt文件夹里。*FileSaveproject，WindowCloseAll*双击CHADMol3GeomOpt文件夹中的CHA.xtd文件。,沸石CHA的结构显示在3D模型窗口中；这是一个包含了几何优化进程的轨迹文件。我们可以用Animation工具栏里的控制工具来浏览几何优化的历史进程。如果Animationtoolbar不是可见的，可以从菜单栏的ViewToolbarsAnimation选中。,*为了较清楚地观看沸石几何结构的优化过程，对动画播放做了设置，见右图。,在Animationtoolbar上，点击Play按钮来演示优化进程。结束演示时，点击Stop按钮,终止演示。,最终的总能量可以在.outmol文件中看到。,使用ProjectExplorer，换到CHA.outmol界面，向下拖动文件寻找总能量。由CHAEnergies.xcd文件或前面的动画可知，结构优化了7步。所以要找第7次优化后的能量。,*.xsd文件包含了优化结构，双击ProjectExplorer中的CHADMol3GeomOpt里的CHA.xsd文件。,*要使得计算所得的性质可视化，就要对输出文件进行分析。点击工具栏上的DMol3按钮，选择Analysis，或者从菜单栏里选择ModuleDMol3Analysis，打开DMol3分析对话框。,或,*选择ElectronDensity选项，确定显示在Results栏里的是正确的.outmol文件。设置DensityField设为TotalDensity，并且Viewisosurfaceonimport被勾选。点击Import。,*总电子密度显示在沸石上。放大和旋转以观察其结构。,换到Isosurface标签。我们可以使用Transparency标尺改变所显示的等密度面的透明程度，通过拖拽标尺改变Transparency。,我们可以改变等密度面的显示方式，之前我们需要改变原子的显示方式。在窗口上右击，在DisplayStyle中，点击Atom栏内的Stick。,*也可以通过选择Dots或者Solid选项来改变等密度面的显示方式。点击Dots按钮，然后换回到Solid按钮。把Transparency改回到初始值。,*还可以控制等密度面的大小把Isovalue的值改为0.1，按下TAB键；然后把Isovalue的值改回0.2，再按一回TAB键。,一旦得到了电子的等密度面，接下来我们可以在其上把其它的性质绘出来。从Module工具栏里打开DMol3Analysis对话栏；选择Potential；确定Viewisosurfaceonimport选项没有被勾选；点击Import按钮。当数据输入完成，关闭对话栏。,右击鼠标，选择DisplayStyle；换到Isosurface标签，点击与MapField区域相关的选项箭头并选中DMol3electrostaticpotential，关闭对话框。,电静势被绘在电子等密度面上。我们可以使用ColorMap对话栏改变所绘图的颜色。右击鼠标，选择ColorMaps。使用ColorMaps对话栏，我们可以改变颜色方案，所绘图的值以及显示的带数目。,改变Spectrum为Blue-White-Red。点击与From相关的right-arrow，选择MappedMinimum。点击与to相关的right-arrow，选择MappedMaximum，关闭对话框。等密度面上的颜色的改变反映出在电子密度图上所绘图的最大值和最小值。,我们可以从体积可视化工具栏内得到更多的体积可视化工具。从菜单栏里选择View|Toolbars|VolumeVisualization。,我们可以给结构加一个切面。点击与CreateSlicesbutton相关的选项箭头，选择BestFit。,在ChooseFieldstoSlice对话栏里，选择DMol3electrostaticpotential选项，点击OK。,一个穿越电静势的切面就显示出来了。这个切面使用了全区的最大值和最小值。如果要查看最大值和最小值在切面上的分布和具体区域，就需要再次使用ColorMap对话栏。,点击工具条上的ColorMaps按钮,将鼠标移到切面上，点击切面并选择。当你选择切面的时候，一条黄色虚线显示在切面的边沿，而在选中区域的中间有一个黄色虚线十字（若看不见，可旋转3D图）。,点击与From相关的right-arrow并选择MappedMinimum。点击与To相关的right-arrow并选择MappedMean。改变Spectrum为Blue-White-Red。有一些较大的最大值在这个网格区域，但这是不需要看到的，所以没有选择MappedMaximum，而选择了MappedMean。,我们还可以使用ColorMaps对话栏对这些值进行选择性显示。点击与颜色条相关的Redright-arrow，所有大于平均值的点都被从切面上删除。,我们还可以通过操作删除特定的颜色、改变颜色等等。点击按钮，平均值以上的颜色就重新显示出来了。关闭ColorMaps对话栏。,我们可以用DisplayStyle对话栏里的Slice工具改变切面的显示方式。右击鼠标，打开DisplayStyle对话栏。换到Slices标签，拖动Transparency标尺到右边。,当拖动Transparency标尺的时候，切面的透明度也随着变化。当结束对切面和等密度面等的分析工作之后，我们可以删除它们。将鼠标移动到切片上，选中Slice，则切片的边框变黄。按Delete，则删除切片。,将鼠标移动到等值面上，单击选中isosurface，则等值面变成黄色虚点。按Delete，则删除等值面。,WorkingwithVisualizer中的Creatingandcustomingslices(可在help中搜索)有slices的更多介绍。FileSaveProject，WindowCloseAll,