| 查看: 2236 | 回复: 3 | |||
[求助]
Gromacs水中的溶菌酶npt平衡出错
|
|
RT，nvt平衡可以正常运行，但是到npt平衡时，会出错，grompp终止运行 下面是命令及错误信息： [wanghj@c0103 donglixue]$ grompp_mpi -f npt.mdp -c em.gro -t nvt.cpt -p topol.top -o npt.tpr :-) G R O M A C S (-: Grunge ROck MAChoS :-) VERSION 4.5.5 (-: Written by Emile Apol, Rossen Apostolov, Herman J.C. Berendsen, Aldert van Buuren, P盲r Bjelkmar, Rudi van Drunen, Anton Feenstra, Gerrit Groenhof, Peter Kasson, Per Larsson, Pieter Meulenhoff, Teemu Murtola, Szilard Pall, Sander Pronk, Roland Schulz, Michael Shirts, Alfons Sijbers, Peter Tieleman, Berk Hess, David van der Spoel, and Erik Lindahl. Copyright (c) 1991-2000, University of Groningen, The Netherlands. Copyright (c) 2001-2010, The GROMACS development team at Uppsala University & The Royal Institute of Technology, Sweden. check out http://www.gromacs.org for more information. 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_mpi (-: Option Filename Type Description ------------------------------------------------------------ -f npt.mdp Input grompp input file with MD parameters -po mdout.mdp Output grompp input file with MD parameters -c em.gro Input Structure file: gro g96 pdb tpr etc. -r conf.gro Input, Opt. Structure file: gro g96 pdb tpr etc. -rb conf.gro Input, Opt. Structure file: gro g96 pdb tpr etc. -n index.ndx Input, Opt. Index file -p topol.top Input Topology file -pp processed.top Output, Opt. Topology file -o npt.tpr Output Run input file: tpr tpb tpa -t nvt.cpt Input, Opt! Full precision trajectory: trr trj cpt -e ener.edr Input, Opt. Energy file Option Type Value Description ------------------------------------------------------ -[no]h bool no Print help info and quit -[no]version bool no Print version info and quit -nice int 0 Set the nicelevel -[no]v bool no Be loud and noisy -time real -1 Take frame at or first after this time. -[no]rmvsbds bool yes Remove constant bonded interactions with virtual sites -maxwarn int 0 Number of allowed warnings during input processing. Not for normal use and may generate unstable systems -[no]zero bool no Set parameters for bonded interactions without defaults to zero instead of generating an error -[no]renum bool yes Renumber atomtypes and minimize number of atomtypes Ignoring obsolete mdp entry 'title' Back Off! I just backed up mdout.mdp to ./#mdout.mdp.11# Generated 332520 of the 332520 non-bonded parameter combinations Generating 1-4 interactions: fudge = 0.5 Generated 332520 of the 332520 1-4 parameter combinations Excluding 3 bonded neighbours molecule type 'Protein_chain_A' turning all bonds into constraints... Excluding 2 bonded neighbours molecule type 'SOL' turning all bonds into constraints... Excluding 1 bonded neighbours molecule type 'CL' turning all bonds into constraints... Analysing residue names: There are: 129 Protein residues There are: 10824 Water residues There are: 8 Ion 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 Protein is 3895.83 Number of degrees of freedom in T-Coupling group non-Protein is 64965.17 WARNING 1 [file npt.mdp]: You are using pressure coupling with absolute position restraints, this will give artifacts. Use the refcoord_scaling option. Reading Coordinates, Velocities and Box size from old trajectory Will read whole trajectory Last frame -1 time 6.300 Using frame at t = 6.3 ps Starting time for run is 0 ps Largest charge group radii for Van der Waals: 0.163, 0.160 nm Largest charge group radii for Coulomb: 0.172, 0.165 nm Calculating fourier grid dimensions for X Y Z Using a fourier grid of 45x45x45, spacing 0.156 0.156 0.156 Estimate for the relative computational load of the PME mesh part: 0.25 This run will generate roughly 399 Mb of data There was 1 warning ------------------------------------------------------- Program grompp_mpi, VERSION 4.5.5 Source code file: grompp.c, line: 1584 Fatal error: Too many warnings (1), grompp_mpi terminated. If you are sure all warnings are harmless, use the -maxwarn option. For more information and tips for troubleshooting, please check the GROMACS website at http://www.gromacs.org/Documentation/Errors ------------------------------------------------------- "I Caught It In the Face" (P.J. Harvey) 但是在在命令后面加-maxwarn 3 时，即 grompp_mpi -f npt.mdp -c em.gro -t nvt.cpt -p topol.top -o npt.tpr -maxwarn 3就可以正常了，请问这个错误该怎么解决呢？ |
回复此楼» 猜你喜欢
071000生物学求调剂，初试成绩343
已经有7人回复
-
0703化学调剂，求导师收
已经有8人回复
-
0856材料化工调剂 总分330
已经有11人回复
-
本科新能源科学与工程，一志愿华理能动285求调剂
已经有9人回复
-
332求调剂
已经有4人回复
-
085600，专业课化工原理，321分求调剂
已经有3人回复
-
材料工程294求调剂
已经有3人回复
-
295材料工程专硕求调剂
已经有7人回复
-
305求调剂
已经有5人回复
-
食品工程专硕一志愿中海洋309求调剂
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 磷酸亚铁锂的平衡电势曲线
已经有6人回复
- gromacs模拟
已经有4人回复
- 化学计量平衡与热力学平衡的关系
已经有17人回复
- 求助碳纳米管对BSA的吸附平衡时间
已经有4人回复
- 求助：推荐相平衡专家
已经有3人回复
- 请问溶菌酶和溶壁酶是一种酶么？谢谢各位大侠啦
已经有6人回复
- 液相预热、柱平衡问题
已经有7人回复
- 在gromacs中运行mdrun命令时出错
已经有10人回复
- 求Gromacs中Glycam力场数据包。。非常感谢
已经有5人回复
- 【求助】如何使用gromacs中的constraint
已经有8人回复
- 【求助】这是我安装gromacs过程，但是出错阿，这个怎么办呢？谢谢
已经有11人回复
- 【求助】Gromacs和NAMD比Amber快很多吗？如何选择？
已经有7人回复
- 【求助】NAMD和Gromacs哪个简单？
已经有5人回复
2楼2012-04-26 10:00:34
3楼2012-04-26 10:54:29
andylau286
铁虫 (小有名气)
- 应助: 3 (幼儿园)
- 金币: 696.1
- 红花: 1
- 帖子: 125
- 在线: 64.8小时
- 虫号: 1621704
- 注册: 2012-02-16
- 性别: GG
- 专业: 生物信息学
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
zh1987hs: 金币+2, 谢谢 2012-04-26 21:15:03
604758762: 金币+3, ★有帮助 2012-05-03 21:11:50
感谢参与，应助指数 +1
zh1987hs: 金币+2, 谢谢 2012-04-26 21:15:03
604758762: 金币+3, ★有帮助 2012-05-03 21:11:50
|
WARNING 1 [file npt.mdp]: You are using pressure coupling with absolute position restraints, this will give artifacts. Use the refcoord_scaling option 这个不算错误，只是系统提示你，你正在用压力耦合，这会给系统带来人工成分，可能会影响到严谨性。 但这个错误是无害的，可以在命令的最后加 “ -maxwarn 1 "就可以运行了， |
4楼2012-04-26 21:11:19
