YuYuqi的个人博客分享 http://blog.sciencenet.cn/u/YuYuqi

博文

Tao package, QM/MM文件准备时辅助工具

已有 3225 次阅读 2018-1-23 15:56 |系统分类:科研笔记

原文教程:http://chem.wayne.edu/schlegel/Software/oniomtoolTAO/TAOtutorial.html

做QM/MM时,有时结构的原子太多,用gaussview处理时,会出现长时间卡顿,着实恼人,Tao是用命令行操作,可以与gaussview配合使用,提高效率。

程序功能:

1.将PDB文件转成gjf文件(pdb2oniom --help)

2.将系统的中一部分截取(pdbcore)

3.检查不正确的连接方式,前提是gjf里写了connectivity信息(checkconnect)

4.计算高层和低层所带电荷(chargesum)

5. 使用gaussian报错log文件自动补全缺失的参数如键长和键角等(parmlookup),但缺失的原子范德华半径需要手动查找补全

 例: parmlookup -g mmp2_partial_r_06.log -o mmp2_partial_r_06_parm.txt

mmp2_partial_r_06_parm.txt文件即可复制到gjf文件中的末尾。注意,有些无法完全匹配amber力场文件的原子,不能补全,参数以0.000表示,需要手动补全。

amber力场文件里的原子类型多种多样,其含义参考

http://ambermd.org/antechamber/gaff.html#atomtype    (GAFF力场小写原子)

http://www.chem.cmu.edu/courses/09-560/docs/msi/ffbsim/B_AtomTypes.html (AMBER力场大写原子)

补全参数时应搞清楚各原子杂化方式,对应于相似的原子类型。







                   



https://blog.sciencenet.cn/blog-3375268-1096395.html


下一篇:选择小分子附近的完整残基
收藏 IP: 202.127.19.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-4-20 07:42

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部