<u>We are trying to do simulation of lysozyme in water.</u><br><br>step with problem : energy minimization.<br><br>command given: <b style="color:rgb(255,0,0)">[root@localhost gromacs-4.0.5]# grompp -f minim.mdp -c 1AKI_solv_ions.gro -p topol.top -o em.tpr</b><br>
<br> :-) G R O M A C S (-:<br><br> Getting the Right Output Means no Artefacts in Calculating Stuff<br><br> :-) VERSION 4.0.5 (-:<br><br><br> Written by David van der Spoel, Erik Lindahl, Berk Hess, and others.<br>
Copyright (c) 1991-2000, University of Groningen, The Netherlands.<br> Copyright (c) 2001-2008, The GROMACS development team,<br> check out <a href="http://www.gromacs.org">http://www.gromacs.org</a> for more information.<br>
<br> This program is free software; you can redistribute it and/or<br> modify it under the terms of the GNU General Public License<br> as published by the Free Software Foundation; either version 2<br>
of the License, or (at your option) any later version.<br><br> :-) grompp (-:<br><br>Option Filename Type Description<br>------------------------------------------------------------<br>
-f minim.mdp Input, Opt! grompp input file with MD parameters<br> -po mdout.mdp Output grompp input file with MD parameters<br> -c 1AKI_solv_ions.gro Input Structure file: gro g96 pdb tpr tpb tpa<br>
-r conf.gro Input, Opt. Structure file: gro g96 pdb tpr tpb tpa<br> -rb conf.gro Input, Opt. Structure file: gro g96 pdb tpr tpb tpa<br> -n index.ndx Input, Opt. Index file<br> -p topol.top Input Topology file<br>
-pp processed.top Output, Opt. Topology file<br> -o em.tpr Output Run input file: tpr tpb tpa<br> -t traj.trr Input, Opt. Full precision trajectory: trr trj cpt<br> -e ener.edr Input, Opt. Energy file: edr ene<br>
<br>Option Type Value Description<br>------------------------------------------------------<br>-[no]h bool no Print help info and quit<br>-nice int 0 Set the nicelevel<br>-[no]v bool yes Be loud and noisy<br>
-time real -1 Take frame at or first after this time.<br>-[no]rmvsbds bool yes Remove constant bonded interactions with virtual<br> sites<br>-maxwarn int 0 Number of allowed warnings during input processing<br>
-[no]zero bool no Set parameters for bonded interactions without<br> defaults to zero instead of generating an error<br>-[no]renum bool yes Renumber atomtypes and minimize number of<br>
atomtypes<br><br><br>Back Off! I just backed up mdout.mdp to ./#mdout.mdp.1#<br>checking input for internal consistency...<br>processing topology...<br>Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaa.itp<br>
Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaanb.itp<br>Opening library file /usr/local/gromacs/share/gromacs/top/ffoplsaabon.itp<br>Generated 332520 of the 332520 non-bonded parameter combinations<br>
Generating 1-4 interactions: fudge = 0.5<br>Generated 332520 of the 332520 1-4 parameter combinations<br>Opening library file /usr/local/gromacs/share/gromacs/top/spce.itp<br>Opening library file /usr/local/gromacs/share/gromacs/top/ions.itp<br>
Excluding 3 bonded neighbours molecule type 'Protein_A'<br>Excluding 2 bonded neighbours molecule type 'SOL'<br><br>-------------------------------------------------------<br>Program grompp, VERSION 4.0.5<br>
Source code file: toppush.c, line: 1641<br><br><b style="color:rgb(255,0,0)">Fatal error:<br>No such moleculetype NA</b><br><br>what should we do????<br><br>