Multiwfn 3.0波函数分析程序的意义、功能与用途

文/Sobereva  2013-Apr-29

前言：Multiwfn波函数分析程序的功能极其丰富，用途十分广泛，用法十分灵活，所以绝大多数用户可能不容易完全了解Multiwfn的全貌。虽然在J. Comp. Chem. 33, 580-592 (2012)笔者曾经对Multiwfn的各个方面进行过系统的介绍，但那篇文章对应于2.1.2版，而撰文时最新的Multiwfn版本，即3.0版，已经比2.1.2版有了许多十分重要的扩展和改进。为了帮助用户全面认识Multiwfn，本文将对Multiwfn 3.0版的功能和用途进行一个既系统又简明扼要的概述。注意本文涉及的个别功能是在将要发布的3.0.1版里才有的。


1 波函数分析与Multiwfn的意义

量子化学本质上就是用基于量子力学的方法研究化学问题。通过现成的量子化学程序，在各种近似条件下求解薛定谔方程，就可以得到体系的波函数。波函数包含了体系中一切信息，可以视为是个黑箱，通过对波函数及其衍生信息（如电子密度）进行各种方式的分析，就可以从这个黑箱中提取各种具有化学意义的信息，使体系的内在特征能够被我们透彻地认知，诸如原子间键的强度与本质、电子在体系中的行为和分布方式等，同时也可以预测当前体系如何与其它体系相互作用。波函数分析的意义可以用下图来形象地描述



主流量子化学程序也可以给出体系的一些信息，如能量、稳定结构/过渡态结构、极化率、溶剂化能、NMR/IR/Raman/UV-Vis/ECD/VCD谱等，但这些信息不是基于直接对波函数分析而给出，故它们的计算一般不纳入波函数分析范畴。

波函数分析是一大类分析方法的集合。实际上，众所周知的前线轨道理论就是波函数分析方法的一种，它直接利用了前线轨道的波函数来分析化学反应性。量子化学界比较知名的两大类分析方法，即Natural bond orbital(NBO)分析和Atoms in molecules(AIM)分析，都是波函数分析方法当中重要的组成部分。各种形式的布居分析、电荷转移分析、电子定域性分析、键级分析、静电势分析等等，以及概念密度泛函理论的诸多内容也都属于波函数分析的范畴，这些后文将会具体列举。

波函数分析无疑是非常重要的。量子化学研究，特别是应用型的研究，决定研究价值的关键往往不在于计算过程，而在于分析得是否深入透彻，是否能给出足够的分析依据来支持自己的结论，是否能通过详细的分析发现新的问题并进而予以合理的解释。然而，波函数分析的普及程度目前相当令人沮丧，这在我国尤甚。从很多文章中，以及国内的量子化学有关的论坛上的讨论中可以明显看出，很多量子化学工作者对波函数分析知之甚少，导致很多人研究体系时感觉没得可分析，从而导致文章深度、档次上不去。也有很多研究者则想要解释体系的一些特点或论证一些设想（比如想说明哪个键比较强、电子跃迁到激发态后电荷是如何转移的，想解释为什么反应容易发生在某个位点），想去分析可是又苦于没有手段。还有不少人，看过他人的文章，知道应该做什么方式的波函数分析，然而却找不到现成、好用的程序来实现。

之所以意义重大、研究者们迫切需求的波函数分析的普及程度非常低，应用价值受到严重制约，主要有两个方面的原因：
(1)长期以来量化教材当中很少有书对波函数分析进行介绍，书籍内容的取材往往比较陈旧和保守，也限于大多数作者的知识水平，波函数分析方法方面往往顶多也就是介绍一下半个世纪前就提出的Mulliken布居，导致量子化学初学者们连接触波函数分析的机会都甚少。另外，虽然和波函数分析有关的文章数目甚巨，相关专著有不少，这个领域的研究从几十年前起就一直十分活跃，然而迄今国内外没有任何一本书籍对各种波函数分析方法进行全面而且深入浅出的介绍，资料过于零散，使得广大量化工作者学习波函数分析比较困难。
(2)现有的波函数分析程序基本都存在一种或多种弊病，导致广大研究者们很难将波函数分析方法付诸于实际研究。这些弊病包括：1 收费、不开源 2 计算速度慢而且没有并行化 3 使用、操作十分复杂难懂，输入文件编写麻烦，输出信息冗杂，不人性化，也没有像样的教程，更不提供技术支持 4 不够灵活，无法通过调整设置满足特殊研究需求 5 功能单一，每个程序往往只能实现很少的功能，缺乏整合性 6 作者别有用心，程序不公开发布，只给合作者用。另外值得一提的是，虽然各种半经验/从头算/第一性原理程序自身往往都带有一定的波函数分析功能（特别是Mulliken布居分析几乎都有），但是支持的分析甚少，而且往往不好用。

笔者开发的Multiwfn程序着眼于解决上述第(2)类当中的全部问题。虽然Multiwfn的开发目的有许多，但最主要的目的，是提供给广大科研工作者一个免费开源，并且是世界上功能最全面最强大、使用最方便、最容易上手、计算速度最快、最灵活的波函数分析程序。Multiwfn程序从2009年11月发布的1.0版开始一直不断开发至今，目前最新的3.0版已经很成熟，基本达到了预期目标，给波函数分析在实际研究中的普及扫除了重大障碍。如今Multiwfn的用户已经遍及世界各地，除我国外还包括美国、墨西哥、欧洲各国、俄罗斯、印度乃至伊朗等国，已经被百篇高水平SCI期刊上的研究文章以及书籍所引用，其中有不少发表在JACS、Angew。

Multiwfn的手册编写十分详尽，260多页的手册中不仅介绍了程序的具体使用，还将程序各种功能涉及到的所有的理论进行了介绍，还包含了约60个实际应用例子。结合笔者以往陆续在网上发布的一系列介绍各种波函数分析方法的帖子，其实已经解决了不少上述第(1)类问题。笔者有意日后专门撰写一本名为《量子化学波函数分析--原理与实践》的书来彻底解决这类问题，其中将对各种主流的波函数分析方法进行全面系统的讲解，尽可能压低阅读的门槛，同时结合Multiwfn等程序的具体操作步骤给出研究实例，使读者明白波函数分析方法都有哪些、原理是什么、能解决什么问题、在实际研究中如何具体使用。想必此书出版后将会对波函数分析方法在我国的推广有重大的意义。另外笔者在不久的将来还将要开办“波函数分析暨Multiwfn程序培训班”来亲自对波函数分析和Multiwfn的使用进行讲解。



2 Multiwfn的主要亮点

Multiwfn相对于其它种类繁多的波函数分析程序主要有以下几大亮点：

(1) 功能非常全面。Multiwfn是迄今功能最全面最强大的波函数分析程序，它将各种波函数分析方法高度、有机地整合到了一起。除了几个和NBO有关的分析外，所有最重要的波函数分析方法都可以仅靠Multiwfn一个程序来完成。Multiwfn的推出已经使得好几十个波函数分析程序（不少还是收费的）彻底失去了存在价值。同时Multiwfn的许多功能都是独一无二的，无法靠任何其它现有程序实现。

(2) 非常易于使用。Multiwfn程序是交互式程序，界面设计得十分周到并且高度智能化，用户完全不需要像使用大多数量子化学程序那样需要编写复杂的输入文件（关键词还容易忘，老得查手册很麻烦），而只要仔细按照屏幕提示一步步进行就行了，使用非常方便。大多数分析只需要敲几个键就能轻松完成。另外，Multiwfn的输出的信息简明易懂，不会输出一大堆只有内行甚至开发者才能搞明白的信息。并且手册里还提供了60个左右的例子帮助用户使用，其中涵盖了多数Multiwfn的最主要的功能。

(3) 高效。Multiwfn的开发过程特别注重提高运算效率以及降低内存的使用。并且几乎所有运算量大的部分都用OpenMP方法进行了并行化。对于计算较为耗时的波函数分析，Multiwfn的运算速度往往高于其它功能类似的程序一个甚至两个数量级。

(4) 结果可直接可视化。对于需要可视化分析的波函数分析方法，Multiwfn会自动调用DISLIN图形库直接将结果绘制成图像，这使分析过程大大简化，这对于实空间函数的分析尤为重要，并且很多绘图设定在Multiwfn中都允许由用户自行调节。而大多数量子化学程序都做不到这一点，用户往往得把数据进行处理后导入到第三方可视化程序里才能观看，特别耗时、繁琐。Multiwfn生成的图像都可以直接导出到图像文件，分辨率可随意设定，支持的格式丰富，包括ps、eps、pdf、wmf、gif、tiff、bmp。



3 Multiwfn的使用方法简介

初学者使用Multiwfn之前强烈建议先阅读一遍《Multiwfn入门tips》（http://hi.baidu.com/sobereva/item/896ee3a19f6d7d3c020a4d76），其中对Multiwfn的基本信息和基本使用方法进行了完整的概括，下面只给出一些关键性的信息。另外在使用程序前务必阅读程序手册的第二页，即“All users must read”中的信息。

Multiwfn的主页是http://multiwfn.codeplex.com，所有Multiwfn版本的源程序、编译好的可执行文件和手册pdf文档皆可从此网站免费下载，下载不需要注册或申请。由于服务器在国外，所以国内用户需确保能访问国外网站。可执行文件的压缩包很小巧，仅不到5MB，下载、传播都很便利。

Windows版Multiwfn支持从WinXP到Win8的各个版本，32/64bit系统皆可兼容。程序不需要安装，只要解压后即可直接使用。源代码可在CVF6.x或IVF12.x下顺利编译。对于初学者强烈建议使用Windows版，Multiwfn对Windows支持是最完美的。
Linux版的Multiwfn支持Red Hat Enterprise Linux 6.x 64bit，对于其它Linux发行版一般来说也都可用，但不保证能完美支持。使用Linux版之前需要根据手册对系统进行配置，并且可能需要自行补全所需的库文件。Linux版Multiwfn源代码可以在ifort 12.x下顺利编译（gfortran等编译器一般无法通过）。

Multiwfn支持各种类型波函数，包括闭壳层/非限制性开壳层/限制性开壳层的HF/DFT波函数，以及闭壳层/开壳层的后HF波函数。支持高斯函数角动量最高到g。Multiwfn对于分析的体系尺度没有限制，原则上只要内存够大，就可以分析任意大的体系。

启动Multiwfn后首先需要载入输入文件。Multiwfn的功能十分丰富，由于不同种类的波函数分析所需要的信息不同，所以要用的输入文件也有所不同。Multiwfn可以载入的文件类型很多，大多数功能要用到.wfn文件，可以由Gaussian, GAMESS-US/UK, Firefly, Q-Chem, deMon2k等程序产生，也可以通过Zork开发的molden2aim程序将Molden输入文件进行转换来得到。由于Gaussian的.fch文件包含的信息远多于.wfn，用.fch作为输入时可以启用Multiwfn里更多的功能。还有的功能则要用到cube/grd文件、Gaussian或NBO程序输出文件等，详见《Multiwfn入门tips》的说明。另外顺带一提的是，在将要发布的3.0.1版中，Molden输入文件将可以直接被Multiwfn所读取，效果与fch文件等同，而Molden输入文件可以被更广泛的量子化学程序所生成，如ORCA、Molpro、deMon2k、BDF等。

Multiwfn的所有操作都是交互式的，大部分功能只要按照屏幕上的提示一步步进行即可，要输入什么在屏幕上都会看到明确的说明。遇到问题时应先参阅手册对功能的介绍，以及手册中相关的例子，并举一反三。手册的内容很详尽，简单易懂，查阅比较方便，章节结构在《Multiwfn入门tips》里有介绍。

虽然Multiwfn操作是交互式的，但也可以方便地通过重定向方式以命令行模式来运行。也可以直接编写shell脚本来实现Multiwfn的批处理。方法在手册附录里有具体介绍。

使用Multiwfn发表文章时应引用分析相应方法的原文（在手册的理论介绍部分里都会给出），并引用Multiwfn的介绍文章Tian Lu, Feiwu Chen, J. Comp. Chem. 33, 580-592 (2012)。如果定量分子表面分析功能被用到，还应同时引用笔者的介绍这一算法在Multiwfn中具体实现的文章Tian Lu, Feiwu Chen, J. Mol. Graph. Model., 38, 314-323 (2012)。

笔者的电子邮箱是Sobereva@sina.com，关于任何使用上的问题或者建议都可以通过这个邮箱与笔者联系。也可以在思想家公社QQ群里咨询，群号码是18616395，加入时应注明是Multiwfn用户。也可以在小木虫的量子化学区的Multiwfn板块进行讨论。如果长时间未得到回复建议直接发电子邮件咨询，并注明操作过程、所用版本，最好将输入文件压缩后作为附件一起发送。



4 Multiwfn的基本功能

4.0 Multiwfn支持的实空间函数

Multiwfn的众多功能都与实空间函数密切相关，实空间函数在这里指分布在三维空间中的标量函数，它们在波函数分析中有重要价值，所以在介绍Multiwfn的功能前先将Multiwfn支持的实空间函数列举一下。如果对它们不熟悉，可以参看手册2.6节的简要介绍。Multiwfn支持计算几乎所有最重要的实空函数，包括：
(1)电子密度
(2)电子密度的梯度
(3)电子密度的拉普拉斯函数
(4)轨道波函数。轨道可以是正则分子轨道、定域化轨道、自然轨道、NTO、NBO、NHO、NAO等，取决于输入文件类型
(5)电子自旋密度或自旋极化参数
(6)电子的哈密顿动能密度K(r)
(7)电子的拉格朗日动能密度G(r)
(8)由原子核产生的静电势。如果用.chg文件（记录了原子电荷数值）作为输入文件，则对应的是原子电荷产生的静电势
(9)电子定域化函数(Electron localization function, ELF)。包括Becke和Tsirelson各自提出的形式。开壳层体系的Becke的ELF所用的具体形式是由笔者导出的，见《电子定域化函数的含义与函数形式》（物理化学学报,27,2786-2792）
(10)定域化轨道指示函数(Localized orbital locator,LOL)。包括Becke和Tsirelson各自提出的形式。关于ELF、LOL和电子密度拉普拉斯函数的简介可参见《电子定域性的图形分析》（http://hi.baidu.com/sobereva/item/4a8aeae2291d8d0f8c3ea8d8）
(11)局部信息熵
(12)总静电势(ESP)
(13)约化密度梯度(RDG)
(14)Promolecular近似下的RDG
(15)Sign(lambda2)*rho（电子密度Hessian矩阵第二大本征值的符号乘以电子密度）
(16)Promolecular近似下的Sign(lambda2)*rho
(17)电子的交换穴、相关穴、相关因子
(18)平均局部离子化能(ALIE)
(19)源函数(Source function)
Multiwfn还支持很多其它的实空间函数，由于不很常用，或者可以由上述函数经过简单运算获得，所以就不一一列举了，比如Weizsacker泛函、局部温度、平均局部静电势、形状函数（Shape function）、势能密度、能量密度、键的金属性（Bond metallicity）、近似形式的DFT线性响应核等。

出于灵活性考虑，Multiwfn提供了一个自定义函数，用户只要在程序源代码的userfunc子程序里面简单地填上代码然后编译，就可以让Multiwfn直接分析和绘制用户自行定义的实空间函数。

下面将Multiwfn的比较主要的功能依次进行简要说明。限于篇幅，有大量功能、选项无法在此一一提及。


4.1 观看分子结构并绘制轨道

载入含有分子坐标的文件（pdb、xyz等），在Multiwfn的主功能0里就能观看分子结构。如果载入的文件包含波函数信息（wfn/wfx/fch/molden等），还可以观看轨道等值面。如果载入的是格点数据文件（.cub或Dmol3的.grd文件），则可以直接观看格点数据的等值面。



特别值得一提的是，Multiwfn显示轨道等值面的速度飞快，远胜于Gaussview产生图像的速度。而且观看轨道特别方便，想看哪个轨道直接在列表里一点就出现，调整等值面对应的数值（isovalue）只需要直接拉动滑条，远比Gaussview省事。Multiwfn还能同时显示两个轨道，分析轨道相互重叠问题的时候很有用。

Multiwfn可以方便地观看NBO程序产生的各种轨道，包括NBO（自然键轨道）、NHO（自然杂化轨道）、NAO（自然原子轨道）、NLMO（自然定域化分子轨道），以及前头带P的类型，如PNBO（初自然键轨道）等，详见《使用Multiwfn绘制NBO及相关轨道》（http://hi.baidu.com/sobereva/item/b0226b3a534a4b637d034bc5）。而其它能观看这些轨道的程序，如NBOView、Chemcraft都是收费的，并且功能和易用性上都不如Multiwfn。

等值面的显示方式支持：不透明表面、网格状表面、点状表面、不透明表面+网格状表面、透明表面。等值面和网格的颜色都可以任意设定。显示分子结构时的原子尺寸、是否成键的判断阈值、标签大小等参数也可以自由调节。


4.2 计算实空间函数在一个点、一条线、一个面上和一个空间区域内的数值并绘图

这几种功能是Multiwfn最初版本就加入的，也是Multiwfn最核心、最强大的功能之一。

在主功能1里，用户只要输入一个点的坐标就能立刻将所有支持的实空间函数在此处的数值，连同指定的函数的梯度、Hessian矩阵一起输出出来。

在主功能3里，用户只要定义一条线的两个端点的坐标，就能迅速将指定的实空间函数在这条连线上的变化作出曲线图，并且还可以给出曲线上的极大极小点位置、特定函数值处的坐标位置。曲线的数据点也可以导出，供用户在第三方程序如Origin里直接作图。



在主功能4里，用户可以对指定的函数在指定的面上作图。支持的图形种类丰富，包括填色图、等值线图、地形图（可以包含填色和投影效果）、梯度线图、向量场图。作图不仅特别方便、速度很快，而且效果很好，有丰富选项可以用于调节作图参数。还可以在图上直接画出分子范德华表面，显示出拓扑分析的临界点、拓扑路径、盆分界线。平面上的数据也可以导出，用于用户自行在sigmaplot、matlab、surfer等程序上绘制。下面是几个Multiwfn直接绘制的平面图的例子，未经任何加工：

两个NBO轨道的等值线图，可以清楚地考察二者的重叠程度


ClF3的静电势等值线图，实、虚线分别代表正、负值。蓝色粗线表示的是Bader定义的分子范德华表面（电子密度=0.001等值面）


尿嘧啶的电子密度梯度场图


镁卟啉的电子密度等值线图+AIM临界点、键径和原子盆间面


B13+团簇的LOL填色图


Li6团簇的ELF填色地形图+投影。（Multiwfn的Logo就是来自Li6的ELF填色图）


在主功能5里，可以计算指定实空间函数的格点数据，并且绘制等值面图。格点的定义非常灵活方便，计算速度也远胜具有相同功能的其它各种程序。格点数据既可以导出到文本文件也可以导出到Gaussian型cube文件中。也可以之后进入Multiwfn的其它功能对刚算完的格点数据进行进一步利用和操作。在主功能5里也可以计算指定一批数据点上的实空间函数，数据点坐标通过外部文本文件导入。

在计算直线、平面上以及空间区域内的实空间函数时，用户可以设定“自定义运算”操作，来让多个波函数文件对应的实空间函数相互运算。比如作密度差图，只需先载入整体的波函数文件，然后设定它与哪些片段的波函数文件相减，然后在计算密度的时候就会自动依次算它们的密度并做差值。因此用Multiwfn计算诸如Fukui函数、双描述符特别方便，详情可参见《使用Multiwfn作电子密度差图》（http://hi.baidu.com/sobereva/item/96eaeb5dd1b12309e7c4a5c1）。当然，不只是密度差，也不光只是减法运算，Multiwfn可以直接让任何支持的实空间函数对多个波函数文件进行加减乘除运算。另外，Multiwfn还可以直接得到promolecular属性和deformation属性。前者就是原子在自由状态属性的叠加，后者是指实际体系的属性与promolecular属性的差值。这里的属性是指任意一种实空间函数，如果选的是电子密度，deformation属性具体就是指电子变形密度，这对于分析原子形成分子时电子密度如何聚集、转移十分有用。

如果产生相同体系不同坐标下的一批波函数文件（如对应IRC过程的每个点），以批处理方式调用Multiwfn对它们依次产生图像，并结合做动画的程序，就可以得到描述电子结构变化的动态图像，对于分析反应过程内在机制很有帮助，详见《制作动画分析电子结构特征》（http://hi.baidu.com/sobereva/item/c69893149c743a6e70d5e8c3）。


4.3 拓扑分析

在AIM分析中，通常要对电子密度进行拓扑分析，包括寻找各种临界点、产生拓扑路径。Multiwfn不仅可以对电子密度做拓扑分析，还可以对ELF、LOL、电子密度拉普拉斯函数等实空间函数也进行拓扑分析。搜索临界点和产生拓扑路径速度非常快，显著胜于其它有类似功能的程序。结果可以直接在图形界面里可视化（如下图所示），十分方便。还可以绘制出盆间分界面(interbasin surface)。每个临界点上各种实空间函数的数值都能直接得到。临界点和拓扑路径可以灵活地自行修改、导入导出。Multiwfn还提供了测量工具来获得临界点/原子核间的几何参数。关于Multiwfn的拓扑分析功能的介绍具体可参见《使用Multiwfn做拓扑分析以及计算孤对电子角度》（http://hi.baidu.com/sobereva/item/c84fa5089156806dd55a11c3）。




4.4 检查和修改波函数

Multiwfn提供了许多选项对波函数进行检查和修改。包括将当前波函数输出为新的wfn文件；输出GTF、基函数和轨道的信息；输出系数矩阵、重叠矩阵、密度矩阵；交换、设定GTF的信息（所属中心、系数、指数等）；设定轨道的占据数和轨道类型；平移或平移+复制体系（后者可以用来近似地将单胞波函数扩展为超胞波函数）；将某些原子从波函数信息中去除等。

对波函数进行修改后会直接影响实空间函数的计算结果，这对于一些特殊目的的分析十分有用，比如可以只考察某些原子、某些轨道（如pi轨道）上的电子产生的实空间函数。


4.5 布居分析

Multiwfn提供了丰富的布居分析功能，用于计算原子电荷或基函数/壳层上的电子占据数。支持的方法包括Hirshfeld布居、Voronoi变形密度布居(VDD)、Mulliken布居、Lowdin布居、修改的Mulliken布居（包括三种方法：SCPA、Stout-Politzer和Bickelhaupt）、Becke布居、原子偶极矩校正的Hirshfeld布居(ADCH)、CHELPG和Merz-Kollmann(MK)拟合静电势方法。

Mulliken布居分析功能特别灵活，不仅可以得到基函数/壳层/原子上的电子布居数（开壳层体系会同时给出自旋布居），还可以得到原子间、基函数间的重叠布居以及它们自身的局部布居，并且结果还可以再进一步分解为每个轨道的贡献。

虽然主流的程序如Gaussian也能做CHELPG和Merz-Kollmann(MK)拟合静电势计算得到原子电荷，但是Multiwfn的这两个功能更为灵活，用户可以任意设定拟合点的分布，还可以任意设定额外的拟合中心（为了准确描述分子间静电相互作用，有时必须引入额外的非原子中心的点电荷）。

上述功能中ADCH布居方法是由笔者在J. Theor. Comput. Chem., 11, 163 (2012)中提出的。笔者在《原子电荷计算方法的对比》（物理化学学报,28,1-18）中对许多种计算原子电荷方法相比较，其中ADCH展现出了优秀的性质，特别是能精确重现分子偶极矩，同时也有很好的静电势重现性。Multiwfn是唯一能实现ADCH方法的程序。


4.6 轨道成份分析

Multiwfn提供了多种方法进行轨道成份分析来得到各个原子轨道、壳层和原子对轨道产生的贡献，并可以灵活地定义片段。支持的方法包括Mulliken、Stout-Politzer、SCPA、NAOMO和Hirshfeld。笔者写的《分子轨道成分的计算》（物理化学学报,69,2393）是迄今第一个专门讨论轨道成份分析方法的文章，其中对这些方法的原理进行了介绍并通过实例做了对比分析。NAOMO和Hirshfeld方法做轨道成份分析也是在这篇文章里首次提出的，并且目前只有Multiwfn能够实现这两种方法（注：NAOMO分析需利用NBO程序的输出信息）。另外在《谈谈轨道成份的计算方法》（http://hi.baidu.com/sobereva/item/d52876160110ca8d89a956c5）当中也对这些方法进行了粗略介绍，并给出了在Multiwfn中操作的实例。


4.7 键级分析

Multiwfn支持的键级的计算方法众多，包括：
(1)Mayer键级。这是目前最常用的键级。
(2)轨道占据扰动的Mayer键级。此功能用于分析各个轨道的电子对Mayer键级的影响（扰动）。
(3)Wiberg键级。此键级原始定义只能用于正交基，由于目前主流量化程序使用非正交基，Multiwfn会自动先做Lowdin正交化。
(4)多中心键级。这对于分析多中心键、芳香性极有用。Multiwfn最多支持到10中心，且可以自动搜索多中心键。
(5)Mulliken键级。也称Mulliken重叠布居数。Multiwfn可将它分解为轨道的贡献。
(6)模糊键级（Fuzzy bond order）。相当于模糊空间下计算的离域化指数。对于离子性不很强的键其结果和Mayer键级结果相仿佛，虽然计算会多耗时一些但是对基组的敏感性显著降低。
(7)Laplacian键级(Laplacian bond order, LBO)。这是笔者在J. Phys. Chem. A, 117, 3100 (2013)中提出的，基于对成键区域（模糊重叠空间）内的电子密度Laplacian的负值区域积分来获得。LBO与键的强度有极好的对应关系，并且能清楚地反映出键的极性。

对于上述部分键级，Multiwfn还可以定义片段，来直接得到两个片段的原子间的键级总和。


4.8 绘制态密度图

虽然态密度(Density-of-states, DOS)一般是对于周期性体系来讨论的，但是对于分子体系，将能级进行适当地展宽得到态密度，对于讨论电子结构是很有益的。Multiwfn可以产生Total DOS (TDOS)、Partial (PDOS)、Overlap population DOS (OPDOS)。通过TDOS可以直观地考察轨道在各个区域能量分布状况，通过PDOS图可以直观地考察指定的片段对各个能量区域内的轨道的贡献，而借住OPDOS图可以直观地了解不同能量区域的轨道对两个片段间的结合能力（以Mulliken键级来表征）产生的影响。

对于PDOS和OPDOS，Multiwfn提供了十分灵活方便的界面用于定义片段。片段可以由基函数、壳层、原子以任意方式组合而成。Multiwfn也提供了许多选项用于用户自行调节作图设定，诸如展宽函数、半高宽(FWHM)等。

下面是Multiwfn产生的二茂铁的TDOS+PDOS+OPDOS图（做了些许ps以使内容更清晰）




4.9 绘制光谱图

虽然绘制光谱图不属于波函数分析，但是是理论化学研究中常涉及到的，所以Multiwfn也具备这样的功能，可以绘制IR/Raman/UV-Vis/ECD/VCD图。激发能和强度数据可以从Gaussian的输出文件中直接读取，也可以通过文本文档导入。作图的可调参数十分丰富，可以满足专业用户的需求（比诸如Swizard更灵活更好用）。




4.10 定量分子表面分析

这里的“分子表面”是指电子密度等值面，比如0.001 e/bohr^3的表面常被用来作为分子的范德华表面，这是最初由Bader定义的。“定量分析”就是指的在这样的表面上，对映射到上面的实空间函数进行定量的分析，诸如统计平均值、方差、不同数值范围所占的面积大小，并且得到所映射的函数在分子表面上的极小值和极大值。虽然原则上各种实空间函数都可以作为映射到分子表面的函数，但通常只将静电势和平均局部离子化能(ALIE)映射上去。对静电势的定量分子表面分析对于预测分子间相互作用、通过QSPR（定量结构属性关系）方式预测分子凝聚态性质十分有用，而对ALIE的分析则是一种很重要且很成功的预测亲电反应位点的方法。例如下图是Multiwfn产生的ALIE在苯酚分子表面上的极值点分布，极小点（蓝球）清楚地表现了亲电位点是在邻对位。



定量分子表面分析的详细介绍和一些应用实例可参见《使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用》（http://hi.baidu.com/sobereva/item/6601fff28c4fd9d643c36acc），算法细节见J. Mol. Graph. Model., 38, 314。


4.11 处理格点数据

Multiwfn的处理格点数据的功能源自于笔者早前开发的GsGrid程序，但它的开发早已终止，因为GsGrid的全部功能已被整合进了Multiwfn，同时做了不少扩充和完善。使用Multiwfn的这个功能前需要先载入cube或Dmol3的.grd文件，它们记录了格点数据；也可以先用Multiwfn自身计算格点数据的功能产生格点数据，然后进入这个功能直接对其进行处理。

Multiwfn的处理格点数据功能包含众多子功能：
(1)观看格点数据等值面。
(2)将当前格点数据导出为cube文件。
(3)将一些格点的坐标和数值导出到文本文件。决定格点范围的标准可以为：(a)所有格点 (b)接近某个平面的格点（特定的X/Y/Z平面或通过三个点来定义） (c)数值在指定范围内的格点。
(4)设定远离某些原子的格点的数值，可以由此屏蔽不感兴趣的区域的等值面；设定远离两个片段间重叠区域的格点的数值，可以由此只保留重叠区域的等值面，对于RDG方法靠等值面展现弱相互作用区域比较有用，见手册4.13.4节的例子。
(5)对当前格点数据进行运算（例如取绝对值、求对数、求幂），或者从cube/grd文件中载入另一套格点数据与当前格点数据相运算（相互加减乘除、求二者的平方和、取平均等等），由此产生新的格点数据。
(6)从cube文件中载入另一套格点数据并映射到当前格点数据中具有指定数值的格点上，并将坐标和数值输出为文本文件。
(7)对格点数据范围进行scale，以使得格点数据的范围落在指定区间内。
(8)对指定数值范围或空间范围的格点数据进行统计，给出最大/最小值、平均值、均方根、总和、积分值、标准偏差、正值/负值/整体区域内数值的重心。
(9)绘制X/Y/Z方向的格点数据的积分曲线。如果当前的格点数据是片段密度差，那么这积分曲线也叫电荷转移曲线，对于分析电荷转移很有益，如JACS,130,1048。


4.12 自适应自然密度划分(Adaptive natural density partitioning, AdNDP)

AdNDP在某种意义上是对NBO方法的扩展，可以得到“半定域化”轨道，它既不像分子轨道那样一般会离域到整体，又不像NBO或传统的定域化轨道那样只分布在1~3个原子上，而是可以显示出局部区域的离域特征，特别适合分析、展现多中心键，在研究团簇芳香性上已有大量应用。

Multiwfn是目前唯一一个公开的AdNDP分析程序。虽然AdNDP方法本身有一些人为性而不得不由用户来干预AdNDP轨道的搜索，但是Multiwfn精心设计的界面已经将这个过程最大程度地便利化，并且可以直接观看AdNDP轨道图形。例如下图是菲的边缘的苯环上的三个六中心AdNDP轨道，展现了两侧的苯环的芳香性。



关于AdNDP方法在理论方面的介绍和讨论，以及在Multiwfn中的实际操作，在《使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键》（http://hi.baidu.com/sobereva/item/02b67c3ca3aa46c7382ffac6）当中有详述。


4.13 模糊空间分析

对体系中原子所属空间的划分有两类，一类是离散式划分，即原子间有明确的边界，例如AIM的原子盆就是典型的一种；另一类是模糊式划分，原子间没有明确的边界，而用平滑的权重函数来表示原子空间，越接近原子核则此原子的权重越大。模糊空间定义有很多，如Becke和Hirshfeld所定义的。在模糊空间下分析的好处之一是积分原子空间来得到原子的性质比较容易，这是由于边界平滑所以容易积分得比较准确。Multiwfn的模糊空间分析就是用来计算各个模糊原子空间里的性质，来展现原子在分子中的特征。具体可以做以下分析：
(1)计算模糊空间以及模糊重叠空间内的任意实空间函数的积分值（模糊重叠空间是指两个模糊原子空间之间重叠的空间）
(2)计算模糊空间内的电多极矩（即得到原子单极矩、原子偶极矩、原子四极矩）
(3)计算模糊空间内各个分子轨道间的重叠积分矩阵（也叫Atomic overlap matrix）
(4)基于模糊空间计算定域化指数(Localization index)、离域化指数(Delocalization index)和简缩线性响应核(condensed linear response kernel)
(5)计算对位离域化指数(Para-delocalization index, PDI)、芳香波动指数(Aromatic fluctuation index, FLU)和对位线性响应指数(Para linear response index, PLR)用于研究芳香性
(6)计算多中心离域化指数用于研究多中心键


4.14 电荷分解分析(Charge decomposition analysis, CDA)

CDA用来将体系内两个片段之间的电荷转移进行分解，可以获知各个分子轨道（或自然轨道）的形成导致电荷在片段间以什么方向转移、转移的量是多少，而非仅仅只得到一个整体转移量。CDA对于分析金属-配体间的电荷转移极为有用。CDA也同样可以分析共价键相连的片段间的电荷转移。

Multiwfn还支持扩展的CDA (Extended CDA)。CDA没有明确将电荷转移和电荷极化分开考虑，而ECDA则将它们分离开，排除了电荷极化效应对计算电荷转移量的影响。

Multiwfn的CDA模块还可以绘制轨道相互作用图（如下所示）、给出各个片段轨道在复合物轨道中所占成份、输出片段轨道之间的重叠矩阵。



关于CDA的详细讨论以及在Multiwfn中的使用，见《使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图》（http://hi.baidu.com/sobereva/item/178fbd02969cb690a2df4375）


4.15 盆分析

Multiwfn的盆分析功能极其灵活、方便、强大和高效，所用的数值算法是基于立方格点的算法（见J. Phys.: Condens. Matter,21,084204）。不仅可以对电子密度寻找吸引子（也叫(3,-3)临界点，即局部极大点）并产生相对应的盆，还可以对所有其它的实空间函数寻找吸引子并产生盆，诸如ELF、LOL、静电势，甚至于密度差、福井函数等。每个盆都是体系中的一个局部空间，比如电子密度的盆对应于AIM原子空间，而ELF和LOL的每个盆则对应于一个电子结构的特征区域，如共价键区域、孤对电子区域、内核区域等。在产生盆之后，就可以直接观看吸引子的位置和盆的范围，例如下面的绿色区域代表了HCN中H的电子密度的盆：



Multiwfn可以对盆进行以下分析：
(1)在盆内积分任意实空间函数。对于在原子核附近波动大的函数如源函数，可以使用动态调节积分格点精细度的积分方法以保证精度。
(2)计算盆内的电单极矩、偶极矩和四极矩。
(3)计算定域化和离域化指数并输出盆区域内的轨道间重叠矩阵。

Multiwfn出于灵活、便利考虑，还提供其它一些选项，诸如测量吸引子/原子核间的几何参数、将盆导出到cube文件并将吸引子导出到.pdb文件便于用外部程序作图。

关于盆分析的基本理论、算法和Multiwfn中的分析实例见《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》（http://hi.baidu.com/sobereva/item/47427a1487ce95633e87ce25）。


4.16 其它实用功能

以上介绍的功能属于Multiwfn的主功能。在Multiwfn当中还包括很多计算化学中会用到的实用工具以及与波函数分析有一定关联但不是“大型”的实用功能：

(1)同时计算两个实空间函数的格点数据并绘制它们之间的散点图。这对于基于RDG的弱相互作用分析特别有用。
(2)基于蒙特卡罗方法计算分子范德华体积。此功能在《谈谈分子体积的计算》（http://hi.baidu.com/sobereva/item/2707923863a0fd5b80f1a7df）当中对原理和操作有介绍。
(3)对任意实空间函数在整个空间范围内进行积分。积分方法用的是Becke提出的对DFT交换相关泛函进行积分的方法，积分可以达到相当高的精度。结合Multiwfn提供的自定义函数，此功能十分灵活有用。比如自己提出一个新的相关泛函，那么在源代码的自定义函数中填上自己的泛函的公式（计算密度、密度的梯度、动能密度在Multiwfn里都有现成的函数可以直接调用），利用此功能对自定义函数进行积分就能得到这个相关泛函所对应的相关能。
(4)计算非限制性计算的Alpha和Beta轨道间的重叠积分，这对于研究自旋极化程度是有用的。
(5)监视Gaussian的SCF收敛过程并将收敛过程绘图。这有助于了解收敛状态，判断是否应该停掉任务重设关键词。
(6)产生带有片段初猜的Gaussian输入文件。这对于产生高质量初猜波函数、计算自旋极化单重态体系很有用，基于此还可以进行比较简单的能量分解，可以将“极化能”与“静电能+交换能”二者从总相互作用能中分离开。
(7)分析电荷转移。这个方法主要用于激发态过程的电荷转移分析，可以很清楚地考察电荷转移距离、转移电荷的分离程度等信息。详见J. Chem. Theory Comput., 7, 2498 (2011)及手册相关章节。
(8)计算两个轨道的模的重叠积分。
(9)绘制Gaussian输出的跃迁密度矩阵，这可以图形化直观地展现出电子跃迁涉及到的区域和程度以及原子间的相干性，这对于分析大尺寸、共轭分子的跃迁十分有用。详见《绘制跃迁密度矩阵平面图分析电子跃迁》（http://hi.baidu.com/sobereva/item/fe16b92f0c57eaf951fd87c6）
(10)计算HOMA和Bird芳香性指数。HOMA是非常流行的基于结构的度量芳香性的指标。
(11)计算LOLIPOP (Localized Orbital Locator Integrated Pi Over Plane)。这是一种基于对芳香环上方的LOL值进行积分来衡量芳香体系pi-pi堆积能力的指标，详见Chem.Commun.,48,9239。
(12)计算分子间轨道重叠积分。这个积分与分子间电荷转移能力密切相关。详见《分子间轨道重叠的图形显示和计算》（http://hi.baidu.com/sobereva/item/7cc644e165d9b10c65db00ae）。
(13)Yoshizawa电子转移路径分析。此功能基于Yoshizawa公式分析电子在平面共轭体系中可能的转移路径，原理见Acc. Chem. Res., 45, 1612 (2012)。
(14)对平面体系自动检测出pi轨道及其它轨道，批量设定它们的占据数。占据数设为0相当于此轨道在计算中被忽略，由此可以在各种Multiwfn的分析中方便地去除pi电子或其它电子的贡献。
(15)将函数分布拟合为原子的数值。这与MK拟合静电势方法很相似，但是被拟合函数可以是任何其它实空间函数。
(16)辅助计算非平面体系的NICS_ZZ。虽然NICS_ZZ对于平面体系的计算很容易，但是对于非平面体系不易得到要计算的点的坐标，也难以得到垂直于环平面方向上的磁屏蔽值的分量，而Multiwfn的这个功能可以帮助解决这些问题。

Multiwfn还有不少其它辅助功能，包括将当前结构保存为pdb及Gaussian输入文件、计算原子配位数、计算指定一批原子的几何中心、质心、转动惯量、主轴、回转半径等信息。



5 Multiwfn的用途

Multiwfn功能很多、用法灵活，所以有非常广泛的实际用途，不可能一次性全面概括，这一节仅列举一些比较常见的应用，使读者了解前述众多的功能能发挥什么实际用处。此节所讨论的分析方法都可以直接靠Multiwfn实现，细节在手册中都有介绍，重要的应用在手册或相关帖子里也都有实例，请自行查阅。很多文献里的研究方法虽然在本文没有提到，但是只要对Multiwfn的原理和操作领悟透彻，仅需稍作变通，自行写一点处理脚本或稍微修改下Multiwfn的代码就能轻易地实现。

5.1 观看轨道

这已经在4.1节提到过了。Multiwfn观看轨道又方便速度又快，还能支持NBO程序产生的轨道。不仅可以绘制轨道的等值面图，还可以方便地作轨道波函数的各种类型的平面图、曲线图。

作等值面图时如果想得到更好的显示效果，还可以生成轨道波函数的cube文件，用第三方软件如VMD来观看。而平面图、曲线图则允许将数据导出，可放到诸如sigmaplot、origin、surfer等程序里重新作图并更灵活细致地调节作图选项。


5.2 分析化学键

化学键的分析是Multiwfn的主要用处之一。首先可以考虑图形分析。例如可以直接作ELF、LOL图，考察电子主要定域的区域，由此判断键的极性与类型、指出共享电子与孤对电子存在区域等。也可以作电子密度拉普拉斯图，分析电子聚集区域。还可以做变形密度图（分子的密度与组成它的各个原子在自由状态密度之差），展现原子成键后电子分布如何变化，比如通常电子密度会在共价键形成区域有显著增加。简要讨论可参见前述的《电子定域性的图形分析》，更多的讨论可参见笔者整理的“ELF综述和重要文献合集”里的文章（下载地址http://pan.baidu.com/share/link?shareid=465261&uk=1074012119）。ELF/LOL/拉普拉斯/变形密度分析往往有很大的相互关联，可以得到类似结论，即原子间电子定域性高的区域往往和电子凝聚区域相对应(详见JPC,100,15398)，也往往是变形密度图中密度显著增加的区域；但由于不同函数的物理本质有所差异，不同函数着重展现不同的侧面，所以可以放到一起共同来讨论，取长补短。能量密度的分析对于判断化学键的类型（是闭壳层作用还是共价作用）也往往比较有用，成键区域能量密度为负往往对应于共价作用。可参见《AIM键临界点处电子密度拉普拉斯值符号判断相互作用类型失败原因的图形分析》（http://hi.baidu.com/sobereva/item/65c98f0c122375036c9048ab）文中的一些讨论。

拓扑分析和盆分析可以使得实空间函数的分析定量化、严格化。例如AIM理论认为键临界点(BCP)是表征化学键特征的关键的点，所以BCP上的电子密度、拉普拉斯值、能量密度等性质在大量文章中被用来探讨化学键强度和本质问题，如Inorg.Chem.,39,2360、The Quantum Theory of Atoms in Molecules一书的13章等。键径的轨迹（键径不总是直的，可以是弯曲的）被用来讨论键的张力等问题，而键径的存则被一些人当做原子间有化学键存在的必要标准。对于ELF及LOL的拓扑分析也常被用来研究化学键的特征，如Nature,371,683、JCC,30,1093。ELF/LOL的盆与电子结构的特征区域相对应，积分与化学键相对应的ELF/LOL盆是研究化学键特征的重要途径之一，例如可以获知有多少电子涉及成键、定域程度如何、成键区域内电子密度的极性等，前述的《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》文中有实例，在上述的“ELF综述和重要文献合集”里也有很多文章对此有探讨。

键级分析对研究化学键极为有用。例如Mayer键级是很常用的键级度量标准，Dalton Trans.,2001,2095给出了大量应用实例。模糊键级则更适合存在弥散函数的情况。Mulliken键级虽然在定量上不好，但它可以精确分离为轨道贡献之和，其中贡献为负值的轨道体现出此轨道的电子对于相应原子间成键存在削弱效果。另外可以作OPDOS图来直观地考察不同能量范围的轨道对于指定两个片段间成键（以Mulliken键级来表征）的影响。拉普拉斯键级(LBO)则既能展现键的强度又和键的极性相密切关联，在讨论电子密度拉普拉斯函数时可以借助LBO将讨论定量化。多中心键级在探讨多中心键的强度时极为有用，利用AdNDP方法则可以将多中心轨道可视化。

CDA对于研究原子间电荷转移的细节很有用，诸如反馈pi键之类的概念靠CDA都可以很清楚地展现，可参见JPC,99,9352等文章和前述的《使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图》。


5.3 研究化学反应过程

上述分析化学键的方法如果用到分析化学反应过程上，可以显著加深对化学反应过程本质的认识。具体来说，就是取IRC路径上的一些点，都绘制ELF、LOL、密度差图等实空间函数图形，讨论这些函数随反应过程的变化。也可以把IRC上所有点对应的函数图形都作出来组合成动画，见前述的《制作动画分析电子结构特征》文中的例子。

同时也很建议将IRC上各个点的键级计算出来，并绘制出曲线图（纵坐标是键级，横坐标是IRC路径上的坐标，同时把能量变化曲线也画在上面），其中要考虑的键是反应过程中要断裂的键和要形成的键。从图上将可以清楚地看出随着反应的进行，要断的键的键级逐渐趋近于0，而新的键的键级逐渐增加，而过渡态的位置往往就在这两套键级相交的位置附近。很多反应过程在过渡态附近都会出现多中心键特征，比如DA加成反应过程中就有一定的六中心键特征，对这种情况很建议也把多中心键级数值附在曲线图上进行讨论。另外还可以把诸如原子电荷数值、盆积分数值等性质在IRC上的变化曲线也绘制出来加以讨论。


5.4 弱相互作用分析

弱相互作用指的是诸如范德华作用、氢键、卤键、二氢键、pi-pi堆叠等形式的作用，相较化学键要明显更弱，但是却和化学键有着等同的重要意义。Multiwfn在分析弱相互作用上有很重要的用处。杨伟涛课题组在JACS,132,6498中提出的基于约化密度梯度(RDG)的分析方法是一种十分直观有效的表征相互作用区域和相互作用类型的方法，从小分子、固体到生物体系几乎无不适用，虽然提出时间不长但已经得到了很广泛的应用。其原理、在Multiwfn中的操作以及在VMD中绘图的过程在《使用Multiwfn图形化研究弱相互作用》（http://hi.baidu.com/sobereva/item/bb1162d881b8d5ee3cc2cbdc）中有详细介绍。最近还有文章主张可以将这种RDG分析方法和ELF方法结合起来用，同时显示二者的等值面，从而同时展现化学反应过程中强相互作用和弱相互作用的变化，见JCTC,8,3993。

氢键、卤键、二氢键、pi-氢键的主要成份是静电效应（虽然色散作用也参与其中），也就是显正电的区域与显负电的区域相作用，用化学的语言来说则是Lewis酸性区域与Lewis碱性区域的作用。由于静电势与静电相互作用密切相关，因此可以通过定量分子表面分析方法分析单体范德华表面的静电势，特别是静电势极大点和极小点的位置和及数值，来考察这些弱键在什么位置容易以什么朝向形成，也可定性地估计强度的大小。Peter Politzer等人靠这种分析讨论了大量氢键、卤键体系，如JMM,13,305、IJQC,107,3046、PCCP,12,7748。另外，有很多其它种类的复合物的形成也是静电作用所主导的，至少静电作用可以决定复合物形成的朝向，如H3N--GeF4、FH2As--AsH2F等，也同样可以通过分子表面的静电势分析来探讨，见JMM,14,659、JMM,18,541、JMM,15,723等文章。由于电子密度拉普拉斯函数可以表明哪里电子凝聚和疏散，也因此对应哪里易于亲核和亲电，因此对于讨论以静电主导的弱相互作用也是颇有用的，例如JMM,19,2035通过分析拉普拉斯函数的临界点性质很好地解释了几种卤键复合物的相互作用能的大小。

对于静电主导的弱相互作用，也可以计算原子电荷来分析。显然一个电荷越正的原子与一个电荷越负的原子间形成的静电相互作用越强。但是原子电荷只是原子附近电荷分布最简单、抽象的描述，如果想要研究得更为透彻，建议讨论与这些弱相互作用有关的原子的电多极矩（一般考虑原子偶极矩和原子四极矩），它是原子附近电子密度分布的更细致的定量化描述。分析例子可参考Bader的著作Atoms in molecules-A quantum theory的7.4.3节。在Multiwfn中既可以用模糊空间分析模块计算原子多极矩，也可以靠盆分析模块计算原子多极矩。

对于很短很强的氢键（如PNAS,95,12799），由于共价成分不再可以忽视，此时前述的分析化学键的方法，如键级、ELF图形等分析方式也同样适用。

复合物的形成会或多或少地导致电荷的转移（特别是以配位键的方式形成复合物时），可以通过计算原子电荷的变化值检验电荷是如何转移的，也可以让复合物的密度与单体的密度相减得到片段密度差来考察电子的转移方式，做法见前述的《使用Multiwfn作电子密度差图》。Multiwfn还提供积分密度差绘制成电荷转移曲线的功能，见前文4.11节的(8)，由此可以既直观又定量地分析不同X/Y/Z坐标区间内电荷的转移量，尤为适合分析分子与一个平面分子或与固体表面相结合导致的电荷转移。另外还可以用Multiwfn的盆分析功能积分密度差图中的特征区域，由此讨论区域内转移电子量。

基于弱相互作用原子间的BCP的性质的AIM分析是传统的分析弱相互作用的方式。例如对于氢键，Lipkowski在JPCA,108,10865提出了一种常用的判断氢键存在与否的标准，即BCP处密度应当在0.002~0.04之间，电子密度拉普拉斯值应该在0.02~0.15之间。而在一些文章中，如CPL,285,170，给出了BCP处动能密度、势能密度与氢键键能的关系式。近年来也有人利用BCP处的源函数来讨论氢键，如Struct Bond,146,53。

分析范德华表面穿透距离也是分析弱相互作用强度的一种很有用的方法。假设A、B原子核距离分子范德华表面（用电子密度0.001等值面表示）的最近距离分别为rA和rB，那么两原子相互作用越强，显然它们的间距比rA+rB小得就越多，这个偏差值就是范德华表面的穿透距离。在使用Multiwfn定量表面分析功能分析之后，就能看到选项用于输出指定原子核与范德华表面最近距离。

如果研究的是平面共轭体系，那么根据Chem.Commun.,48,9239提出的LOLIPOP指数，可以预测分子pi-pi堆叠能力的大小。


5.5 电荷分布分析

最简单的表现体系中电荷分布的方式就是计算原子电荷，Multiwfn支持很多计算原子电荷的方法，各有特点，可参见物理化学学报,28,1的对比分析。对于Mulliken方法，Multiwfn还可以分析不同类型原子轨道上的电子占据数。原子电荷只能体现原子附近有多少电子，而更细致的研究就是计算前述的原子偶极矩、四极矩，前者可以展现出原子附近电子密度的极化方向和程度，后者可以展现出原子附近密度偏离球对称的方式和程度。

通过盆分析功能得到ELF、LOL函数的盆，则可以将整个空间划分为对应于不同电子特征的局部区域，如内核区域、孤对电子区域、共价键区域等，再对盆进行积分，就可以得到这些区域内的各种性质，如盆内的电子数、偶极矩、四极矩等等。这对于深入分析化学键和弱相互作用的特征也都是很有益的。特别是对应于成键的ELF盆内的偶极矩可以视为键偶极矩，在不少文章中都有讨论）。可以参看前述的《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》的例子。另外还可以讨论诸如体系不同区域的电子对分子总偶极矩的贡献等问题，可参考JCC,29,1440文中的讨论。

直接绘制电子密度图、电子密度拉普拉斯图显然也是表征电荷分布的直观的方法。

分析电子分布的方法也完全可以用来单独分析alpha电子、beta电子或者自旋密度的分布。比如可以直接作自旋密度图、通过Mulliken分析得到特定原子轨道上的自旋布居数、做盆分析/模糊空间分析得到特定区域内自旋电子数。这对于研究磁性体系、自由基体系是很有用处的。


5.6 反应性、反应位点预测

预测反应位点的方法很多，比较重要的包括分析前线轨道成份、计算原子电荷、考察分子范德华表面上静电势以及平均局部离子化能的极值点分布、计算福井函数、计算双描述符，还有人通过分析电子密度拉普拉斯函数的特征来讨论被亲核、亲电进攻的位点和难易程度，如JPC,95,4698。这些分析方法在Multiwfn里全部能实现。有关的原理和操作可参见前述的《谈谈轨道成份的计算方法》、《使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用》、手册的4.5.4及4.7.3等章节中的例子。


5.7 芳香性分析

Multiwfn支持的用来芳香性的方法多达十几种，如AdNDP、ELF-sigma/pi、FLU/FLU-pi、PDI、多中心键级、HOMA等等。原理与操作方法的介绍见《衡量芳香性的方法以及在Multiwfn中的计算》（http://hi.baidu.com/sobereva/item/0a4fb6fb77c4a648922af224）。另外，Multiwfn还提供了辅助计算非平面体系的NICS_ZZ的功能，见手册3.100.24节。


5.8 电子跃迁分析

Multiwfn对于分析电子跃迁有很大帮助。电子激发可以通过分子轨道或自然跃迁轨道(NTO)的跃迁来表示，后者在《跃迁密度分析方法-自然跃迁轨道(NTO)简介》（http://hi.baidu.com/sobereva/item/16a9e1d833658019d68ed0c3）一文中有介绍。利用Multiwfn分析主要涉及的轨道的成份，就可以清楚地讨论跃迁主要涉及到哪些原子，哪些原子的电子被激发到哪些原子上。Multiwfn可以绘制出跃迁密度矩阵图，可以将跃迁过程涉及的原子范围以及跃迁过程中原子间的相干性直观地表现出来，见前述的《绘制跃迁密度矩阵平面图分析电子跃迁》。Multiwfn可以方便地对任意的态之间做密度差图，因此可以对激发态和基态作密度差来图形化研究激发过程中电子是怎么转移的。同时，Multiwfn的主功能100里的功能10是专门用于分析激发过程的电子转移的模块，原文见JCTC,7,2498，实例见手册4.100.4节。它基于密度差数据，经过一番变换后可以将激发过程的电子转移以更简明的方式用图形化表现出来，电子转移距离、转移电子分布的广度、转移电子的分离程度等特征还可以被定量化地表征出来。另外，也可以计算基态和激发态的原子电荷并求差值，以得到不同范围内的电荷转移的确切的量。


5.9 分子性质预测

Multiwfn可以计算分子体积、分子面积（由电子密度等值面定义）、分子表面上静电势的最大和最小值/平均值/方差/电荷分离程度/电荷平衡程度等分子描述符，利用它们来构建QSPR关系可以很好地预测分子凝聚相性质。如JPCA,110,1005就构建出一套关系可以成功地用来预测蒸发焓，即：delta_H_vap=a*sqrt(S)+b*sqrt(v*sigma2)+c，其中S是分子表面积，v是电荷平衡程度，sigma2是分子表面静电势的方差，a/b/c都是拟合得到的参数。


5.10 导电分析

如4.8节所述，Multiwfn可以通过基于对分子轨道能级进行展宽来得到TDOS，还能得到PDOS和OPDOS。从TDOS图可以考察不同能量区域轨道分布的疏密，特别是Fermi能级位置，如果TDOS曲线有明显的凹陷则意味着有相对较大的gap，导电性也相对较弱。通过PDOS可以考察TDOS的组成，在前线轨道区域有较大PDOS的片段内的原子与导电/电子激发有较密切的关系。

对于分子晶体或寡聚物体系，体系的导电能力和单体间电荷转移积分有密切关系，而电荷转移积分又与单体间的轨道重叠积分有很强的相关性。Multiwfn可以利用Gaussian输出文件计算出单体轨道间重叠积分值。而且利用Multiwfn还可以以一种新颖的图形化的方式考察单体间的交叠程度和交叠方式，重叠积分数值大小的成因可以一目了然地搞清楚，见前述的《分子间轨道重叠的图形显示和计算》。

Multiwfn的主功能100里的功能18可以帮助使用者基于Yoshizawa的公式(Acc.Chem.Res.,45,1612)获知共轭体系内电子易于在哪些位点间转移，见手册4.100.9节的例子。


5.11 绘制光谱图

如4.9节所述，Multiwfn可以基于量子化学计算数据产生可以与实验图谱相对比的IR/Raman/UV-Vis/ECD/VCD图。虽然很多计算化学可视化程序也附带了这一功能，但Multiwfn提供的作图参数的灵活程度是绝大多数程序不能及的。


5.12 结构化学教学

给本科生讲述结构化学的时候，将问题用直观的图形来展现远比用复杂抽象的数学公式来讨论更容易让学生们理解、激发学习兴趣。目前的结构化学的教学大纲普遍十分陈旧，笔者建议适当将简单易懂又有重要实际意义的波函数分析理论适当引入课程。

Multiwfn在结构化学教学上有很多用处。例如讲述原子壳层结构、讨论原子轨道时，就可以利用Multiwfn方便地作出相应的图形展现它们的特征。这个应用在《使用Multiwfn绘制原子轨道图形、研究原子壳层结构及相对论效应的影响》（http://hi.baidu.com/sobereva/item/4b3b0f083d82e7103b53eec8）一文中有详述。讨论前线轨道的时候也可以用Multiwfn作出各种形式的图来。

在讨论共价键/离子键，非极性键/极性键的时候，就可以向学生们展示电子密度拉普拉斯图或密度差图，其物理意义对于量子化学外行人来说也很容易理解。在讲述VSEPR模型的时候，可以展示电子密度拉普拉斯函数的(3,+3)临界点位置（即电子凝聚程度局部最大的位置），临界点位置和VSEPR模型描述的杂化方式有很好对应性。讲述多中心键时建议结合多中心键级、AdNDP轨道图以及ELF图形来讨论。在讲述金属配合物电子结构的时候，建议将CDA分析结果、CDA模块绘制出的轨道相互作用图连同轨道图形展现出来。在讲述氢键、范德华作用的时候，可以将上述RDG分析方法给出的图形展现出来，既漂亮又易懂。

涉及到化学反应过程时，将前述的IRC路径上的键级变化曲线，以及反应过程中ELF或密度差图变化的动画向学生们展示，必然会使学生们豁然开朗。

笔者打算以后专门在《大学化学》上写一篇文章讨论Multiwfn在化学教学上的应用。


5.13 其它

Multiwfn还有很多很多其它的应用，比如作为fch或molden输入文件->wfn/pdb/Gaussian输入文件的格式转换器、编辑波函数文件、监控SCF收敛、做简单的能量分解、作为格点文件统计/处理以及相互运算的工具（诸如Gaussian的cubman程序完全可以被Multiwfn所替代）等等等等，请结合手册自行尝试。



6 总结与展望

Multiwfn的推出大大拉低了波函数分析的门槛。以往许多十分重要的但只有专家们才会做的分析，如今变得可以被广大量子化学工作者所方便地利用，在实际问题的研究中已产生了巨大价值。也正因为拉低了门槛，使不少原先被束之高阁少有人问津的却又有重要实际价值的波函数分析方法进入了广大量子化学研究者们的视野。波函数分析的概念、意义想必也因此更加地深入了人心。Multiwfn的开发不会终止，而是会不断发展、完善，变得更强大、更高效、更好用，成为量子化学领域公认的世界第一的波函数分析程序。

Multiwfn后续版本中目前预计将要加入的新功能或进行的改进有以下这些，其中前三项是优先考虑的。
1 支持动态RDG分析方法（dx.doi.org/10.1021/ct4001087），也就是可以用RDG方法分析分子动力学轨迹，描述出动力学过程中的平均弱相互作用方式
2 加入直接耦合法和格点能修正法(SEC)计算电荷转移积分计算功能，并且能根据Marcus理论直接输出电荷转移速率常数。
3 AOSB方法，见J.Mol.Struct.(THEOCHEM),870,1。此方法可以将Mayer键级直接分解为sigma、pi、delta的贡献。
4 增加轨道定域化功能，优先考虑Pipek-Mezey定域化方法。
5 支持分布多极分析(DMA)。可以将体系电荷分布展开为大量点多极矩来描述。
6 支持更多的程序，如Crystal09、ADF等，还可能将支持基于平面波的程序，使Multiwfn可以在周期性体系的分析上大展拳脚。
7 提高等值面图的显示质量、提升静电势计算速度
如果读者从文献中发现了或者自行提出了一些很有用的波函数分析方法但Multiwfn尚不能实现，可以联系笔者(sobereva@sina.com)来将这些方法引入Multiwfn。

Multiwfn也有开发在线版本的可能（只打算实现最重要的核心功能）。目前设想是运行在公开的服务器上，用户直接画个分子结构（或者上传波函数文件），设定一些简单的选项，任务就会提交到后台进行，然后返回到用户的浏览器中。基于这样的在线版本，甚至连根本没用过量化程序的人，比如非理论化学专业研究者、学生们也能借助于Multiwfn轻松地分析体系的电子结构特征，使得Multiwfn的受众人群能扩展到所有化学研究者。另外，有了在线版本，用户甚至也可以在手机、平板电脑上来做波函数分析。

[ Last edited by sobereva on 2013-4-29 at 23:12 ] 返回小木虫查看更多