“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register

使用Multiwfn+ORCA自动批量计算静电势与偶极矩的脚本

查看数: 5068 | 评论数: 17 | 收藏 Add to favorites 12
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-5-7 03:23

正文摘要:

本帖最后由 Stardust0831 于 2025-5-9 20:08 编辑 静电势(ESP)和偶极矩作为非常常用的数据,批量计算的需求很旺盛。 其中,静电势对于考察分子间静电相互作用、预测反应位点、预测分子性质等方面有重要意义, ...

回复 Reply

yubozou@163.com 发表于 Post on 2026-2-23 10:35:34
Stardust0831 发表于 2026-2-22 02:11
tachyon_WIN64.exe是在vmd安装文件夹下的,需要确保1.使用的确实是64位1.9.3的vmd,而非官网上那个32位的 ...

感谢感谢,原来我的是32位的
Stardust0831 发表于 Post on 2026-2-22 02:11:30
yubozou@163.com 发表于 2026-2-18 16:53
你好,我用上面新的脚本运行后vmd窗口正常,但是.bmp格式的文件没办法输出,窗口提示couldn't execute "tac ...

tachyon_WIN64.exe是在vmd安装文件夹下的,需要确保1.使用的确实是64位1.9.3的vmd,而非官网上那个32位的版本2.已经把vmd所在目录加到了系统环境变量里

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
wal + 5

查看全部评分 View all ratings

yubozou@163.com 发表于 Post on 2026-2-18 16:53:12
你好,我用上面新的脚本运行后vmd窗口正常,但是.bmp格式的文件没办法输出,窗口提示couldn't execute "tachyon_WIN64.exe": no such file or directory,想问下要怎么操作,是我哪里有问题?谢谢

202602181648102745..png (17.34 KB, 下载次数 Times of downloads: 10)

cmd命令窗口

cmd命令窗口

202602181648423955..png (19.13 KB, 下载次数 Times of downloads: 9)

esp文件夹内容

esp文件夹内容
CM234 发表于 Post on 2026-2-12 18:45:49
Stardust0831 发表于 2026-1-29 16:51
当前文件夹下有哪些文件呢?

感谢回复,已经排除了,我有一个路径写错了。
Stardust0831 发表于 Post on 2026-1-29 16:51:35
CM234 发表于 2026-1-23 17:32
请问我将vmd,orca和Multiwfn都加入系统变量,也修改了ESP_orca.sh文件中的路径,运行后出现如下报错,请问 ...

当前文件夹下有哪些文件呢?
CM234 发表于 Post on 2026-1-23 17:32:45
本帖最后由 CM234 于 2026-1-23 18:22 编辑

请问我将vmd,orca和Multiwfn都加入系统变量,也修改了ESP_orca.sh文件中的路径,运行后出现如下报错,请问要怎么排除呢?附上我的输入文件

屏幕截图 2026-01-23 182004.png (365.37 KB, 下载次数 Times of downloads: 24)

屏幕截图 2026-01-23 182004.png

Achun-B-Al.gjf

6.2 KB, 下载次数 Times of downloads: 0

wwwwwt 发表于 Post on 2025-8-15 08:39:00
Stardust0831 发表于 2025-8-14 21:29
g16用不了GFN2-xTB这种半经验级别,所以我自己的工作流不会用g16。
简单给你讲一下怎么批量运行:
(假 ...

感谢感谢!
wwwwwt 发表于 Post on 2025-8-15 08:38:51
Stardust0831 发表于 2025-8-14 21:07
gjf是gaussian的格式,我会用gjf当输入是因为我建模的主力软件是gview。
你文件里写的事实上是orca的格 ...

感谢回复!
Stardust0831 发表于 Post on 2025-8-14 21:29:28
本帖最后由 Stardust0831 于 2025-9-24 13:53 编辑
wwwwwt 发表于 2025-8-13 11:08
请问有针对Gaussian16 Linux版本的批量计算静电势脚本吗?

g16用不了GFN2-xTB这种半经验级别,所以我自己的工作流不会用g16。
简单给你讲一下怎么批量运行:
(假设gaussian已经被正确安装,在环境变量里了)
1.批量执行gjf文件
  1. for inf in *.gjf
  2. do
  3. echo Running ${inf} ...
  4. g16 < ${inf} > ${inf//gjf/out}
  5. done
复制代码
这5行单独当一个sh脚本,可以执行当前目录下所有以.gjf结尾的文件。

2.批量生成波函数文件:
  1. for inf in *.chk
  2. do
  3. formchk ${inf}
  4. done
复制代码
3.生成用于绘图的格点文件:
  1. for inf in *.fchk
  2. do
  3. echo Export cub file from ${inf} ...
  4. Multiwfn.exe $inf << EOF > log.txt
  5. 5
  6. 1
  7. 3
  8. 2
  9. 0
  10. 5
  11. 12
  12. 1
  13. 2
  14. 0
  15. r
  16. totesp.cub
  17. 13
  18. 11
  19. 5
  20. 627.5095
  21. 0
  22. totesp.cub
  23. -1
  24. q
  25. EOF
  26. mv -f density.cub ${inf//.fchk/_density.cub}
  27. mv -f totesp.cub ${inf//.fchk/_ESP.cub}
  28. echo
  29. done
复制代码
4.绘制esp,此时基本上类似《使用Multiwfn+VMD快速地绘制静电势着色的分子范德华表面图和分子间穿透图(含视频)》里的流程,只是.vmd文件要换成这个:
ESP_stardust0831.vmd (2.53 KB, 下载次数 Times of downloads: 4)
和原版的差异是cub文件命名是按“分子名_ESP.cub”这种格式的,可以直接看cub认出是哪个分子,并且有我在二楼提到的一些自定义函数,比如“xr 文件名”直接按特定文件名导出并渲染图片等。


Stardust0831 发表于 Post on 2025-8-14 21:07:49
本帖最后由 Stardust0831 于 2025-8-14 21:11 编辑
wwwwwt 发表于 2025-8-14 09:47
楼主你好,我在运行esp_sh时,一直出现输入文件无效的错误,附件是我的输入文件,请问是哪里出错了么?

gjf是gaussian的格式,我会用gjf当输入是因为我建模的主力软件是gview。
你文件里写的事实上是orca的格式和关键词(并且关键词是过时的),这显然和gjf这个后缀不符,multiwfn当gaussian输入文件来读取自然认不出来。
把脚本里的
  1. File_extension="gjf"
复制代码
改成
  1. File_extension="in"
复制代码
并且把这个orca格式的文件的后缀改成.in,这样multiwfn就可以读取了。
wwwwwt 发表于 Post on 2025-8-14 09:47:05
楼主你好,我在运行esp_sh时,一直出现输入文件无效的错误,附件是我的输入文件,请问是哪里出错了么?

test.gjf

274 Bytes, 下载次数 Times of downloads: 6

wwwwwt 发表于 Post on 2025-8-14 09:43:38
wal 发表于 2025-8-13 14:59
显示在哪个原子/基团上有点难吧。顶多告诉你index是多少,具体属于哪个原子/基团得自己判断

感谢
wal 发表于 Post on 2025-8-13 14:59:54
wwwwwt 发表于 2025-8-13 14:31
请问显示静电势极大值和极小值时,是否能显示其极大极小值在哪个原子或基团上么?

显示在哪个原子/基团上有点难吧。顶多告诉你index是多少,具体属于哪个原子/基团得自己判断
wwwwwt 发表于 Post on 2025-8-13 14:31:29
请问显示静电势极大值和极小值时,是否能显示其极大极小值在哪个原子或基团上么?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-4-13 14:12 , Processed in 0.184971 second(s), 27 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list