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

博文

qm/mm 小分子计算电荷

已有 4580 次阅读 2015-8-20 14:54 |系统分类:科研笔记

antechamber 功能比较弱,只能处理H,C,N,O,S,P and halogens这些元素。

如果分子中含有其他元素如硼,可通过multiwfn计算resp电荷。

然后再转成力场文件。



测试后发现

g09 A 版本 不能产生gesp 文件

g09 D 版本 可以长生gesp文件

-----------

antechamber in  amber12 不能直接识别高斯的log文件,所以输入文件必须是gesp文件

antechamber in amber14 可以直接识别高斯的log文件



antechamber -i ../*log -fi gout -o *_lig.prep -fo prepi -c resp

-----------

gaussian中的关键词关键词是pop=mk  iop(6/33=2,6/42=6)  这样可以直接使用log文件。 产生log文件,可被amber14中的命令使用。

G09D 中的关键词是pop=mk  iop(6/50=1);文件的末尾(一般在坐标后面空一行写上  )xxx.gesp,   产生log文件以及gesp文件,其gesp文件可被amber12中的命令使用。


antechamber -i xxx.gesp -fi gesp -o xxx.prep -fo prepc -c resp



------

        3)antechamber -i ../*log -fi gout -o *_lig.prep -fo prepi -c resp

4)parmchk  -i *_lig.prep -f prepi -o ligand.frcmod -a y

#  选项和值之间要有空格 -f   prepi




- Check Unusual Elements --
Warning: Unusual element (B) for atom (ID: 7, Name: B1).
GAFF does not have sufficient parameters for molecules having unusual
         elements (those other than H,C,N,O,S,P and halogens).
         To ensure antechamber works properly, one may need to designate
         bond types for bonds involved with unusual elements.
         To do so, simply freeze the bond types by appending "F" or "f"
         to the corresponding bond types in ac or mol2 files.




https://wap.sciencenet.cn/blog-950202-914487.html

上一篇:bpython 更好的python shell
下一篇:QMMM计算简单流程
收藏 IP: 202.127.19.*| 热度|

0

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

数据加载中...

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

GMT+8, 2024-5-3 14:47

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部