| 查看: 2855 | 回复: 8 | |||
[求助]
topol文件不识别小分子 已有1人参与
|
|
最近刚刚学gromacs，按照网上的tutorial步骤走遇到一个问题就是我要添加我的小分子到pdb2gmx产生的conf.gro文件，我按照网上的tutorial添加需要的文件之后，还是不识别我的小分子，怎么办？我把topol文件里面的分子数手动修改之后，就出现下面的错误了 NOTE 1 [file topol.top, line 10673]: System has non-zero total charge: 6.000001e+00 Warning: atom name 1694 in topol.top and solv.gro does not match (CAA - AA) Warning: atom name 1695 in topol.top and solv.gro does not match (CAG - AG) Warning: atom name 1696 in topol.top and solv.gro does not match (CAH - AH) Warning: atom name 1697 in topol.top and solv.gro does not match (CAJ - AJ) Warning: atom name 1698 in topol.top and solv.gro does not match (CAF - AF) Warning: atom name 1699 in topol.top and solv.gro does not match (HAG - AG) Warning: atom name 1700 in topol.top and solv.gro does not match (CAD - AD) Warning: atom name 1701 in topol.top and solv.gro does not match (HAE - AE) Warning: atom name 1702 in topol.top and solv.gro does not match (CAC - AC) Warning: atom name 1703 in topol.top and solv.gro does not match (HAD - AD) Warning: atom name 1704 in topol.top and solv.gro does not match (CAE - AE) Warning: atom name 1705 in topol.top and solv.gro does not match (HAF - AF) Warning: atom name 1706 in topol.top and solv.gro does not match (CAI - AI) Warning: atom name 1707 in topol.top and solv.gro does not match (OAB - AB) Warning: atom name 1708 in topol.top and solv.gro does not match (HAL - AL) WARNING 1 [file topol.top, line 10673]: 15 non-matching atom names atom names from topol.top will be used atom names from solv.gro will be ignored NOTE 2 [file topol.top]: The largest charge group contains 11 atoms. Since atoms only see each other when the centers of geometry of the charge groups they belong to are within the cut-off distance, too large charge groups can lead to serious cut-off artifacts. For efficiency and accuracy, charge group should consist of a few atoms. For all-atom force fields use: CH3, CH2, CH, NH2, NH, OH, CO2, CO, etc. Analysing residue names: There are: 163 Protein residues There are: 1 Other residues There are: 10445 Water residues Analysing Protein... Analysing residues not classified as Protein/DNA/RNA/Water and splitting into groups... Number of degrees of freedom in T-Coupling group rest is 67791.00 Largest charge group radii for Van der Waals: 0.354, 0.287 nm Largest charge group radii for Coulomb: 0.354, 0.287 nm Calculating fourier grid dimensions for X Y Z Using a fourier grid of 72x72x72, spacing 0.109 0.109 0.109 Estimate for the relative computational load of the PME mesh part: 0.37 This run will generate roughly 383 Mb of data There were 2 notes There was 1 warning 这个怎么办？ |
回复此楼» 猜你喜欢
一志愿华中师范071000，325求调剂
已经有6人回复
-
070303一志愿西北大学学硕310找调剂
已经有6人回复
-
311求调剂
已经有4人回复
-
0703化学调剂
已经有8人回复
-
070303 总分349求调剂
已经有5人回复
-
309求调剂
已经有5人回复
-
070300化学学硕求调剂
已经有5人回复
-
277材料科学与工程080500求调剂
已经有3人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有5人回复
-
环境工程调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 分子进化
已经有60人回复
- 有人遇到过MS文件读写错误的情况吗？
已经有4人回复
- 求助蛋白纯化的高手
已经有5人回复
- 求可以说明主客体分子识别的实验
已经有10人回复
- 15个金币求助投哪个期刊杂志（分子识别方面）
已经有4人回复
- QMolSymmetry:分子对称性识别，调整，优化（C++源码）
已经有6人回复
- 使用towhee进行GEMC模拟相平衡的MC运动定义问题..
已经有15人回复
- 如何入门有机中的分子识别或者离子识别或者阴阳离子受体等方向大讨论
已经有4人回复
- 【求助】玻璃化转变温度的识别
已经有20人回复
- 【转贴】识别图片PDF文件中的文字【已搜索无重复】
已经有16人回复
- CTEX2.8编辑的文件TEX/LATEX2.0识别吗
已经有6人回复
- 【求助】namd和vmd就只能识别.pdb.psf.coor.vel等文件吗
已经有6人回复
- 专题调查Understanding Molecular Simulation: from Algorithm to Application
已经有39人回复
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
jiaoyixiong(金币+5): 鼓励交流 2011-12-30 22:28:14
蓝紫色玻璃心(金币+1): 2011-12-31 08:53:56
感谢参与，应助指数 +1
jiaoyixiong(金币+5): 鼓励交流 2011-12-30 22:28:14
蓝紫色玻璃心(金币+1): 2011-12-31 08:53:56
|
Warning存在的问题：你在top文件里面的分子的排例顺序与gro文件出现的顺序不一样，所以会这样。举个例子，如果你的体系有两种物种A与B.A有10个原子组成（A1,....A10），B有5个原子组成(B1,....B5). 在拓扑文件里面，若你是这样写的话： [ molecules ] A 10 B 10 那么，你的gro文件里面的顺序也必须是：先写10个A分子的坐标，然后再写10B分子的坐标。这样，才能一一对应。 如果你的gro文件是先写10B分子，然后再写10A分子的话，那么，前50个原子都会对不上号，即A1--B1,这样的mismatching的现像。假如这个时候，你忽略这个warning(-maxwarn 来忽略)，那么，就按top文件这个顺序走下去，显然会出现错误。 结论是：你再检查一下你的top与gro有没有一一对应。 |
2楼2011-12-30 21:39:53
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
蓝紫色玻璃心(金币+1): 2011-12-31 08:59:40
御剑江湖(金币+5): 谢谢 2011-12-31 17:51:48
感谢参与，应助指数 +1
蓝紫色玻璃心(金币+1): 2011-12-31 08:59:40
御剑江湖(金币+5): 谢谢 2011-12-31 17:51:48
|
There was 1 warning 这个warning可以忽略不管，因为： atom names from topol.top will be used atom names from solv.gro will be ignored 我在做模拟的时候，也遇到过类似的问题，这个warning忽略掉不会影响。在后续的模拟中，就不会提示了。 There were 2 notes 第一个notes的意思是你的体系没有电中性，你的 系统带有6e个电荷 因此你需要给系统添加 -6e个电荷 |
3楼2011-12-30 22:32:15
4楼2011-12-31 08:59:14
5楼2011-12-31 19:59:32
juanjuan0618
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 179.6
- 帖子: 94
- 在线: 23.6小时
- 虫号: 1160553
- 注册: 2010-12-02
- 专业: 实验地球化学和计算地球化
6楼2014-03-07 19:44:22
fangsteel
木虫 (正式写手)
- 应助: 96 (初中生)
- 金币: 1213.6
- 散金: 32
- 红花: 22
- 帖子: 456
- 在线: 192.1小时
- 虫号: 2248690
- 注册: 2013-01-18
- 性别: GG
- 专业: 原子和分子物理
7楼2014-03-07 19:52:02
juanjuan0618
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 179.6
- 帖子: 94
- 在线: 23.6小时
- 虫号: 1160553
- 注册: 2010-12-02
- 专业: 实验地球化学和计算地球化
8楼2014-03-07 21:18:30
fangsteel
木虫 (正式写手)
- 应助: 96 (初中生)
- 金币: 1213.6
- 散金: 32
- 红花: 22
- 帖子: 456
- 在线: 192.1小时
- 虫号: 2248690
- 注册: 2013-01-18
- 性别: GG
- 专业: 原子和分子物理
9楼2014-03-08 13:48:51
20