基于Multiwfn的分子机器学习助手：MultiwfnMLhelper

yumingsuxmu

近期，自己的工作里面需要快速提取Multiwfn分子描述符进行机器学习分析，我编写了一个名为“基于Multiwfn的分子机器学习助手”的Python脚本工具MultiwfnMLhelper。它的使用非常简便，给大家分享一下

【原料】任意能产生高斯fchk文件的任务即可，产生后将所有数据集中fchk文件放在某文件夹即可

【使用方法】
自己安装相应python包并运行MultiwfnMLhelper.py.建议将文件清单中中的内容全部复制到包含高斯输出的多个.fchk文件的文件夹中，运行MultiwfnMLhelper.py，选择当前文件夹的路径，然后点击确定，程序就会自动一键生成Multiwfn分子描述符,保存到一个名为 "Multwfn_analysis_feature_matrix3.csv" 的 CSV 文件中。

【文件清单】
(1) Multiwfn_helper.py：主要的Python脚本，用于提取Multiwfn描述符
(2) batchspec.bat：用于运行Multiwfn分析的批处理文件
(3) Step_1.txt1：Multiwfn分析时需要的参数文件
(4) settings.ini：务必isilent= 1(5) Multiwfn ：记得把最新的Multiwfn加到文件夹里

【参考】
http://sobereva.com/601
http://sobereva.com/612

【如果有什么建议或者感受，请告诉我】
（本来生成了一个exe文件，但受限于文件太大32M，没法上传了，有需要的可以私聊我）

paperz

yumingsuxmu 发表于 2023-4-12 10:11
多谢社长提醒，windows可直接运行程序.exe链接如下：
链接：https://pan.baidu.com/s/12mduHZQPwGR4-cq--Q ...

这些参数的单位都是什么呢？有的参数可能有多个单位，怎么确认输出的是哪个单位的数值？

ffcc

非洲黑炭 发表于 2025-5-17 14:17
你好 ！想请问下计算后 只有txt输出文件，没有描述符信息汇总是为啥

请问解决了吗

ffcc

lemon2223 发表于 2024-1-2 23:13
**** 作者被禁止或删除 内容自动屏蔽 ****

请问是怎么解决的

paperz

张晓婷 发表于 2023-4-17 21:12
你好，我下载了网盘链接文件，里面含有的和一楼发的文件清单是不一样的里面有两个exe，我运行了第二个exe ...

可能跟电脑硬件有关系，硬件和软件不适配，我换了一个电脑就可以了

paperz

张晓婷 发表于 2023-4-17 21:12
你好，我下载了网盘链接文件，里面含有的和一楼发的文件清单是不一样的里面有两个exe，我运行了第二个exe ...

你好，我也遇到这个问题，fchk和操作都没问题，就是生成空的txt，请问你解决了吗

Eternal

楼主提供的Step_1.txt1是不是可以完善一下，Step_1.txt1末尾改为

2
4
0
-1
-1
q
不改的话会提示forrtl: severe (24): end-of-file during read, unit -4, file CONIN$

〇66

sobereva 发表于 2025-10-20 00:16
把我前面提到的601博文完整看一遍，99%的描述符是什么含义、怎么算的，你应该就能了解了，剩下的极少数糊 ...

好的老师 我完整看一遍 谢谢

sobereva

〇66 发表于 2025-10-19 11:35
是的老师，我是在这个帖子http://sobereva.com/426和手册里ctrl+F搜索“Mol_Size_2”这个关键词，但是没 ...

把我前面提到的601博文完整看一遍，99%的描述符是什么含义、怎么算的，你应该就能了解了，剩下的极少数糊涂的可以发帖问楼主

〇66

Uus/pMeC6H4-/キ 发表于 2025-10-19 12:01
完全找错地方了，这种变量名并不是社长的multiwfn直接输出的而是楼主自定义于脚本的，得到MultiwfnMLhelp ...

原来是这样，感谢你耐心详细的解答，我目前还没有python基础，那后续找别的描述符的物理意义感觉也是怪麻烦的==只能一个帖子一个帖子的盯对了。

Uus/pMeC6H4-/キ

〇66 发表于 2025-10-19 11:35
是的老师，我是在这个帖子http://sobereva.com/426和手册里ctrl+F搜索“Mol_Size_2”这个关键词，但是没 ...

完全找错地方了，这种变量名并不是社长的multiwfn直接输出的而是楼主自定义于脚本的，得到MultiwfnMLhelper.py里面找并结合python知识理解。实际上第1页就有列表但格式不工整且解释不完善。这个关键词对应的部分代码是

1.     for line in lines:
   
2.     [omitted]
   
3.         elif ' Length of the three sides:' in line:
   
4.             mol_size = list(map(float, re.findall(r'([\d.]+)', line)))
   
5.             sorted_mol_size = sorted(mol_size)
   
6.             mol_size_short = sorted_mol_size[0]
   
7.             mol_size_2 = sorted_mol_size[1]
   
8.             mol_size_l = sorted_mol_size[2]
   
9.             length_ratio = mol_size_l / sum(mol_size)
   
10.             len_div_diameter = mol_size_l / (2 * mol_radius)

复制代码

这相当于定位到multiwfn输出文件中有“Length of the three sides:”的行，存入并排序该行所有浮点数（一共3个）；由于sorted()函数按升序排列，故最小的边长mol_size_short在第1个元素sorted_mol_size[0]，最大的边长mol_size_l在第3个元素sorted_mol_size[2]，中间的边长mol_size_2在第2个元素sorted_mol_size[1]。（变量命名实在好乱，又缺乏清晰文档，所以出这种问题并不怪用户。）

〇66

sobereva 发表于 2025-10-19 05:39
我的博文里有，手册里也有
使用Multiwfn计算分子的长宽高以及显示分子的主轴
http://sobereva.com/426 ...

是的老师，我是在这个帖子http://sobereva.com/426和手册里ctrl+F搜索“Mol_Size_2”这个关键词，但是没有找到对应的地方。想着后续其他描述符在不知道物理意义的情况下，查找起来可能不是很方便。谢谢您的回复。

sobereva

〇66 发表于 2025-10-18 17:29
感谢卢老师和楼上老师的整理，我这里还想请教一下，每个描述符都是英文字母，可能不太清楚他的物理意义是 ...

我的博文里有，手册里也有
使用Multiwfn计算分子的长宽高以及显示分子的主轴
http://sobereva.com/426（http://bbs.keinsci.com/thread-10226-1-1.html）

下文也提到了
Multiwfn可以计算的分子描述符一览
http://sobereva.com/601（http://bbs.keinsci.com/thread-23729-1-1.html）

Multiwfn能计算分子的长、宽、高，见《使用Multiwfn计算分子的长宽高以及显示分子的主轴》（http://sobereva.com/426）。

〇66

sobereva 发表于 2023-4-12 11:42
我已把此文链接加入此帖了
Multiwfn可以计算的分子描述符一览
http://sobereva.com/601（http://bbs.kein ...

感谢卢老师和楼上老师的整理，我这里还想请教一下，每个描述符都是英文字母，可能不太清楚他的物理意义是什么，想要具体了解的话，有没有一个简便的形式可以看到他对应解释的链接。比如16. Mol_Size_2这个描述符，我在这个帖子中对应的每个帖子和multiwfn手册中进行检索还是没找到，查找起来也比较麻烦一点。问问有没有什么简便的方式能一下知道对应描述符多物理意义或者解释，谢谢老师。

facia

你好，请问这个exe文件是不是没有办法直接生成CSV表格呀？我使用exe文件只生成了txt文件