简介：Multiwfn波函数分析程序资源包，面向量子化学与计算化学领域的学生和科研人员，提供波函数分析、电荷分布、自然键轨道、多体相互作用能量等功能，并配有大量可直接运行的示例文件，可帮助使用者快速上手这款主流后处理工具。压缩包共883个文件，大小约18.95MB，涵盖rad、gjf、wfn、fch等量子化学输入输出格式，以及txt说明、out日志、vmd可视化脚本和sh批处理脚本，便于对照不同计算场景练习。已有729人学习下载。资源内示例分子丰富，包括丁二烯、甲酰胺、菲、乙醇等体系，涉及HOMO/LUMO、Mulliken电荷、振动光谱、氢键与电荷转移分析等典型任务。通过逐项拆解这些文件，使用者既能理解Multiwfn的操作流程，也能掌握从基础性质计算到反应机理探索的完整方法，适合系统自学或教学辅助。 做量化计算这几年，我电脑里装得最勤快的程序，不是Gaussian，不是ORCA，而是Multiwfn。每次算完一个体系，拿到能量和轨道之后，真正的重头戏才开始——分析波函数。Multiwfn就是把这件事做到极致的工具。它免费、跨平台、功能密度极高，从基础的Mulliken布居分析到复杂的ETS-NOCV成键分析、从头算分子动力学轨迹分析，几乎覆盖了你对波函数分析的所有想象。这篇内容不打算做成官方文档的复读，我想从实际使用的角度，把Multiwfn的选型逻辑、核心功能、实操套路和踩坑记录都摊开来讲，希望能帮你少走点弯路，尤其是刚接触量化计算的朋友。
1. 先说清楚：Multiwfn到底是干嘛的
1.1 从“算完能量之后”说起
做量子化学计算，最典型的流程是：选一个方法（比如B3LYP/6-311G**）、建模、提交任务、等计算结果。很多人拿到的结果就是一个能量值、一组轨道能级，或者一张HOMO/LUMO的图片。但实际研究中，能量只是冰山一角。电荷怎么分布？哪里亲电、哪里亲核？两个分子之间的作用力本质是什么？芳香性到底强不强？这些问题，都需要对波函数本身做深入分析。
Multiwfn就是一个专门干这个的通用波函数分析程序。它读入量化程序输出的波函数文件（比如Gaussian的.fchk、.wfn、.molden，ORCA的.molden，CP2K的.molden等），然后在这个基础上执行各种分析任务。它不是一个量子化学计算引擎，它不负责求解薛定谔方程，它只负责在你已经算完的基础上，把波函数里藏着的化学信息一点点挖出来。
我刚接触它的时候，感觉这就是个“命令行黑箱”，界面全是字符菜单，没有图形界面，第一印象确实劝退。但用顺手之后会发现，这种设计反而带来极高的效率：所有操作都可以脚本化，批处理非常方便，而且程序本身极其轻量，单文件可执行程序，放在U盘里都能跑。
1.2 为什么不用Gaussian自带的分析功能？
这是个很实际的问题。Gaussian、ORCA自带的分析功能并不少，比如Mulliken电荷、NBO分析、轨道成分等等。但自带的工具通常有个问题：每类分析都是独立模块，数据格式不统一，跨程序对比需要花大量时间做格式转换，而且很多高级分析根本不在自带功能里。
举一个我亲历的例子。之前我需要做IGMH（一种弱相互作用可视化方法）分析，同时对比三个不同泛函的结果。用Gaussian自带工具做，得先分别重新计算、然后手动提取数据、再用第三方软件画图，流程很割裂。而用Multiwfn，一次性载入三个.fchk文件，用批处理脚本统一跑IGMH分析，输出结果格式完全一致，画出来的图放在一起对比，立刻就能看出趋势。这个体验上的差距，用一次就很难回去了。
再看一个维度——功能覆盖面。Multiwfn官方的功能列表有几十项，从最基础的布居分析、Mayer键级，到学术界前沿的IQA能量分解、ADCH电荷、完全态空间平均方法，基本覆盖了当前理论化学期刊里可能用到的绝大部分分析方法。这意味着你只需要学一个软件的交互逻辑，就能完成多种分析任务，学习成本省下来的时间，足够你多算两三个体系了。
提示：Multiwfn的原理和功能细节，在卢天老师发表的论文（“Multiwfn: A multifunctional wavefunction analyzer”）里有非常详尽的描述，做研究需要引用时一定要引用这篇论文。
2. 核心功能地图：你能用它做什么
2.1 成键分析三件套：键级、布居、轨道成分
按我自己的使用频率排序，最常用的是成键分析相关功能。新手入门的时候，最想回答的问题通常是“这个键是单键还是双键？”“两个原子之间有没有成键相互作用？”这时候键级和布居分析就是最直接的答案。
Multiwfn里可以做Mulliken布居、Lowdin布居、自然布居（NPA）、Hirshfeld布居，以及基于这些布居导出的Mayer键级、Wiberg键级等。我一般默认用Hirshfeld布居，因为它的数值对基组的依赖性相对较小，而Mulliken在加弥散函数时经常出现不合理的负占据数（这种问题说起来都是泪）。
轨道成分分析也是高频使用功能。你只需要把HOMO、LUMO的轨道序号告诉程序，再用Hirshfeld方法或者SCPA方法计算，就能得到每个原子/片段在这些轨道里的占比。这个方法在分析配位化合物的时候极其好用，比如判断金属d轨道和配体π轨道的相互作用程度，输出结果可以定量说事，比对着轨道图“目测”要靠谱得多。
2.2 静电势分析与弱相互作用分析
如果说成键分析是“里子”，那静电势分析就是“面子”，因为它直接关系到分子识别、晶体堆积、药物分子与靶点结合这类实际问题。静电势（ESP）描述的是分子周围一个单位正电荷感受到的静电作用，正值区域代表缺电子、容易被亲核攻击，负值区域代表富电子、容易被亲电攻击。
Multiwfn可以非常方便地计算并输出分子范德华表面上的静电势极值点，还能导出适合用VMD、PyMOL等可视化的ESP映射图。做药物设计相关性研究的朋友应该知道，配体分子表面静电势分布与结合亲和力之间常常存在结构性关联，用Multiwfn算一下ESP，就能从定量角度验证这些关联。
弱相互作用分析方面，我之前提过的IGMH方法就是Multiwfn率先实现并推广的功能，比传统的RDG（约化密度梯度）方法给出更清晰的原子对相互作用图景。这个方法分析π-π堆积、氢键、卤键都非常直观，输出结果用VMD一画，相互作用的区域一目了然，审稿人看了也容易理解。
2.3 激发态与芳香性分析：进阶方向
再往上走，Multiwfn还能做电子激发分析。比如计算TDM（跃迁密度矩阵）、空穴-电子分析、电子-空穴重叠积分等，这些是研究激发态性质的重要工具。以前做这些分析要么自己写代码，要么用不太方便的商业软件，Multiwfn把所有步骤都收纳进来了，配合Gaussian的TD-DFT计算输出就能直接分析。
芳香性分析方面，Multiwfn实现了ACID（各向异性感应电流密度）、NICS（核独立化学位移）等计算所需的波函数分析部分，还可以做电子定域化函数（ELF）、LOL等实空间函数分析，判断芳香环的电子定域化特征。如果做的是平面共轭体系、杂环化合物这类研究，这些功能基本绕不开。
2.4 拓扑分析：QTAIM的省心实现
Bader的“分子中的原子”（QTAIM）理论是实空间成键分析的经典方法。用Multiwfn做QTAIM拓扑分析非常成熟，临界点的搜索和分类、键径绘制、原子积分都能完成。我印象最深的是它能自动批量搜索分子中所有临界点，不需要手动给初始猜测，这一点比当年的AIMAll要省心很多。
3. 实操案例：从Gaussian输出到关键分析
3.1 准备工作：用合适的方法算、存对文件格式
Multiwfn本身不做量化计算，所以第一步必须用量化程序生成波函数文件。以Gaussian为例，一个典型的高斯计算输入文件长这样：
%chk=water.chk #p B3LYP/6-311G** opt freq water 0 1 O 0.0 0.0 0.0 H 0.0 0.0 0.96 H 0.0 0.93 -0.12算完以后，用formchk water.chk water.fchk把检查点文件转成文本格式的fchk文件，Multiwfn就能直接读取。注意，先用freq做振动分析得到稳定点确认是没问题的，但如果你只需要波函数分析，可以不加freq，省时间。
如果你是ORCA用户，一般会在ORCA输入文件里加上! def2-TZVP这类方法关键词，同时确保输出molden文件。我常用的组合是：orca_2mkl会把.mko文件转成.molden，然后直接在Multiwfn里载入。
注意：
fchk文件保存的是波函数的所有信息，包括基函数、轨道系数、格点信息等，所以文件会比较大，这是正常的。不要因为大就删掉，后面很多分析都需要它。
3.2 静电势极值点分析实操
假设我现在已经拿到了water.fchk，想看水分子的静电势分布。启动Multiwfn，输入文件名：
./Multiwfn water.fchk进入主菜单之后，选择静电势相关选项。在Multiwfn的主菜单中是选项3（“Electrostatic potential analysis”），然后选“Evaluate ESP on vdW surface and find surface minima/maxima”。它会自动计算范德华表面的静电势，然后统计出极值点。
输出里会给出极值点的坐标和ESP值，以及对应的原子。像水分子，氧原子附近会出现一个负的极小值，这就是孤对电子区域的富电子特征；氢原子那侧会出现正的极大值，对应缺电子的酸性氢。这些数据放在论文里，可以直接支撑你关于亲核/亲电位点的讨论。
再配合导出表面ESP映射数据：选“Export electron density and ESP on vdW surface to cube file”，导出.cube文件，再用VMD导入，用Color Scale标记ESP数值，就能得到那张经典的分子表面着色图。图的配色建议：红=负、蓝=正，这是领域内的通用惯例。
3.3 ELF/LOL实空间函数分析实操
做键的性质分析时，ELF（电子定域化函数）是特别直观的工具。在主菜单中选择实空间函数分析，再选电子定域化函数（ELF），设置好网格精度，程序会生成一个.cube文件。用VMD打开，调整等值面值（通常从0.7-0.8开始），就能看到电子在原子间区域的定域化特征，这对判断共价键、金属键、离子键的差异很有帮助。
网格精度设置是我特别想提醒的：默认网格有时候偏粗，画出来的等值面会有锯齿，做图不够光滑。我一般把格点间距调到0.05 Bohr左右，虽然需要的时间长一些，但出图质量好很多。对中大型分子（比如几百个原子），可以在关键区域缩小范围来提高效率，不需要全局都用细网格。
3.4 批处理玩法
Multiwfn有一个-silent模式，可以完全无人值守运行。把需要执行的命令按顺序写进文件，然后：
./Multiwfn water.fchk < input.txt > output.txt这样就能把整个分析流程脚本化。我做IGMH对比分析的时候，就是一次构建一批脚本，循环处理几十个分子构型，输出结果全部统一格式，再用自动化脚本提取关键数据。第一次配好之后，后面相同类型的分析只需要改改文件名，效率提升巨大。
4. 常见问题与排查技巧实录
4.1 文件格式导致的“读取失败”
Multiwfn对输入文件的版本兼容性总体是做得不错的，但偶尔还是会出现读取失败的情况。最典型的，Gaussian的.fchk文件如果包含非常规的泛函定义（比如用了自定义泛函），或者文件路径里有中文空格，可能导致读取异常。解决办法：先跑一遍测试小分子确认程序正常，然后再排查文件本身。
另一个高频坑是：直接从Windows拖文件到Linux终端，Windows格式的换行符（CRLF）可能会导致Multiwfn读取文件时出错。用dos2unix water.fchk转换一下再载入，基本就能解决。
提示：如果程序提示基函数信息不全（比如缺少
Shell types信息），检查一下是否用了pop=full。Gaussian默认可能不输出某些轨道信息，加上pop=full可以保证输出完整波函数信息。
4.2 泛函和基组怎么选才能“分析得准”
这是新手最容易搞混的一件事：计算级别直接决定了波函数质量，波函数质量决定分析结果可靠性。做定性分析（比如判断成键/反键、弱相互作用的有无），B3LYP/6-31G*级别已经够用；但如果你想定量比较不同体系的某个指标（比如ESP极值、键级数值），建议至少用到B3LYP/def2-TZVP或者ωB97XD/def2-TZVP级别。
这里有一个和Multiwfn本身没有直接关系、但困扰很多人的问题：Mulliken布居分析对基组的依赖非常大，同一个分子用6-31G*和6-311+G**算出来的Mulliken电荷可能差出0.3以上。我在自己的实践里，如果要做电荷分析并和文献对比，优先选择Hirshfeld电荷或ADCH电荷，这两个方法的结果对基组变化相对不敏感，更稳健。
4.3 画图软件和显示问题
Multiwfn自己不带三维可视化功能，它生成.cube文件或.pdb文件交给VMD、PyMOL等工具。VMD是我最常用的，但VMD的默认设置画出来的等值面可能需要手动调参数。记得在Graphics → Representations里把Drawing Method设为Isosurface，选择对应字段，再调整Isovalue，这样才能显示出你想要的等值面。
另外，Multiwfn导出的.txt格式的静电势极值点列表，直接复制到Excel里做表非常方便。不过要注意分隔符，有时候是空格分隔，有时候是制表符，如果你粘贴后数据错列，试试点选“文本分列”功能。
4.4 性能优化与超大体系
如果做的是上千个原子的体系（比如蛋白质片段、MOF模型），Multiwfn的实空间函数计算会非常耗时。这里有几个优化思路：
- 缩短格点范围，只计算你感兴趣的区域，而不是整个盒子；
- 降低格点精度，先跑出一个粗略结果确认趋势，最后做图时再提高精度；
- 利用对称性。如果分子有对称性，很多计算量可以简化，Multiwfn里有多极矩加速选项可以开启。
我做过一个约800个原子的体系，直接默认网格计算ELF，几乎跑了一个通宵。后来把计算范围裁剪到中心活性位点的10 Å范围内，同样的分析不到半小时就完成了，图像质量也没有肉眼可见的差别。所以大体系一定要学会裁剪范围。
4.5 不同量化程序输出的差异
Multiwfn支持的程序很多，但不同程序输出的波函数在轨道相因子、轨道排序规则上存在差异，尤其是涉及周期性计算（如VASP、CP2K）的时候，要特别注意：这些程序输出的波函数通常不是严格的分子轨道，而可能是Bloch函数相关数据，Multiwfn虽然提供了相应接口，但适用分析种类是受限的。如果遇到分析选项报错或者结果明显异常，先去手册里查一下该程序的支持范围再动手。
5. 一些顺手的小技巧
最后再分享几个我平时用得比较多的小技巧，这些细节不一定写在官方FAQ里，但确实能省下不少时间。
第一个：Multiwfn支持从命令行直接传入选项，例如：
./Multiwfn water.fchk -ES如果想跳过交互式菜单，直接运行某个特定的分析模块，可以在官方手册里查对应功能的命令行参数。配合Shell脚本做批量分析非常顺手，尤其适合需要对几十个文件做相同分析的场景。
第二个：导出的.cube文件通常比较大，如果你只是做视觉检查，可以用cube文件精简工具，或者直接在VMD里用volmap重建一个低密度的网格来显示趋势，视觉差异不会太大，但文件加载速度会快很多。
第三个：Multiwfn的“定义片段”功能特别强大。计算片段间的电荷转移、片段间的相互作用能分解等，都要用到它。我第一次做分子间CT分析时，费了不少劲才搞懂片段划分的逻辑。具体操作是：主菜单找“定义片段”，把原子列表按你需要的分组填进去，然后程序就会按片段输出各类数据。这个功能对弱相互作用体系（如氢键二聚体）很有用，建议认真读一下手册对应章节。
第四个：Multiwfn输出的大量数据，如果嫌手动记录麻烦，可以用tee命令或者直接重定向保存：
./Multiwfn water.fchk | tee result.log这样终端显示和文件保存同时进行，回头写论文的时候翻result.log比重新跑一遍省事太多。
Multiwfn这个工具，每次研究到一个新方向，需要一种新的分析方法时，我第一反应都是先去它的功能列表里搜一下。绝大多数时候，它都已经实现了，而且实现得足够好。如果你也刚接触波函数分析，或者正被某个分析需求卡住，不妨打开Multiwfn的手册，找到对应章节，按着例子跑一遍。它的官方手册写得很认真，绝大多数坑都能在里面找到答案。用顺手之后你就会发现，波函数分析不再是挡在研究面前的门槛，而是一把顺手好用的钥匙。
本文还有配套的精品资源，点击获取
