24小时热门版块排行榜    

查看: 2722  |  回复: 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
这个怎么办?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangsteel

木虫 (正式写手)

引用回帖:
6楼: Originally posted by juanjuan0618 at 2014-03-07 19:44:22
你好,请问下你,这个NOTE2是怎么解决的呢?  
我现在也遇到这个问题,多谢啦

看看你的top是怎么写的,有没有把几个原子写入一个char group
gromacs
7楼2014-03-07 19:52:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fangsteel

木虫 (正式写手)


月只蓝: 金币+1, 鼓励交流! 2014-03-09 10:10:59
引用回帖:
8楼: Originally posted by juanjuan0618 at 2014-03-07 21:18:30
我写的丙烷的全原子分子,总共有11个原子,我把它写成一个charge group和3个charge group  都出现这个问题...

试试把每个原子写成一个char group
gromacs
9楼2014-03-08 13:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 蓝紫色玻璃心 的主题更新
信息提示
请填处理意见