Hi,<br> I encountered the following remarks during normal mode analysis of pure water. Can someone please suggest a solution.<br> mdrun_d -v -s water.tpr -o traj.trr -c confout.gro -e ener.edr -g md.log -mtx nm.mtx<br>Getting Loaded...<br>
Reading file water.tpr, VERSION 4.0.4 (double precision)<br>Loaded with Money<br><br>Non-cutoff electrostatics used, forcing full Hessian format.<br>Allocating Hessian memory...<br><br>starting normal mode calculation 'Pure Water'<br>
648 steps.<br><br>Maximum force: 2.74758e+03<br>Maximum force probably not small enough to ensure that you are in an<br>energy well. Be aware that negative eigenvalues may occur when the<br>resulting matrix is diagonalized.<br>
Finished step 648 out of 648<br><br>Writing Hessian...<br><br>gcq#258: "In the End Science Comes Down to Praying" (P. v.d. Berg)<br><br>g_nmeig_d -f nm.mtx -s water.tpr -of eigenfreq.xvg -ol eigenval.xvg -v eigenvec.trr<br>
<br>Reading file water.tpr, VERSION 4.0.4 (double precision)<br>Reading file water.tpr, VERSION 4.0.4 (double precision)<br>Reading double precision matrix generated by Gromacs VERSION 4.0.4<br>Full matrix storage format, nrow=1944, ncols=1944<br>
<br>Diagonalizing to find vectors 1 through 50...<br><br>One of the lowest 6 eigenvalues has a non-zero value.<br>This could mean that the reference structure was not<br>properly energy minimized.<br>Writing eigenvalues...<br>
Writing eigenfrequencies - negative eigenvalues will be set to zero.<br><br>Writing average structure & eigenvectors 1--50 to eigenvec.trr<br><br>gcq#132: "Shit Happens" (Pulp Fiction)<br><br><br>The eigenfreq is 0 which is the problem<br>
<br>The .mdp file used was;<br>; User spoel (236)<br>; Wed Nov 3 17:12:44 1993<br>; Input file<br>;<br>cpp =<br>define = -DFLEXIBLE<br>constraints = none<br>integrator = nm<br>
nsteps = 10000<br>;<br>; Energy minimizing stuff<br>;<br>emtol = 0.01<br>emstep = 0.01<br>coulombtype = PME<br>nstcomm = 1<br>ns_type = grid<br>
rlist = 0.9<br>rcoulomb = 0.9<br>rvdw = 0.9<br>Tcoupl = no<br>Pcoupl = no<br>gen_vel = no<br><br><br><br><br>