| 查看: 4579 | 回复: 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 |
回复此楼» 猜你喜欢
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
能否退出参与的面上项目解除限项
已经有17人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
-
放榜前的不淡定
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gromacs中加入离子
已经有5人回复
- [新手求教]安装Gromacs后，查看pdb2gmx命令，显示无此命令
已经有9人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】gromacs出现了一个问题
已经有10人回复
- 【求助】求一个gromacs的中文教程，包括安装的部分，小弟新手，从零学起
已经有8人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
告别流水线爽文！108万字完结仙侠《女娲石之灵石奇缘（仙界篇）》值得细品
+1/5787
-
小木虫认识的老公在一起11年，去世快100天了。
+2/1986
-
2026年工程管理与智慧城市国际会议 (EMSC 2026)
+1/949
-
基金申请小调查！祝申请必定中！！祈福！！！
+5/395
-
基金祈福
+1/161
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
+2/86
-
武汉，诚征男友，不异地，年龄大了
+1/39
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/22
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/21
-
这些国际期刊正在招募青年编委
+1/15
-
安徽医科大学周可成教授课题组诚聘博士后 2026年
+1/14
-
安徽医科大学周可成教授课题组诚聘博士后 2026年
+1/11
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/4
-
中科院大连化学物理研究所 招聘 催化剂研发方向 科研助理 2名
+1/4
-
南京大学 智能驱动与感知材料实验室 诚招申请考核博士生/科研助理/博士后
+1/4
-
祈祷老公中国基
+1/3
-
上海大学王亮课题组诚招博士后（碳材料、光催化、电催化方向）
+2/2
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/2
-
新审查指南施行半年，AI方向的专利申请是不是更看"技术味"了？
+1/2
-
【产品测评分享】体外转录产量不理想、dsRNA 副产物高？T7 RNA 聚合酶实测体验
+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):谢谢参与
