计算化学公社

标题: 将当前目录下所有gbw转化为同名fch文件的bash脚本 [打印本页]

作者
Author: sobereva 时间: 2020-1-19 20:16
标题: 将当前目录下所有gbw转化为同名fch文件的bash脚本
运行附件里的bash脚本，即可将当前目录下所有gbw转化为同名的fch或fchk。
(, 下载次数 Times of downloads: 127)

机子里必须已经正常安装ORCA和Multiwfn，从而能直接输入orca_2mkl和Multiwfn命令启动。

有些人可能用得着此脚本。比如Multiwfn可以很容易地计算RESP电荷，看《RESP拟合静电势电荷的原理以及在Multiwfn中的计算》（http://sobereva.com/441），并且其中可以调用cubegen降低静电势计算耗时，但此时只支持fch/fchk作为构象列表里设的文件的格式；而对于ORCA做量化计算的用户，就可以先用这个脚本把所有构象的gbw文件一键转成fch，再给Multiwfn用于计算考虑多构象的RESP电荷的目的。

作者
Author: zsu007 时间: 2020-1-20 10:04
谢谢社长的分享！

作者
Author: alexanderqi 时间: 2020-2-4 12:05
多谢社长分享！非常实用！

作者
Author: 447951397 时间: 2020-6-17 21:55
社长，做完频率计算后，能不能将gbw文件中的Hessian矩阵一并转换到fch文件中，因为Compliance算柔性力常数时必须用Gaussian的fch才能读到

作者
Author: sobereva 时间: 2020-6-21 03:10

447951397 发表于 2020-6-17 21:55
社长，做完频率计算后，能不能将gbw文件中的Hessian矩阵一并转换到fch文件中，因为Compliance算柔性力常数 ...

Multiwfn本身读不了gbw，而且Multiwfn是先把gbw用orca_2mkl弄成molden文件才转化成的fch，而molden格式自己并没有记录Hessian
所以还是自行把.hess文件里的内容拷到fch里吧

作者
Author: 447951397 时间: 2020-6-21 10:32

sobereva 发表于 2020-6-21 03:10
Multiwfn本身读不了gbw，而且Multiwfn是先把gbw用orca_2mkl弄成molden文件才转化成的fch，而molden格式自 ...

好的，谢谢社长

欢迎光临计算化学公社 (http://bbs.keinsci.com/)