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

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

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

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

sobereva

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

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

〇66

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

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

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$