| 查看: 4595 | 回复: 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 |
回复此楼» 猜你喜欢
国自然面上复盘~欢迎讨论
已经有9人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
基金未中，这种答复是模板吗？
已经有5人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有5人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gromacs中加入离子
已经有5人回复
- [新手求教]安装Gromacs后，查看pdb2gmx命令，显示无此命令
已经有9人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】gromacs出现了一个问题
已经有10人回复
- 【求助】求一个gromacs的中文教程，包括安装的部分，小弟新手，从零学起
已经有8人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
发点金币玩
+5/715
-
面上得中
+1/493
-
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/253
-
中国基，散金币！！！！！！
+5/204
-
东华大学材料科学与工程学院丁奇老师招收2026级硕士研究生
+1/100
-
北京化工大学何静教授团队诚聘青年英才和博士后
+1/100
-
中国科学院苏州纳米技术与纳米仿生研究所 仿生自适应智能系统实验室招收联培生
+1/42
-
现在生物医药相关有什么兼职可以做吗？
+1/36
-
新审查指南堵死了"一案双申"的老路，实用新型和发明只能二选一了？
+1/23
-
国家蛋白质科学中心（北京），招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/22
-
国家蛋白质科学中心（北京），招收分析化学、药学和蛋白质组代谢组学博士研究生
+1/20
-
浙江大学医学院黄浩杰教授团队招聘4-5名博士后
+1/9
-
细胞治疗全流程关键检测技术与研究策略
+1/5
-
项目助理招聘启事
+1/5
-
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/4
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/2
-
量子物理团队招收博士生与博士后
+1/1
-
重大突破！酶法合成破解生物法DNA合成原料“卡脖子”难题
+1/1
-
【2027年9月入学】宁波东方理工大学联培博士生(中科大/上交大)招生（物理/材料方向）
+1/1
2楼2010-12-22 15:18:04
3楼2010-12-22 16:10:49
4楼2010-12-22 16:54:22
5楼2010-12-22 16:58:39
|
static void check_gbsa_params(t_inputrec *ir,gpp_atomtype_t atype) { int nmiss,i; /* If we are doing GBSA, check that we got the parameters we need * This checking is to see if there are GBSA paratmeters for all * atoms in the force field. To go around this for testing purposes * comment out the nerror++ counter temporarily */ nmiss = 0; for(i=0;i if (get_atomtype_radius(i,atype) < 0 || get_atomtype_vol(i,atype) < 0 || get_atomtype_surftens(i,atype) < 0 || get_atomtype_gb_radius(i,atype) < 0 || get_atomtype_S_hct(i,atype) < 0) { fprintf(stderr,"\nGB parameter(s) missing or negative for atom type '%s'\n", get_atomtype_name(i,atype)); nmiss++; } } |
6楼2010-12-22 17:10:30
★
zyj8119(金币+1):谢谢参与
zyj8119(金币+1):谢谢参与
7楼2010-12-22 21:23:43
8楼2010-12-22 23:32:22
9楼2010-12-23 10:28:51
10楼2010-12-25 11:42:26
11楼2011-01-01 00:08:59
12楼2011-01-01 00:36:42
13楼2011-01-01 00:43:38
15楼2011-01-01 14:21:47
17楼2011-01-06 10:36:12
18楼2011-01-06 17:23:05
19楼2011-01-07 17:24:55
20楼2011-01-07 17:25:09
21楼2011-03-09 13:19:10
22楼2011-03-09 15:40:58
23楼2011-03-09 18:34:30
简单回复
matao01114楼
2011-01-01 01:09
回复
zyj8119(金币+1):谢谢参与
iceem16楼
2011-01-01 23:59
回复
zyj8119(金币+1):谢谢参与
