计算化学公社

 找回密码 Forget password
 注册 Register

求助在top文件中添加其他小分子力场时报错Invalid order for directive atomtypes

查看数: 4524 | 评论数: 7 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-1-16 17:59

正文摘要:

本帖最后由 sj12138 于 2022-1-17 13:29 编辑 前辈们好,我是做的水合物分子动力学的。我在ATB网站上下载的氢气的itp文件以及其自带的力场,将H的itp文件include到top文件中后报错。 我的命令是: gmx_mpi_d gro ...

回复 Reply

sobereva 发表于 Post on 2022-2-23 18:46:35
494085639 发表于 2022-2-23 15:29
老师可以请问下include到主top文件时的格式是什么样的吗

语义不明
搞清楚gromacs的各个字段都该有什么内容,该以什么顺序出现,自然就明白了
PS:我在GROMACS培训里讲得特别透彻 http://www.keinsci.com/workshop/KGMX_content.html,搞不明白这些基本知识的话最好系统性学学
494085639 发表于 Post on 2022-2-23 15:29:14
sobereva 发表于 2022-1-17 17:12
这种体系完全不需要用pdb2gmx
直接把小分子itp文件include到主top文件就完了,折腾rtp纯粹是自找麻烦

老师可以请问下include到主top文件时的格式是什么样的吗
sobereva 发表于 Post on 2022-1-17 17:12:23
sj12138 发表于 2022-1-17 14:07
老师,我在54A7的rtp文件中加上了19P的残基和原子之后采用gmx_mpi_d pdb2gmx -f II-zj-222.pdb -o II-zj- ...

这种体系完全不需要用pdb2gmx
直接把小分子itp文件include到主top文件就完了,折腾rtp纯粹是自找麻烦
sj12138 发表于 Post on 2022-1-17 14:12:42
牧生 发表于 2022-1-16 18:55
看我以前犯的相同错误
http://bbs.keinsci.com/thread-19761-1-1.html
其次,力场不建议混用,opls和grom ...

帖子我看了,我尝试之后还是没有跑通才发帖请教的。由于我的H.itp并没有atomtype的字段,所以我在ffnonbonded.itp中的内容贴到了top中但是依然报最开始的错误
sj12138 发表于 Post on 2022-1-17 14:07:53
sobereva 发表于 2022-1-17 07:16
怎么能两个力场的forcefield.itp都include,没有基本常识

ATB提供的修改版54A7力场包里又不是没有spce.i ...

老师,我在54A7的rtp文件中加上了19P的残基和原子之后采用gmx_mpi_d pdb2gmx -f II-zj-222.pdb -o II-zj-222.gro -ff gromos54a7_atb生成top文件并grompp之后跑通了,谢谢sob老师,是我之前没想到不需要oplsaa。

但是老师我想问一下,因为我的小分子是氢气很简单所以在rtp文件中可以很容易补充,但是如果是一个比较复杂的分子应该如何在rtp中确定需要补充的内容那?

谢谢sob老师
sobereva 发表于 Post on 2022-1-17 07:16:07
怎么能两个力场的forcefield.itp都include,没有基本常识

ATB提供的修改版54A7力场包里又不是没有spce.itp,include OPLS-AA力场下的spc.itp作甚
牧生 发表于 Post on 2022-1-16 18:55:54
看我以前犯的相同错误
http://bbs.keinsci.com/thread-19761-1-1.html
其次,力场不建议混用,opls和gromos不能混用

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

GMT+8, 2026-2-23 21:14 , Processed in 0.187769 second(s), 25 queries , Gzip On.

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