sobereva 发表于 2022-4-19 20:55
Multiwfn 3.7不支持周期性
目前的Multiwfn支持考虑周期性，然而你的xyz文件里的盒子信息又是完全错的，Mul ...

原来是这样的，了解一二了。谢谢社长指教。

snljty 发表于 2022-4-19 20:50
你xyz第2行那个晶胞参数明显有问题啊。如果去掉当孤立体系跑，中等格点质量，8核电脑5秒钟算完，用2022.3.2 ...

老哥666，把晶胞参数删掉之后，IGM分析果然就正常了。

snljty 发表于 2022-4-19 20:50
你xyz第2行那个晶胞参数明显有问题啊。如果去掉当孤立体系跑，中等格点质量，8核电脑5秒钟算完，用2022.3.2 ...

我用 10 Set box of grid data visually using a GUI window 这个菜单定义的计算范围和格点密度，按说应该和晶胞关系不大了呀。如果正常的话，确实是两三分钟。

Multiwfn 3.7不支持周期性
目前的Multiwfn支持考虑周期性，然而你的xyz文件里的盒子信息又是完全错的，Multiwfn自然没法正常计算。进主功能0，把盒子显示出来，一看就知道完全离谱


不想当周期性来考虑就根本别在xyz里写盒子信息

你xyz第2行那个晶胞参数明显有问题啊。如果去掉当孤立体系跑，中等格点质量，8核电脑5秒钟算完，用2022.3.21更新的Multiwfn 3.8 (dev)。如果不去掉，我是不知道你晶胞用的多大...不过8核电脑应该3分钟左右能算完。

sobereva 发表于 2022-4-19 17:22
把输入的具体命令贴出来，便于别人重复

3.7不代表稳定。现在也强烈不建议用老掉牙的3.7

谢谢社长关心。
操作步骤和命令是这样的。
启动Multiwfn.exe （当前最新版3.8），把structure_1238_atoms.xyz 文件拖到 Multiwfn打开。
20  //弱相互作用可视化分析
10  //IGM analysis
12 //定义12个片段  
817-952 //第一个片段里的序号，对应xyz文件转化为pdb文件，使用Gauss View查看得到的原子序号
1-136  //第二个片段里的序号
545-680  //第三个片段里的序号
273-408  //第四个片段里的序号
409-544  //第五个片段里的序号
681-816  //第六个片段里的序号
137-272  //第七个片段里的序号
953-1088  //第八个片段里的序号
1089-1160  //第九个片段里的序号
1161-1232  //第十个片段里的序号
1236-1238  //第十一个片段里的序号
1233-1235  //第十二个片段里的序号
10    更改Grid spacing (Bohr)至0.17 //Grid spacing (Bohr)设置为其他值时，比如0.25 和0.7，也会出现相同的Bug
点击右上角RETURN 返回Multiwfn窗口。
然后就是几个小时的等待。计算完成后
3   //导出cub文件
然后使用VMD 或Pymol作图，都只能看到结构，看不到任何等值面。

后改用Multiwfn 3.7版本进行相同的操作，则一切顺利正常，且整个过程也不用几个小时，几分钟就OK了。

附：
vmd作图方法为（参考 http://sobereva.com/407）。将dg_inter.cub和sl2r.cub都拷到VMD目录下，把Multiwfn自带的examples目录下的VMD作图脚本IGM_inter.vmd也拷到VMD目录下，然后启动VMD，在命令行窗口输入source IGM_inter.vmd。
Pymol作图方法为（参考 http://bbs.keinsci.com/thread-15756-1-1.html。）
将 Draw_RDG_with_PyMOL.ini 中isovalue设置为0.005。打开Draw_RDG_with_PyMOL_MultipleCube.bat ，依次加载dg_inter.cub, sl2r.cub 和 sl2r.structure.mol2，自动打开PyMol

把输入的具体命令贴出来，便于别人重复

3.7不代表稳定。现在也强烈不建议用老掉牙的3.7
最新版的bug总数总是显著少于上一个正式版

问题解决了。使用旧版 Multiwfn3.7稳定版 计算正常。