请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 472|回复: 3

[ORCA] ORCA QMMM 测试

[复制链接]

4

帖子

0

威望

61

eV
积分
65

Level 2 能力者

发表于 2019-8-18 21:33:33 | 显示全部楼层 |阅读模式
在最近发布的ORCA 4.2.0中新增了对QMMM的支持,支持读取top/pdb文件。
我一直在使用Gaussian做QMMM,这几个月在与chemshell的tcl做斗争,目前我能在Chemshell/orca这一套工具中进行,opt,PES(potential energy scan),TS opt,还差一个frequency。若是ORCA可以读取top+pdb,那真的是极好,这几天测试一下,分享一下其中的坑与我似乎已经能填上的坑,以此抛砖引玉。

1关于ORCAFF.prms

首先,orca说他能读取top。实际上不是所有的格式的top,能直接读取的仅限于charmm的top,并且是VMD(或者QWikMD)生成的这种charmm.psf文件。并不支持amber的prmtop。本人测试,由amber转的charmm.psf也是不可以的。好在另一个选择是生成orca自己的top文件——ORCAFF.prms。于是我打算自己写这个top,但是发现文档里面没有写ORCAFF.prms的格式信息,我只好根据自己的理解模仿orca_mm的输出,如果又那位热心的朋友发现了这部分的官方文档,请告诉我,在这里先谢了。

这里添加了我用于做这个转换的jupyter notebook。其中使用了parmed包,原则上经过少许的修改就可以从charmm/gmx/amber转ORCAFF.prms。附件里也附带了3个测试top文件,欢迎大家下载去玩玩。

2第二个坑在于index
在输入文件里,定义QMatom/ActiveAtom是从0开始数的,但是ORCAFF.prms拓扑文件确实从1开始数的。

3
我有些怀疑在进行QMMM的时候,不是所有的关键词都是原来的意思,比如
%geom
    MaxIter 50          # max. number of geometry iterations
    coordsys redundant  # New redundant internal coords
    Convergence loose   # Default is normal
end
但是这个任务却并没有在50步优化的时候停止,也许是我的理解有误,我对orca的使用其实也才刚刚开始没几个月。

4输出
orca在运行QMMM的时候会把输出信息分成好几个文件,本人觉得这样比gaussian要高明许多。Gaussian把大量的文件输出在一个log文件里,注释让人头秃。ORCA会把轨迹输出为.dcd文件,scf单独一个文件,gradient单独一个文件,能量/deltaE/RMS(Grad)/Max(Grad)/Max(Step)单独输出一个.csv文件。但是这个CSV文件居然不是用都好分割而是用分号;分隔,我估计这是一个BUG。

我的一点点使用体验就是这样了



评分

参与人数 2eV +14 收起 理由
puzhongji + 4
sobereva + 10

查看全部评分

2万

帖子

25

威望

2万

eV
积分
45911

管理员

公社社长

发表于 2019-8-19 10:25:10 | 显示全部楼层
帖子里提到的jupyter notebook没有上传?(虽然我在ORCA论坛里已经看到了)

对于不是着急用QM/MM的人,我倾向于等ORCA这个功能再成熟、完善一些再上手。刚出来的,选项一大堆的较大功能难免很不完善。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

4

帖子

0

威望

61

eV
积分
65

Level 2 能力者

 楼主| 发表于 2019-8-19 16:09:20 | 显示全部楼层
本帖最后由 惠成功QUB 于 2019-8-22 14:22 编辑

第一次发帖,还不熟悉论坛的操作
在这里补上我的jupyter notebook。
确实,ORCA论坛里也是我发的,是一样的内容。我也确实不是很有动力继续尝试ORCA的QMMM,毕竟我的chemshell看上去已经挺好用的了


PS 此处的脚本有误,其中的Atom部分有错误。更正的脚本在后续回复中


prmtop2ORCAff.zip

1.46 MB, 下载次数: 1

4

帖子

0

威望

61

eV
积分
65

Level 2 能力者

 楼主| 发表于 2019-8-22 14:20:11 | 显示全部楼层
本帖最后由 惠成功QUB 于 2019-8-22 14:22 编辑

更正一下我的jupyter notebook。本次修正了之前的Atom部分LJ项的错误

官方给出的ORCAFF.prms文件格式的解释看这个链接
https://orcaforum.kofo.mpg.de/viewtopic.php?f=11&t=5061#p21563

本人新的jupyter notebook,包含LJ项输出的更正,和增加angle UB项

prmtop2ORCAff-08-22.zip

1.53 MB, 下载次数: 3

jupyter notebook可以将Amber.prmtop/Charmm.psf装换位ORCAFF.prms(ORCA MM 拓扑)

评分

参与人数 1eV +5 收起 理由
sobereva + 5

查看全部评分

您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-9-21 09:36 , Processed in 0.155108 second(s), 28 queries .

快速回复 返回顶部 返回列表