注：此文的内容以后有可能更新，在blog中可看到最新版本http://hi.baidu.com/sobereva/item/932f96260c57eaf951fd8716

Multiwfn支持的弱相互作用的分析方法概览

文/Sobereva(3)  2014-Aug-28

这里所谓的弱相互作用，是指强度明显弱于一般化学键的各种形式的相互作用，范德华作用、pi-pi堆积、氢键、卤键、二氢键等等皆属于这一范畴。分析弱相互作用的方法非常多，Multiwfn程序可以方便地实现其中大部分分析方法，在此文将对它们进行提纲掣领的归纳，但并不进行具体讲解，而是给出相关参考资料，多数是已发表的帖子和手册里的章节，看了之后很快就能明白原理和操作。此文内容对应的是高于3.3.5版的Multiwfn的手册和程序，老版本可能不支持其中一些功能或手册缺少相应章节。

在http://multiwfn.codeplex.com可以免费下载到Multiwfn所有版本的程序、源代码和手册。如果对Multiwfn不了解，建议阅读《Multiwfn波函数分析程序的意义、功能与用途》（http://hi.baidu.com/sobereva/item/b35581418db96cea1f19bc09）和《Multiwfn入门tips》（http://hi.baidu.com/sobereva/item/896ee3a19f6d7d3c020a4d76）。

下面开始依次列举Multiwfn支持的弱相互作用分析方法，顺序是随意的。


1 Bader的Atoms in Molecules (AIM)拓扑分析。这是最常用、最经典的弱相互作用分析方法，通过相互作用的原子间的键临界点(BCP)的性质来分析弱相互作用特征。理论介绍可参考AIM入门资料和重要文献小合集（http://pan.baidu.com/share/link? ... p;amp;uk=1074012119）、Multiwfn 2013年培训班的day1里的《Some concepts in wavefunction analysis》（此次培训班所有幻灯片可在这里下载：http://pan.baidu.com/share/link? ... p;amp;uk=1074012119）。Multiwfn的AIM拓扑分析是诸AIM程序中效率最高、使用最方便、最灵活的，使用方法参见《使用Multiwfn做拓扑分析以及计算孤对电子角度》（http://hi.baidu.com/sobereva/item/c84fa5089156806dd55a11c3）或手册4.2.1节的例子。Multiwfn可以直接图形观看拓扑分析结果，但也可以结合VMD来绘图以对图像效果进行更灵活自由的设定，见《Multiwfn结合VMD绘制AIM拓扑分析图》（http://hi.baidu.com/sobereva/item/40c03f072667a931a3332a8e）


2 杨伟涛等人在JACS,132,6498(2010)提出约化密度梯度(RDG)分析是一种极为重要、优秀的弱相互作用分析方法，提出后立刻得到十分广泛的使用。它在某种意义上是对AIM理论的扩展，使分析结果更为清晰直观，而且还能处理许多原本AIM理论不适用的场合，如pi-pi堆积、无BCP出现的情况（如JCTC,9,3263）。RDG方法的理论介绍和使用见《使用Multiwfn图形化研究弱相互作用》（http://hi.baidu.com/sobereva/item/bb1162d881b8d5ee3cc2cbdc）或手册4.100.1节，在Multiwfn 2013年培训班的day3里的《Analyze weak interaction by RDG method》当中也有大量图形示例。RDG分析原本是对特定结构进行分析的，随后还被扩展到动态环境的分析，比如可以分析分子动力学模拟过程中的弱相互作用，见《用Multiwfn研究分子动力学中的弱相互作用》（http://hi.baidu.com/sobereva/item/bef6c53314d400be124b147f）或手册4.200.1节。RDG还能和ELF联用同时展现强相互作用和弱相互作用，见JCTC,8,3993的讨论和《通过键级曲线和ELF/LOL/RDG等值面动画研究化学反应过程》（http://hi.baidu.com/sobereva/item/e839b2d2efd0aadf241f4015）。


3 静电势分析。一般的氢键、卤键、二氢键都是以静电相互作用主导的，都可以归纳为Lewis碱（富电子区域）和Lewis酸（缺电子区域）之间的静电吸引作用。静电势是与静电效应关系最为密切的实空间函数，十分适合分析各种静电主导的弱相互作用。如果对静电势不了解，建议看看静电势综述文章小合集（http://pan.baidu.com/s/1stw8L）和《Multiwfn与波函数分析简介》（http://hi.baidu.com/sobereva/item/4fc0c9d3efd0aadf251f408d）中的相关部分。静电势具体实现方法有很多：
3.1 直接考察静电势填色的分子范德华表面，这可以在定性层面上直观了解体系哪个部位容易发生什么类型和强度的静电相互作用。很多文章都是这么分析的，比如JMM,13,291(2007)。用gview作这种图算是最方便的，而Multiwfn结合VMD作这种图则可以同时在图上展现更丰富的定量层面的信息，即静电势极大点、极小点位置和数值（见下节）。另外Multiwfn还可以给出不同静电势数值区间内的分子表面积，这使得静电势分布方式能被了解得更为深入透彻。操作见《使用Multiwfn结合VMD分析和绘制分子表面静电势分布》（http://hi.baidu.com/sobereva/item/692f9a83d8e9a5c8b071545d），应用实例可参看笔者的DOI：10.1007/s11224-014-0430-6、J.Phys.Org.Chem.,26,473(2013)。

3.2 研究分子范德华表面上静电势极大点、极小点位置和数值。其位置可以用于研究发生相互作用的位点，而数值则和静电相互作用强度密切相关。这种分析方式被广为采用，特别是Politzer和Murray的大量文章都是这么分析的。相关应用实例非常多，诸如PCCP,15,14377(2013)、JMM,13,305(2007)、IJQC,107,3046(2007)、PCCP,12,7748(2010)、JMM,14,659(2008)、JMM,18,541(2012)、JMM,15,723(2009)、Practical Aspects of Computational Chemistry (2009)一书的第六章。获得这些分子表面极值点可以通过Multiwfn的定量分子表面分析模块实现，操作参见《使用Multiwfn的定量分子表面分析功能预测反应位点、分析分子间相互作用》（http://hi.baidu.com/sobereva/item/6601fff28c4fd9d643c36acc）或手册4.12.1节，更详细的操作和算法介绍可参见手册3.15节和笔者的J.Mol.Graph.Model.,38,314(2012)。

3.3 静电势等值线重叠分析。笔者在J.Mol.Model.,19,5387(2013)中提出的这种分析静电相互作用的方法十分直观，将单体的静电势等值线图简单地相互叠加就能一目了然地估出二聚体的不同构型的相对稳定性。在《静电效应主导了氢气、氮气二聚体的构型》（http://hi.baidu.com/sobereva/item/fe74588e6c0223ebe496e088）中对这种方法做了介绍。在Multiwfn中利用功能4可以十分便利地绘制高质量的静电势等值线图，见手册4.4.4节的例子。

3.4 在JPCA,118,1697(2014)中，Mohan和Suresh发现静电势在相互作用原子的原子核位置的数值在单体和复合物环境下的变化量与相互作用能有着奇好的线性关系，这个关系非常有用。在手册4.1.2节对此方法进行了介绍，并举例如何计算。


4 原子电荷分析。原子电荷是位于原子核的点电荷模型，是对化学体系电荷分布最简单、直观的描述。通过原子电荷的符号和大小，可以直观了解哪些原子或哪些片段间有着明显的静电相互作用。计算原子电荷的方法非常多，近100种，笔者在物理化学学报,28,1(2012)中对主流的原子电荷计算方法进行了详细的对比讨论。只有对静电势重现性好的原子电荷，如笔者提出的ADCH电荷（J.Theor.Comp.Chem,11,163(2012)）以及CHELPG、MK、RESP、AM1-BCC等基于拟合静电势得到的电荷才能较好用于分析静电相互作用；而NPA电荷（总被人谬称成为NBO电荷）、Mulliken电荷，尤其是AIM电荷等都不适合此目的。注意原子电荷分析是有局限性的，像是导致卤键出现的本质原因，即σ-hole，是由于卤原子的电子的各向异性分布造成的，而原子电荷相当于将原子附近的电荷球对称化了，故不可能用来讨论卤键。主流的量化程序如Gaussian本身可以算一些类型的原子电荷，在Multiwfn中也能算许多种类的原子电荷，见手册3.9节的介绍和4.7节的一些例子，其中有些是Multiwfn独家支持的，如ADCH电荷。


5 键级和离域化指数分析。各种键级，以及AIM理论里的离域化指数(DI)主要涉及的是共价作用，而弱相互作用中共价成份通常可以忽视，因此键级和DI分析对这类问题通常用处不大。但是，弱相互作用范畴中也有强度较强的，比如低势垒氢键(LBHB)、电荷辅助氢/卤键等，共价作用就有一定贡献了，就可以适当讨论键级和DI了。键级的定义很多，前述的《Multiwfn与波函数分析简介》里进行了简要介绍，这里说的键级具体指Mayer键级和模糊键级。实际上，Mayer键级、模糊键级、DI的物理本质都一样，都是衡量原子间共享的电子对儿数，只不过它们分别是在希尔伯特空间、模糊原子空间和AIM原子盆下计算的，所以实际数值有一些差别，特别是对于极性键而言。各种键级都可以在Multiwfn里计算，见手册4.9节的例子和3.11节的介绍。DI可以在Multiwfn的盆分析功能中计算，参见《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》（http://hi.baidu.com/sobereva/item/47427a1487ce95633e87ce25）和手册4.17.1节的例子。通常还是推荐用Mayer键级来讨论，计算几乎不需要耗时，但不适合基组有弥散函数的情况，此时应当用模糊键级或DI。
（顺带一提，在一些文章中，比如DOI: 10.1039/c3cp51077h，作者用轨道、键级来讨论本质是色散作用的pi-pi堆积，发现达到pi-pi堆积最佳距离时键级达到最高，这绝对不代表pi-pi堆积的本质是共价作用，不要被文中的一些说法所迷惑！实际上，色散作用才是pi-pi堆积的“因”，而轨道和键级的变化只是由于分子间被拉近后出现的“果”）


6 电子定域化函数(ELF)分析。如果对ELF不熟悉，建议看看前述的《Some concepts in wavefunction analysis》中的相应介绍和《定域性的图形分析》（http://hi.baidu.com/sobereva/item/4a8aeae2291d8d0f8c3ea8d8）。一些ELF综述和相关重要文献见http://pan.baidu.com/share/link? ... p;amp;uk=1074012119。更深入的关于定域性本质的讨论见《电子的定域性与相关穴》（http://hi.baidu.com/sobereva/item/0fdcf532597122352f20c4c3）。ELF通常是用来分析强相互作用范畴的问题，但是对于偏强的弱相互作用分析，ELF也是有用的。在Theor. Chem. Acc.,104,13(2000)中Fuster和Silvi基于ELF定义了CVI指数来对氢键强度进行分类，之后在JPCA,115,10078(2011)中作者又利用CVI等方法研究了一大共振辅助氢键，发现CVI指数和其它与氢键强度相关的指标有较好的相关性。只要弄懂了CVI的定义，利用Multiwfn的主功能3在氢键涉及的原子间绘制ELF曲线图，并令Multiwfn找出曲线极小点位置和数值，马上就能得到CVI指数。还有其它一些文章也利用ELF分析氢键，比如Chem.Rev.,111,2597(2011)，它们有的涉及到ELF盆分析，可以参考手册4.17.2节的例子；而有的涉及到绘制ELF等值面图，在手册4.5.1节有例子。另外如果结合Chimeria还可以把ELF图画得更好看，把不同特征区域以不同颜色区分，见《在Chimeria中用不同的颜色显示不同的等值面》（http://hi.baidu.com/sobereva/item/c5164fd0efd0aadf241f40f4）。


7 电荷变化分析。弱相互作用常伴随着电荷的变化（尽管程度往往并不大），研究电荷是如何在分子间转移的，在出现另一个分子后体系的电荷是如何被极化的，这都是值得关注的问题。Multiwfn提供了大量用于研究电荷分布变化的功能：
7.1 绘制复合物与单体间的电子密度差图。这是最重要、直观、常用的研究电荷变化的手段，在《使用Multiwfn作电子密度差图》（http://hi.baidu.com/sobereva/item/96eaeb5dd1b12309e7c4a5c1）和手册4.5.5节有相应的例子。从密度差图上经常会看到复合物形成后某些原子间密度下降了，而有些升高了，表现出相应键的削弱和增强。

7.2 绘制电荷转移曲线。在利用Multiwfn的主功能5生成了复合物与单体间的密度差格点数据后，对它进行指定方向的积分，就能产生所谓的电荷转移曲线。这对于定量分析电子在某一方向上是如何转移的极为有用。例如一个分子吸附在一个界面上，就可以在垂直于界面方向进行积分，在什么位置电子净转移了多少一目了然。详见3.16.14节的介绍和4.13.6节的例子。

7.3 对单体和复合物状态下的原子电荷求差值得到变化量，可以再加和成片段或分子的变化量。这样一来，弱相互作用发生前后每个原子、片段或分子的电子增减了多少就非常清楚了。这个方法极为重要、常用！

7.4 生成复合物与单体间的密度差格点数据后，可以再进一步利用Multiwfn的盆分析功能对密度差的盆进行积分，这样就可以定量获知不同局部特征区域（如σ-hole区域、孤对电子区域）的电子增减情况。做法参见前文提到的《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》或手册4.17.4节。

7.5 电荷分解分析(Charge decomposition analysis, CDA)。通过7.3的做法，我们只能知道电荷净变化是多少，但不知细节，而CDA可以展现内在细节，能获知片段之间是如何通过复合物的各个分子轨道来贡献和反馈电子的，另外Multiwfn的CDA功能还可以给出复合物轨道是如何由单体轨道混合而成的，并能直接输出轨道相互作用图。通常CDA被用于强相互作用范畴的研究，比如金属和配体间的电荷转移，但是如果将弱相互作用体系中的单体定义为片段，CDA对于讨论弱相互作用导致的电荷转移也是比较有用的。CDA的原理介绍和实例参见《使用Multiwfn做电荷分解分析(CDA)、绘制轨道相互作用图》（http://hi.baidu.com/sobereva/item/178fbd02969cb690a2df4375）。

7.6 Multiwfn有个专门基于密度差数据来分析电荷转移激发的功能，能得到很多有用的统计量，比如电子转移距离、转移电子分布的弥散程度等，并且还能进一步将正负值区域繁复交错的密度差数据变换为十分直观简单的描述形式。见手册3.21.3的介绍和4.18.3节的例子。这个方法虽然没有被人们用于单体间电荷转移的研究，但我想没准儿也会有用，特别是如果从密度差图上看电荷有明显的整体转移倾向的话此方法很值得一试，只要把手册例子中计算激发态与基态间的密度差那一步改成计算复合物与片段间的密度差即可。


8 Hirshfeld表面分析。这个方法对于研究分子晶体内的弱相互作用极为有用！综述参见CrystEngComm,11,19 (2009)。此方法可以十分清晰地给出分子间相互作用的界面，并从界面不同区域的颜色上能立刻了解到什么位置发生了什么类型的作用。并且还能绘制指纹(fingerprint)图来定量研究。见手册3.15.5节的理论介绍和4.12.5、4.12.6节的实例。


9 LOLIPOP是Chem.Commun.,48,9239(2012)中提出的一种基于定域化轨道指示函数(Localized orbital locator, LOL)定义的预测分子间pi-pi堆积能力的指数，和相互作用能有较好关系。介绍见手册3.100.14节，实例见4.100.14节。


10 能量分解分析。能量分解分析是研究弱相互作用极为重要的手段，可以将总相互作用能分解成各种有确切物理意义的成份，比如静电、交换互斥、色散、轨道相互作用能各是多少。能量分解分析方法很多，而且分解出的能量项定义也各不相同。最常用的能量分解方法有两类，一类是Morokuma的EDA能量分解，包括后来的进一步发展，如GAMESS-US所支持的苏培峰提出的LMO-EDA和GKS-EDA。另一类是对称匹配微扰理论(SAPT)，以及后来的一些发展，特别是现在流行的DFT-SAPT。SAPT的十分简明扼要的介绍可以参看JMM,19,5387的补充材料，综述见DOI:10.1002/wcms.86。另外还有NBO里的自然能量分解(NEDA)等方法，由于不流行就不多说了。EDA形式的能量分解可以用GAMESS-US、ADF、Q-Chem做，SAPT可以用PSI4、Molpro、Q-Chem或Szalewicz等人开发的独立的SAPT代码来做。不过，ADF、Q-Chem、Molpro、Q-Chem都是收费的，而其它的对于缺乏一定计算机水平的量化初学者可能门槛稍高。利用Multiwfn的片段组合波函数功能，能方便地做一种很简单的能量分解，结合Gaussian和Grimme的DFT-D3程序使用可以把总相互作用能分解为色散、静电+交换互斥、轨道作用这三项，见手册3.100.8节的说明和手册4.100.8节的实例。DFT-D3程序也很容易用，笔者专门写过个帖子《DFT-D色散校正的使用》（http://hi.baidu.com/sobereva/item/772d38c6e13165d497445243）。之所以说Multiwfn做的是“很简单的能量分解”，是因为这种做法没法拆开静电和交换互斥。


11 源函数(Source function)估计知道的人不多，目前用的人也较少，这是Bader和Gatti在AIM框架下定义的一个实空间函数，可讨论不同空间区域对指定的一个点的电子密度的贡献（即讨论密度的源头）。可以参看Gatti写的一个极尽全面的源函数综述Struct.&Bond.,147,193(2010)，特别涉及了氢键的分析。源函数的简明介绍见手册2.6节相应部分，在Multiwfn中的分析实例见4.17.5节。


12 范德华表面相互穿透距离（Mutual penetration distance）分析。弱相互作用会导致分子的范德华表面被彼此穿透。对于同种类型弱相互作用，穿透的距离越大，显然弱相互作用越强（但如果穿透是由于其它外界作用而被迫产生的则不满足这个规律）。假设相互作用的两个原子分别是A和B，最简单的计算穿透距离的方法就是把A的范德华半径加上B的范德华半径减去A-B的距离。但这么做比较粗鄙，范德华半径定义本来就存在任意性，而且简单地用圆球模型也无法很好体现真实的分子范德华表面的形状。最精准、直观的计算范德华表面相互穿透距离的方法参看前文提到的《使用Multiwfn结合VMD分析和绘制分子表面静电势分布》的最后一节，把相互作用的取向都考虑进去了。但是此方法可能操作上稍麻烦。一个简单点的但也比较准确的方法是在Multiwfn中做完定量分子表面分析后选10得到原子的非键原子半径（non-bonded atomic radius），也就是其原子核距离范德华表面最近的距离。将相互作用的两个原子的非键半径相加减去它们的实际距离即是范德华表面相互穿透距离。


13 原子多极矩分析。原子多极矩在手册3.18.3节和前述的《Some concepts in wavefunction analysis》相应部分有介绍。前面提到，原子电荷（相当于原子单极矩）模型把原子附近电子分布给球对称化考虑了，而原子多极矩则正是用来描述原子附近电子分布的各向异性的。原子间静电相互作用可以用多极展开来描述，首项是单极-单极作用，而随后的项就与偶极矩、四极矩相关了（可参看Molecular modelling Principles and applications 2ed的4.9.1节的介绍），因此要更细致地讨论分子间静电相互作用就需要把原子多极矩也进行考虑。一些分析实例参见Bader的著作Atoms in molecules-A quantum theory中的7.4.3节。在Multiwfn中既可以基于模糊原子空间计算原子多极矩，也可以基于AIM原子盆来计算（见手册4.17.1节的例子），最高能算到八极矩。


14 轨道重叠分析。虽然轨道相互作用前面也提了，通常不是弱相互作用的重点，更不是弱相互作用的本质驱动力，但依然可以去讨论，毕竟分子离近了之后势必轨道重叠得更厉害，并因此伴随着轨道混合的出现。轨道重叠越厉害，通常轨道作用往往越强，在《分子间轨道重叠的图形显示和计算》（http://hi.baidu.com/sobereva/item/7cc644e165d9b10c65db00ae）介绍了怎么用Multiwfn计算两个单体轨道间的重叠积分，并且还利用Multiwfn的灵活性实现了一种新颖实用的可视化轨道重叠的方法。NBO框架中的二阶微扰能E(2)分析也和相应两个NBO轨道间的重叠程度有密切关联，较大的重叠是E(2)有较大负值的前提，用Multiwfn也可以很方便地观看NBO轨道以及NBO轨道间的重叠，见《使用Multiwfn绘制NBO及相关轨道》（http://hi.baidu.com/sobereva/item/b0226b3a534a4b637d034bc5）。


15 电子密度拉普拉斯函数通常用于研究强相互作用范畴的问题，它是电子密度的二阶导数，越负（越正）说明电子在局部越凝聚（发散），前述的《Some concepts in wavefunction analysis》和《定域性的图形分析》对之有较多介绍。这个函数有时也被用来讨论静电主导的弱相互作用，因为它的符号和大小能体现原子周围哪里富电子（显局部负电）和缺电子（对核电荷屏蔽不足，故因此表现正电性）。在JMM,19,2035(2013)中，作者在研究卤键的时候把σ-hole位置的电子密度拉普拉斯函数的(3,-1)型临界点，即最缺电子的那个位置的密度与相互作用能很好地关联了起来。JPC,95,4698(1991)讨论的是SN2反应，思路也有些许相似，是基于被进攻的原子的电子密度拉普拉斯函数来讨论局部缺电子的程度，由此考察亲核试剂进攻的方向和难易程度。在前文提到的《使用Multiwfn做拓扑分析以及计算孤对电子角度》中演示了怎么用Multiwfn强大的拓扑分析功能搜索ELF临界点，搜索电子密度拉普拉斯函数的临界点与其过程完全一样，只不过被研究的函数选择成电子密度拉普拉斯函数即可。




弱相互作用分析方法甚多，但是只要把上述分析方法都掌握了，灵活运用，就已经能研究分析得相当深入了。还有一些分析也经常被涉及，但与Multiwfn没直接关系，下面乱谈一下。

NBO E(2)分析：这个方法在弱相互作用分析中被用得很普遍，但实际上总是被滥用！有些人用E(2)讨论相互作用能，甚至以为分子间的E(2)能量加和就能等于相互作用能，这真是大错特错！很多人没有弄懂其本质，就被大量用E(2)的文献给误导了。E(2)分析的本质是轨道相互作用，前面已经提过几次，轨道相互作用在一般强度的弱相互作用中只占极小的角色，更不是复合物中使单体能够保持结合的本质原因。色散和静电相互作用才是弱相互作用的重点，没了色散或静电作用复合物就无法维持。像是一般的氢键，仅凭考虑电多极矩的力场（如杨忠志的ABEEM），尚且没有考虑极化，更没考虑轨道作用，就已经能得到准确的氢键能量，这正表明氢键的本质是静电作用。传统的NBO分析框架里对分子间相互作用的描述根本没有色散也没有静电作用。某人的一句话说得好：“NBO里什么都是超共轭！”，本身不是超共轭的在NBO里也非得用超共轭的形式来描述，超共轭能量也就是常说的E(2)值。拿讨论轨道作用的工具E(2)讨论静电和色散主导的作用，完全就是喧宾夺主，根本没有正确反映相互作用的物理本质。但不是说不能用E(2)讨论，关键是一定要搞清楚因果关系，理解本质，否则必然会对结果给出错误的解释。前面也提到，静电、色散作用使得单体相互靠近、结合，从而造成了轨道重叠的加剧、发生混合，两个单体的电子间也更容易相互离域（共享），随之而来的就是那些衡量共价性的指标，如键级、E(2)数值的少量增加，这是为什么这些衡量共价作用的指标也能或多或少和静电/色散主导的相互作用的能量找出相关性。但如果理解成正是因为有了超共轭，才形成了复合物，作用能就是E(2)值，那就彻底错了。况且，E(2)本来只是一个粗糙的值，计算方式很简单，即便是真正意义的超共轭作用，E(2)对它的作用能的衡量也不是精确的。另外一些人还把涉及相同NBO的多个E(2)数值相加讨论弱相互作用，这样的做法更是严重缺乏物理依据。总之，别再滥用E(2)分析弱相互作用了！

自然布居分析(NPA)：NPA是NBO分析中的核心组成部分。复合物形成前后原子电荷的变化只能反映整个原子的电子得失，而利用NPA分析则可以得知每个原子轨道的电子得失，以及原子杂化状态的改变，这是研究所谓的重杂化效应的关键。

ETS-NOCV：这是Ziegler等人在JCTC,5,962提出的把能量分解和密度差搅合在一起的分析方法，特点是能把片段间的轨道相互作用进行进一步划分，每一部分有它对轨道相互作用能的贡献并能给出与之对应的密度变化图形。ADF支持此方法。

磁感应电流：如果对此概念不熟悉，建议看看《使用AICD程序研究电子离域性和磁感应电流密度》（http://hi.baidu.com/sobereva/item/07d17c23d959de8e9d63d1c6）、《使用GIMIC计算和分析磁感应电流密度》（http://hi.baidu.com/sobereva/item/2fcd6479efa95d346cc37ccc）里的一些基本介绍。在GIMIC综述文章PCCP,13,20500(2011)中提到磁感应电流穿过氢键的强度可以用来衡量氢键强度。其原因就在于氢键越强，相互作用的原子间电子的离域能力越强，感应电流就越容易穿过去。其实这种分析方法，虽然有物理意义，但实际意义不大，计算麻烦，何况都已经有那么多比这好用得多的衡量氢键的方法了。

键长和振动频率变化分析：比如X-H---Y这样的氢键，氢键形成后通常伴随着X-H的键的削弱，因此键长会增长，振动频率也会降低（红移），但也有时相反。键长和频率的变化是氢键文章中经常讨论的。建议同时结合键级来分析。卤键、二氢键等也可以类似地分析。

Mayer能量分解：此方法的原理和使用见《Mayer能量分解及APEX4程序使用简介》（http://hi.baidu.com/sobereva/item/2358dc8bd8e9a5c8b07154de）。此方法和前面提到的能量分解方法不同之处在于它不需要定义片段，而是直接对整体进行分析，得到体系中不同原子间的各种成份的相互作用能（不含色散作用）。Mayer在这方面虽然做出了不少努力，不断发文章改进、探讨，但始终不流行，此方法在弱相互作用上的应用也极少。

Interacting quantum atoms (IQA)：此方法用处和Mayer能量分解类似，也是能给出体系中每对儿原子间的相互作用能的物理成份，但流行度还说得过去。此方法最初在JCTC,1,1096(2005)中提出，强、弱作用都能分析，一篇分析卤键的例子是JPCA,117,8969(2013)。这种方法在原理上比较简单，但算起来却非常困难。因为此方法是用AIM原子盆定义原子空间，本来积分AIM盆就是众所周知的极费时间，而此方法则需要对两个AIM盆之间做二重积分，耗时可想而知。所以此方法目前只能用于非常小的体系，实用性不高，目前公开的程序中仅AIMALL支持。Multiwfn以后可能也会支持，但不是基于AIM盆，而是基于Hirshfeld方法定义原子空间，这样计算量就小多了，结果也不会比IQA差。



方法介绍就说到这里。其实研究弱相互作用还是挺套路化的。寻找或构造一类新的弱相互作用体系，优化一下，算算振动频率，用高精度方法算一下相互作用能（适宜的方法可参见《乱谈DFT-D》http://hi.baidu.com/sobereva/item/4e18b2da4f023317e1f46fde），做个能量分解，然后AIM啊、RDG啊、原子电荷啊、密度差啊、静电势分析啊、范德华表面穿透距离分析啊等等前面提到的方法往上一通招呼，把结果一总结、对比，将得到的一些量与相互作用能拟合个关系，再和已知的其它体系的作用相互对照讨论，一篇文章就出来了，不难，但切勿照搬套路疯狂灌水。另外也不要乱用前面提到的那些方法，而是应该在理解分析方法的物理意义的基础上选择最适合自己体系的方法。

分析弱相互作用的方法肯定还有很多这里没有提到，对于那些涉及到电子结构的分析方法，只要灵活运用Multiwfn的现有功能，大多都能轻松实现。 返回小木虫查看更多