| 查看: 4588 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Gromacs中加入离子
已经有5人回复
- [新手求教]安装Gromacs后，查看pdb2gmx命令，显示无此命令
已经有9人回复
- 【讨论】gromacs电荷不为整数
已经有6人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【原创】在gromacs的gro中插入小分子
已经有12人回复
- 【求助】gromacs出现了一个问题
已经有10人回复
- 【求助】求一个gromacs的中文教程，包括安装的部分，小弟新手，从零学起
已经有8人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
5楼2010-12-22 16:58:39
2楼2010-12-22 15:18:04
3楼2010-12-22 16:10:49
4楼2010-12-22 16:54:22
简单回复
matao01114楼
2011-01-01 01:09
回复
zyj8119(金币+1):谢谢参与
20