Multiwfn从1.0到3.0开发历程回顾.pdfVIP

思想家公社的门口一量化“分子模拟”二次元

此blog全部文章索引：

。计算化学公坛：

。思想家公社（讨论计算化学为主，加入

需注明研究方向）。Sobereva的硬件资料库：

[]杂谈Multiwfn从1.0到3.0版的开发经历

2013-03-2921:48

杂谈Multiwfn从1.0到3.0版的开发经历

文/SoberevaFirstrelease:2013-Mar-29

在Multiwfn3.0发布后，在疲惫、激动之余，回想起从Multiwfn1.0到3.0版这三

年多断断续续且艰苦的开发历程，颇有点感慨。当初只是一个很小的念头，卧在床

上写了几百行的程序，想不到这个小程序如今成为了Multiwfn3.0，有了约三万行

代码，被近百篇文章所使用，用户覆盖、欧洲各国、、甚至伊朗。这

段开发历程是个珍贵的，不禁想写个帖子零碎地闲谈一下，写写自己当时

的想法，等以后很可能就淡忘了。详细介绍某量子化学软件的开发，特别是记

录开发者自己每个阶段想法的文章我未曾见到，但如果有，我认为一定会很有趣很

值得一读。等再过几十年，有人尝试总结一次量化软件发展史的时候，说不定

本文的一些闲杂内容还能当做素材呢。

Multiwfn是波函数分析软件，所以，这篇闲谈自然要上溯到笔者什么时候开始对波

函数分析开始有。

笔者开始量化大抵是从大二的时候。后来听说到有AIM这么个理论，貌似很

有趣，但不明内容。大概在大三的时候，在网上偶然看到一个23页的Bader写的

综述AIM的pdf文档，是总共有五卷售价奇高的EncyclopediaofComputational

Chemistry中的一节。当时无聊的课程、实验多得很，我把那23页打印出来，挤出

时间在上间看、在做实验时等着反应完成的时间看。虽然当时尚没法完全读懂，

但是起码那些读懂的内容，比如键径，令我觉得应该用fascinating来形容，相当令

我着迷。当时还萌生了个念头，想把AIM和前线轨道理论给融合起来来更好地解

释反应如何进行。不过，我现在已经对前线轨道理论有点嗤之以鼻了，大抵是后来

接触了概念密度泛函的缘故吧。接触AIM，应该说是我对波函数分析感的最初

契机。我最早接触的AIM程序则是AIM2000。记得是在一个酷热的夏天，在宿舍

里热得难以入睡，将AIM2000demo版对照着帮助文件把玩了一晚上，基本都玩通

了。好不容易把电子密度的等值线+梯度线图做出来，看着挺欣喜的，还专门

写了个简单的AIM2000中文。Multiwfn的拓扑分析模块的操作方式的设计，

有很多方面都是借鉴了AIM2000。

GsGrid如今早已停止开发，其功能已经全部融入进了Multiwfn，或许有人还记得这

个程序。这个程序，和Multiwfn后来的开发有很微妙的联系。开发GsGrid的动机

是因为，在那个时候，网上有人问怎么做出诸如电子密度之类实空间函数的平面图。

虽说这种图在文献里到处都是，可是，起码对于高斯来说，当时要想做出这种图还

真是没有好法子。虽然gview已经可以直接做出截面的等值线图了，但这也是后来

的事了（当时还是以gview3.0.7版本为主）。虽然我那时候倒不需要做平面图，

但是为了给大家带来方便，用了好像是一个通宵就写出了最初的GsGrid程序（然

后当天发布到网上，中午我就给我亲戚修电脑去了）。它用来高斯生成的

cube文件，然后提取其中指定平面上的数据，以便用sigmaplot等程序作图。

GsGrid中的gs代表Gaussian。（后来在学术搜索中发现已经有其它程序叫

这个名字了，早知道如此可能当初就换个名字了）

后来逐渐给GsGrid加入了更丰富的功能，版本号也一直在提升。其中对格点文件

相互运算的功能很有用，比高斯自带的cubman好很多。小木虫上zhou2009对

GsGrid的开发给予了热情的鼓励，他的的帖子很多都用了GsGrid，特别是用来做

密度差平面图。最初GsGrid只是在上发布，后来建了个主页，不过的免

费网络服务历来都很不靠谱，说停就停，一点责任心都没有，记得是换了两三次主

页空间，最后，在微软