计算化学公社

 找回密码 Forget password
 注册 Register
Views: 4867|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 求助:构建聚合物残基文件运行报错

[复制链接 Copy URL]

67

帖子

0

威望

192

eV
积分
259

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
老师,您好,我构建了一个PVDF的.rtp文件,运行命令gmx pdb2gmx -f PVDF.pdb -o PVDF.gro出现这种错误:Problem with chain definition, or missing terminal residues.This chain does not appear to contain a recognized chain molecule.
If this is incorrect, you can edit residuetypes.dat to modify the behavior.
8 out of 8 lines of specbond.dat converted successfully
Checking for duplicate atoms....
Generating any missing hydrogen atoms and/or adding termini.
Now there are 3 residues with 21 atoms
Making bonds...
Number of bonds was 21, now 19
Generating angles, dihedrals and pairs...
-------------------------------------------------------
Program:     gmx pdb2gmx, version 2019.6
Source file: src\gromacs\gmxpreprocess\toputil.cpp (line 69)

Fatal error:
Increase MAXSLEN in the grompp code to at least 46, or shorten your definition
of bonds like 0.753   2.259   0.000  -3.012  -0.000   0.000 to at most 31
For more information and tips for troubleshooting, please check the GROMACS

website at http://www.gromacs.org/Documentation/Errors
我将PVDF的残基文件和结构文件附上,请老师帮我看一下是什么问题?这种问题该如何修改rtp呢?


PVDF.rtp

4.51 KB, 下载次数 Times of downloads: 7

PVDF.pdb

2.05 KB, 下载次数 Times of downloads: 3

125

帖子

0

威望

2930

eV
积分
3055

Level 5 (御坂)

2#
发表于 Post on 2022-9-23 15:00:41 | 只看该作者 Only view this author
0.753   2.259   0.000  -3.012  -0.000   0.000
这个意思就是说你给的数据太长,可以用一个define的字段来代替,比如:OS_C_CA_CA
然后使用gmx pdb2gmx的时候会在topol.top里面自动用OS_C_CA_CA字段代替你的参数,最后你只用在topol里面#define一下OS_C_CA_CA字段具体的参数值,就行了,注意键参数的函数类型需要删除和修改

具体可以参考TPPMKTOP生成RTP的文件,以及Gromacs自带的力场参数中的RTP文件和ffbonded.itp

本版积分规则 Credits rule

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

GMT+8, 2026-2-25 01:02 , Processed in 0.401727 second(s), 23 queries , Gzip On.

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