24小时热门版块排行榜    

查看: 4614  |  回复: 22
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[交流] 【求助】gromacs的又一个问题

我运行 gromcs,现在是没有原子数不等的问题了,但是又出现了这个问题:
    This program is free software; you can redistribute it and/or
          modify it under the terms of the GNU General Public License
         as published by the Free Software Foundation; either version 2
             of the License, or (at your option) any later version.

                                :-)  grompp  (-:

Option     Filename  Type         Description
------------------------------------------------------------
  -f      minim.mdp  Input, Opt!  grompp input file with MD parameters
-po      mdout.mdp  Output       grompp input file with MD parameters
  -c      total.gro  Input        Generic structure: gro g96 pdb tpr tpb tpa
                                   xml
  -r       conf.gro  Input, Opt.  Generic structure: gro g96 pdb tpr tpb tpa
                                   xml
-rb       conf.gro  Input, Opt.  Generic structure: gro g96 pdb tpr tpb tpa
                                   xml
  -n      index.ndx  Input, Opt.  Index file
-deshuf  deshuf.ndx  Output, Opt. Index file
  -p        sys.top  Input        Topology file
-pp  processed.top  Output, Opt. Topology file
  -o      input.tpr  Output       Generic run input: tpr tpb tpa xml
  -t       traj.trr  Input, Opt.  Full precision trajectory: trr trj
  -e       ener.edr  Input, Opt.  Generic energy: edr ene

      Option   Type  Value  Description
------------------------------------------------------
      -[no]h   bool     no  Print help info and quit
       -nice    int      0  Set the nicelevel
      -[no]v   bool    yes  Be loud and noisy
       -time   real     -1  Take frame at or first after this time.
         -np    int      1  Generate statusfile for # nodes
-[no]shuffle   bool     no  Shuffle molecules over nodes
   -[no]sort   bool     no  Sort molecules according to X coordinate
-[no]rmvsbds   bool    yes  Remove constant bonded interactions with virtual
                            sites
       -load string         Releative load capacity of each node on a
                            parallel machine. Be sure to use quotes around
                            the string, which should contain a number for
                            each node
    -maxwarn    int     10  Number of warnings after which input processing
                            stops
-[no]check14   bool     no  Remove 1-4 interactions without Van der Waals
  -[no]renum   bool    yes  Renumber atomtypes and minimize number of
                            atomtypes

creating statusfile for 1 node...
' for variable optimize_fft, using 'no'
Next time use one of: 'no' 'yes'

Back Off! I just backed up mdout.mdp to ./#mdout.mdp.3#
checking input for internal consistency...
calling /lib/cpp...
processing topology...
Generated 165 of the 1596 non-bonded parameter combinations
Excluding 3 bonded neighbours for drg 64
turning H bonds into constraints...
Excluding 2 bonded neighbours for SOL 218
turning H bonds into constraints...
processing coordinates...
double-checking input for internal consistency...
WARNING 1 [file "sys.top", line 37]:
  For minimization with LINCS constraints, lincs_iter should be 4 to 8.
There was 1 warning
-------------------------------------------------------
Program grompp, VERSION 3.3.1
Source code file: grompp.c, line: 1109

Fatal error:
There were 1 error(s) processing your input
-------------------------------------------------------

"Jesus Not Only Saves, He Also Frequently Makes Backups." (Myron Bradshaw)
我的sys.top是:
;
;       File 'sys.top' was generated
;       By user: zhu (503)
;       On host: apollo2
;       At date: Wed Mar 12 15:48:05 2008
;
;       This is your topology file
;       DRG
;
; Include forcefield parameters
#include "ffG53a6.itp"

; Include surfacetant topology
#include "drg.itp"

; Include water topology
#include "spce.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct       fcx        fcy        fcz
   1    1       1000       1000       1000
#endif

; Include generic topology for ions
#include "ions.itp"

[ system ]
; Name
aok in water

[ molecules ]
; Compound        #mols
drg                 64

SOL               218
,而minim.mdp文件为:
; LINES STARTING WITH ';' ARE COMMENTS
title           = Minimization of AOK    ; Title of run

; The following lines tell the program the standard locations where to find cert
ain files
cpp             = /lib/cpp      ; Preprocessor
include         = -I../top      ; Directories to include in the topology format

; Parameters describing what to do, when to stop and what to save
integrator      = steep         ; Algorithm (steep = steepest descent minimizati
on)
emtol           = 1000.0  ; Stop minimization when the maximum force < 3000.0 kJ
/mol
nsteps          = 100000        ; Maximum number of (minimization) steps to perf
orm
nstenergy       = 10            ; Write energies to disk every nstenergy steps
nstxtcout       = 10            ; Write coordinates to disk every nstxtcout step
s
xtc_grps        = system               ; Which coordinate group(s) to write to d
isk
energygrps      = system               ; Which energy group(s) to write to disk

; Parameters describing how to find the neighbors of each atom and how to calcul
ate the interactions
nstlist         = 5             ; Frequency to update the neighbor list and long
range forces
ns_type         = simple        ; Method to determine neighbor list (simple, gri
d)
rlist             = 1.0 ; Cut-off for making neighbor list (short range forces)
coulombtype              = PME  ; Treatment of long range electrostatic interact
ions
rcoulomb                 = 1.0
fourierspacing       =0.12
pme_order            =4
ewald_rtol           =1e-5
optimize_fft         =yes
rvdw            = 1.0           ; long range Van der Waals cut-off
constraints     = hbonds        ; Bond types to replace by constraints
pbc             = xyz           ; Periodic Boundary Conditions (yes/no)

;freeze
;freezegrps               = SOL
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zyj8119 at 2011-01-07 17:24:55:

此问题解决,我用了最新的gromacs,就什么问题也没有。。。是gmx的版本问题。。

请版主返还金币。。。
20楼2011-01-07 17:25:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答
引用回帖:
Originally posted by zyj8119 at 2010-12-22 14:20:15:
我运行 gromcs,现在是没有原子数不等的问题了,但是又出现了这个问题:
    This program is free software; you can redistribute it and/or
          modify it under the terms of the GNU General Public  ...

重金酬谢。。。
2楼2010-12-22 15:18:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shengxiang

木虫 (著名写手)



zyj8119(金币+1):谢谢参与
引用回帖:
Originally posted by zyj8119 at 2010-12-22 14:20:15:
我运行 gromcs,现在是没有原子数不等的问题了,但是又出现了这个问题:
    This program is free software; you can redistribute it and/or
          modify it under the terms of the GNU General Public  ...

用过一段时间的GROMACS,但是你提到的这个问题还是第一次碰到……无能为了了~
3楼2010-12-22 16:10:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mophyworld

金虫 (正式写手)


★ ★ ★
zyj8119(金币+1):谢谢参与
ghcacj(金币+2):谢谢 2010-12-23 15:52:35
Source code file: grompp.c, line: 1109
你可以查看下这行的代码,找到源头
4楼2010-12-22 16:54:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
matao01114楼
2011-01-01 01:09   回复  
zyj8119(金币+1):谢谢参与
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +13 Kittylucky 2026-08-27 14/700 2026-09-01 11:06 by feng6531
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
信息提示
请填处理意见