使用Multiwfn绘制NBO及相关轨道

1、使用Multiwfn绘制NBO及相关轨道文/SoberevaFirst release: 2012-Apr-4 Last update:2012-Apr-161前言NBO程序输出的轨道包括NBO （自然键轨道）、NHO （自然杂化轨道）、NAO （自然原子轨道）、NLMO （自然定域化分子轨道），以及前头带P的类型，如PNBO （初自然键轨道）。NBO程序输出的plot文件包含了这些轨道的一切信息。可视化这些轨道有很多方法：（1）对于 Gaussian 用户，可以用 pop=saveNBO 或 pop=saveNLMO 将 NBO 和 NLMO 轨道写入 chk文件中，再用Gaussview观

2、看。但是对于其它类型轨道就没辙了。（2）用ChemCraft。这个软件是收费的，而且我曾发现在基组含有高角动量函数时这个程序给出的图形是明显错误的。（3）用NBOView。这个程序也是收费的，而且界面设计非常落伍，操作不便。（4）用Zork编写的NBO2molden程序将NBO plot文件转换为Molden的输入文件，然后用Molden程序去看。然而Molden这个程序我认为很不好用。Multiwfn是一款免费、开源的波函数分析程序，很早就开始支持读入NBO plot文件来绘制 前述各种类型轨道，可以很方便地显示等值面图，也能绘制出各种各样漂亮的截面图，操作 精品文档很方便，可以从下载。从2

3、.3.3版开始还可以同时绘制两 条轨道，对于分析轨道间重叠很有用。在绘制NBO相关轨道方面Multiwfn是最佳选择。本 文将介绍如何用Multiwfn绘制NBO轨道，其过程也同样适用于绘制分子轨道、自然轨道等。2输入文件这一节介绍一下怎样产生能够让Multiwfn显示NBO及相关轨道的输入文件Gaussian的fch文件和NBO plot文件都可以用作Multiwfn的输入文件。对于 Gaussian 用户，可以在 route section 里写上 pop=saveNBO 或 pop=saveNLMO,这样 NBO 轨道和NLMO就会分别代替分子轨道储存在chk文件里。然后用formchk

4、程序将chk文件转 化为fch文件，就可以读入Multiwfn来看NBO和NLMO 了。注意fch文件中的轨道编号和 Gaussian的NBO 3.1模块(L607)显示的轨道编号不全一致，在存入fch文件时会进行重排， NBO模块输出信息末尾的Sorting of NBOs和Reordering of NBOs for storage部分显示 的序号是重排前后的对应关系。如果想看NBO和NLMO以外的轨道，就只能用NBO plot文件了。做法是在Gaussian输入文件中加上pop=nboread关键词，在末尾空一行写上比如$NBO plot file=c:ltwdNH2COH$END，用

5、Gaussian 运行后就会有一批 NBO plot 文件 NH2COH.31、NH2COH.32 . NH2COH.41在c:ltwd下面生成。.31文件储存的是基函数信息，.32.40分别储存的是PNAO/NAO/PNHO/NHO/PNBO/NBO/PNLMO/NLMO的展开系数信息。.41是密度矩阵信息，对绘制轨道没直接用处。如果用的是独立版本的NBO程序，即GENNBO,比如想在h:Yuri目录下生成NBO plot文件的话就将GENNBO输入文件(.47)开头的$NBO和$END中间写上plot file=h:YuriNH2COH然后用GENNBO运行即可。本文的例子都是B3LYP/

6、6-31G*下的甲酰胺，Gaussian输入文件如下 %chk=C:gtestNH2COH.chk# B3LYP/6-31G* opt pop=nboreadB3LYP/6-31G* opt0 1C-0.03549095-0.457814140.00000000H-0.01702195-1.527654730.00000000O-1.286064210.235703360.00000000N1.073468370.208226550.00000000H1.056207621.208077570.00000000H1.94799512-0.276750720.00000000$NBO plot

7、file=c:NH2COH $END3绘制轨道等值面图启动Multiwfn，然后输入.fch文件或.31文件的路径。如果用的是后者，程序会让你再输 入.32.40文件中的一个，输入哪个取决于要绘制哪种轨道，比如要看PNHO，就要输入.34 文件的路径，假设你先输入的是c:Rio_rainbow_gateRioXLinaNH2COH.31，而.34文件就 在这个.31文件所在路径下，那么直接输入数字34就可以了，可以免得输入完整的路径。接下来输入0就会在文本窗口中输出所有原子的坐标，并弹出用于观看分子结构和轨道等值面的窗口，如下所示，图中显示的是12号NBO。此界面的各个物件的用途稍微玩弄一下就

8、明白了，在手册3.2节也有详细说明。右下角是轨 道列表，点击其中一个就会马上在图形窗口中显示相应的等值面，十分方便省事，等值面生 成速度也比起Gaussview等程序快得多。注意，如果使用.fch作为输入文件，如前所述其 中轨道序号已经被重排，所以选择轨道列表中的第X号所显示的未必是NBO输出信息中的第X号轨道的等值面。从NBO模块的二阶微扰能的分析结果看，NBO 12（N的孤对电子轨道）和NBO 56（C-O的反 兀键）之间有很强的相互作用，这能够使体系能量下降62.8kcal/mol。像这种情况，同时 绘制两条轨道的图形来分析它们的相位交叠是很有意义的。为了做这样的图，我们先从轨道 列表中

9、选 12，然后选中 Show+Sel. isosur#2复选框（全称是 Show and select isosurface#2），再从列表中选56,这两条轨道就都显示出来了，如下所示没有选中Show+Sel. isosur#2复选框时选择的轨道所产生的等值面被称为Isosurface #1，选中了此复选框再选的轨道所产生的等值面被称为Isosurface#2。只有已经出现了Isosurface #1时这个复选框才允许选。如果取消选择此复选框，则已出现的Isosurface #2将会消失，若再次选中此复选框，之前消失的那个Isosurface #2又会重新出现。Isosurface #1的正值

10、和负值部分的等值面用绿色和蓝色表示。为了区分，Isosurface #2的正值和负值部分等值面分别用黄绿色和紫色表示。等值面的风格可以通过界面上方Isosur#1 style和Isosur#2 style中的相应选项设定，可选的风格有：不透明面、网、点、 不透明面+网这四种，在下个版本中还会加入透明面的风格。不透明面的颜色以及网/点的颜 色也都可以通过Isosur#1 style和Isosur#2 style中的相应选项进行设定，程序会让你输 入R,G,B（红、绿、蓝）分量的值，每个分量值的范围应在0.01.0以内，例如0,1,0就代 表绿色、1,1,0就代表黄色。R,G,B值需要输入两次，第

11、一次是设定等值面正值部分的颜色， 第二次是设定负值部分的颜色，文本框内直接出现的数值是当前R,G,B值。风格和颜色在设 定后会立刻在屏幕上生效。上面图中的NBO 12和NBO 56的等值面都是不透明的，分析重叠程度比较困难，将两个等值面都设为网状风格后如下所示，交叠区域看起来清晰了。它们之间存在很明显的交叠，是它们之间二阶微扰能很大的主要原因之一。实际上，在显示轨道等值面之前，程序会先计算出一个涵盖整个分子空间范围的轨道波函数的格点数据，这个格点数据的格点密度越大，即相同空间范围内的点数越多，则等值面越精细。为了生成等值面速度比较快，默认的点数并不多。尤其当体系较大时，由于格点密度往往较稀疏，

12、等值面会显得不够光滑，可以通过增加格点数来提高等值面显示效果。方法是选择窗口上方的Isosur. quality，点击其中的按钮后可以输入所期望的格点数。重新设定格点数之后，目前显示的等值面都会自动消除，之后再选择轨道所生成的等值面都会套用刚才输入的设定。4生成轨道波函数格点文件并用VMD绘图虽然直接用Multiwfn来产生NBO及相关轨道等值面图一般够用了，但是为了追求更好的显示效果且不嫌麻烦，可以借助于更专业的能够显示格点数据的程序，这里推荐VMD，可以在 HYPERLINK /Research/vmd/%e5%85%8d%e8%b4%b9%e4%b8%8b%e8%bd%bd%ef%bc%

13、8c%e6%9c%ac%e6%96%87%e7%94%a8%e7%9a%84%e6%98%af /Research/vmd/免费下载，本文用的是 1.9 版。这里还是以同时显示NBO 12和NBO 56的等值面图为例。首先需要用Multiwfn分别生成NBO12和NBO 56对应的.cub格点文件。启动Multiwfn，输入.fch或NBO plot文件名后依次输入5/生成格点数据4/轨道波函数12 /NBO 122中等质量格点数据2/将生成的格点数据保存到当前目录下MOvalue.cub中。然后我们将此文件改名为NBO12.cub0返回主界面5/生成格点数据4/轨道波函数56 /NBO 56

14、中等质量格点数据2/将生成的格点数据保存到当前目录下MOvalue.cub中。然后我们将此文件改名为NBO56.cub现在关闭Multiwfn，启动VMD，将NBO12.cub和NBO56.cub依次拖进VMD主窗口。选主界面 中的 Graphics-representations，将新窗口上方的 Selected Molecule 切换到 NBO 12。将 已有的那个显示方式的Drawing Method改为CPK。然后点Create Rep，将新产生的显示方 式的 Drawing Method改为 Isosurface， Draw 改成 Solid Surface， Show 改成 Iso

15、surface， Isovalue改为0.05，将Coloring Method改为ColorID并在旁边选7 Green。然后再以相 同的方式建一个显示方式，但ColorID选0 Blue，Isovalue改为-0.05。这时NBO 12的等 值面就出现了。之后在Selected Molecule中切换到NBO56，按照处理NBO 12的方法也通 过建立两个显示方式来显示NBO 56的等值面。得到的图像如下在VMD中这两个轨道的等值面还可以用透明的方式显示。将前面建立的四个显示方式中的Material都设为GlassBubble，然后在主界面里面选Graphics-Colors-Displa

16、y-Background-8 white 将背景变为白色，在主界面里选Display-Rendermode-GLSL，就能看到希望的效果。如果想获得更好效果可以用渲染器(当GLSL无法打开时，也只能通过渲染器才能获得透明效果)。例如，在主界面的File-Render里面选 Tachyon (internal, in-memory rendering)然后点 Start Rendering，就会出现下面的图，重叠区域显示得相当清楚。在Multiwfn里不能同时显示多于两条轨道的等值面（因为这种情况一般很少涉及，所以不 打算支持）。如果想同时显示三条及以上轨道的等值面的话，就只能用VMD来实现。方

17、法很 简单，就是将更多的轨道的cub文件拖进VMD里并进行同样的设置即可，想同时显示多少条轨道都没问题。5生成轨道波函数等值线图Multiwfn能够绘制的平面图种类很多，包括填色图、等值线图、地形图、梯度线图、向量 场图。绘制过程很简单，选项也很灵活，在Multiwfn手册的4.4节里给出了很多实例。这 里只介绍一下如何作同时含有两个轨道的等值线图（其它类型的图只能一次做一个轨道的）。 精品文档这次还是作NBO 12和NBO 56的图，对于表现它们的交叠情况，最合适的作图平面应当是垂 直于分子面且穿过N和C的那个平面。这个面不是XY/YZ/XZ面之一，也没法用三个原子坐 标来定义。在Multi

18、wfn里定义这个面最好通过指定的三个坐标点来定义。第一个和第二个 点的坐标就设为C和N的坐标，而第三个点的坐标设为在C或者N的坐标的基础上往Z方向稍微移动一点（Z轴垂直于分子平面）。启动Multiwfn，载入.31和.37文件（或载入.fch文件），然后依次输入4绘制平面图4/轨道波函数12,56两个NBO轨道的编号。如果只输入一条轨道的编号，做出来的图就是一条轨道的直接敲回车，用默认的格点设定5/通过输入三个坐标点来定义作图平面0.0000000.7940890.000000/第一个坐标点，即C的位置。X/Y/Z坐标可以用逗号或空格来分隔，单位是bohr。在进入主功能0的时候屏幕上就会出现各

19、个原子的坐标，可以直接将其拷贝下来粘贴到此处，如果不知道怎么拷贝，可以参见手册5.4节-1.778942-1.0642350.000000/第二个坐标点，即 N 的位置-1.778942-1.0642351.000000/第三个坐标点，是在N的位置的基础上往Z正方向挪了 1 bohr得到的。挪的距离需要反复尝试找到最佳的，不合适的值会使感兴趣的区域不在图的中央。得到的图是如下这样的，实线和虚线代表正值和负值部分这个图看起来比较乱，这是因为在默认的等值线设定下，数值比较小的等值线也显示了出来，然而这些较小数值的等值线的意义并不大。为了让图看起来比较清楚，应该删除一些数值较小的等值线。因此在图上点

20、鼠标右键关闭之，然后输入设定等值线5删除一批等值线1,4/删除1至4号等值线，即分别为0.001, 0.002, 0.004, 0.008的四条5再删除一批等值线28,31 /删除28至31号等值线，即分别为-0.001, -0.002, -0.004, -0.008的四条。如果嫌每次作类似的图都要删等值线比较麻烦，可以接下来选6来将当前等值线设定保存到指定的外部文本文件里，下次再进入这个设定界面时可以选7来从指定的外部文件中读入等值线设定1/退出等值线设定界面2令等值线数值显示在图中25/等值线数值的文字大小设为25-1/重新绘制图像此时得到如下图像，可见绝对值大小低于0.008的等值线都没

21、了，图像也变得十分清楚，很方便分析交叠区域。只有与绘图平面垂直距离小于特定距离的原子的符号才会显示在图中，这个距离阈值以及符号大小分别由settings.ini里的disshowlabel和pleatmlabsize参数调节，修改后需重新启动Multiwfn方可生效。6结合Photoshop同时绘制两条以上轨道的等值线图个别时候需要分析一个NBO轨道与多个NBO轨道的相互作用，这时可能需要同时绘制三条或更多条的轨道的等值线图。虽然Multiwfn不直接支持这种情况，但是通过利用Photoshop（以下简称ps），可以很容易地实现，而且借助于ps强大的功能，在线条风格上可以更自由地控制。这个例子

22、中，我们要将NH2CHO的NBO 4的等值线图利用ps叠加在上一节得到的NBO 12和NBO 56的等值线图上。NBO 4是C-N间的a成键轨道。首先，我们先得到只含NBO 4的等值线图。启动Multiwfn，载入NBO plot文件然后依次输入绘制平面图4/轨道波函数4 /NBO 42/等值线图直接敲回车通过三个点定义绘图平面。这三个点的坐标必须和上一节用的一模一样，只有这样得到的图才能精确地叠加到上一节的图上0.0000000.7940890.000000-1.778942-1.0642350.000000-1.778942-1.0642351.000000接下来还是和上一节一样，进入等值

23、线设定界面，删掉数值较小的等值线，然后退回到上一级菜单，最后选0将图片保存到当前目录下前缀为DISLIN的png格式的图形文件中。所得图像如下所示9.95J|am|I|I|I|1.643.394.936.586.239.8711.51Length unit: Echr& 33 -2840.00-P0.00现在打开ps准备将这个图和上一节的图合并，我这里用的是Photoshop CS2版。先将这两 幅图都拖进ps里，激活NBO4的窗口，按Ctrl+A全选并按Ctrl+C复制，切换到NBO 12+NBO 56的窗口中按Ctrl+V粘贴。这时会有两个图层，NBO 12+NBO 56的图层被NBO 4的图层覆 盖住了，为了能让前者也同时显示出来，就必须把NBO4图层的白色背景删掉以变成透明的 背景。最便捷的方法就是先确保图层列表里已经选定了 NBO 4的图层，然后选ps主菜单的 Select-ColorRange，然后将光标移到图上（会变成取色器形状的指针），点一下图中 的白色部分，确认Color Range窗口中的Fuzziness为0，然后点OK，这时NBO 4的图层的 白色背景就都被选中了，按一下键盘上的del键，背景就透明了，NBO4、12、56的等值线就同时显示出来了。在图上空白处点一下鼠标来取消选择状态，就能看到这样的图:不过，两条轨道以上等值线同时