<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman,new york,times,serif;font-size:12pt"><div>Hi!<br><br>I had the same problem and this is what worked for me:<br>1.&nbsp; You must do an equilibration run WITHOUT any constraints (replace constrained bonds with harmonic bonds (equilibrium bond length = constraint length and force constant suitably high) and increase weight of H atoms to 10.0.<br>2.&nbsp; Do an equilibration run ( 1ns if practical) with this system.<br>3.&nbsp; Use the confout.gro file it generates to start your run WITH constraints and the canonical H atom weights.<br><br>This should work.<br><br>Ashwin<br></div><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><br><div style="font-family: times new roman,new york,times,serif; font-size: 12pt;"><font size="2" face="Tahoma"><hr size="1"><b><span style="font-weight: bold;">From:</span></b>
 Ernesto Andres Román &lt;indio18@hotmail.com&gt;<br><b><span style="font-weight: bold;">To:</span></b> gmx-users@gromacs.org<br><b><span style="font-weight: bold;">Sent:</span></b> Thursday, 28 May, 2009 1:10:54 PM<br><b><span style="font-weight: bold;">Subject:</span></b> [gmx-users] LINCS warnings (6153)<br></font><br>

<style>
.hmmessage P
{
margin:0px;padding:0px;}
body.hmmessage
{
font-size:10pt;font-family:Verdana;}
</style>

Dear gmx-user,<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Im just making my first moves in protein dynamics simulation. Ive done some tutorials and peptide-water simulations. But now im intending to change the solvent to CCL4. Ive got the solvent parameters, made the topology file, and the .gro file with the solvent and peptide together. When I run mdrun commando an error message displays:<br>Fatal error:<br>Too many LINCS warnings (6153)<br>If you know what you are doing you can adjust the lincs warning threshold in your mdp file<br>or set the environment variable GMX_MAXCONSTRWARN to -1,<br>but normally it is better to fix the problem<br><br>My .mdp file is:<br>define&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = -DFLEXIBLE<br>constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
 none<br>integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = steep<br>dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.002&nbsp;&nbsp;&nbsp;&nbsp; ; ps !<br>nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 400<br>nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = grid<br>rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br>coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br>rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.0<br>vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = cut-off<br>rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =
 1.4<br>fourierspacing&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.12<br>fourier_nx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>fourier_ny&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>fourier_nz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>pme_order&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4<br>ewald_rtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1e-5<br>optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = yes<br>;<br>;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Energy minimizing stuff<br>;<br>emtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000.0<br>emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br><br>Do you know what could be happening, and more important, how to solve it or where to look for the problem ?<br>Thanks,<br>Ernesto.<br><br><br><hr>check out the rest of the Windows Live™.
More than mail–Windows Live™ goes way beyond your inbox.
 <a rel="nofollow" target="_blank" href="http://www.microsoft.com/windows/windowslive/">More than messages</a></div></div></div><br>


      <!--3--><hr size=1></hr> Explore and discover exciting holidays and getaways with Yahoo! India Travel <a href="http://in.rd.yahoo.com/tagline_Travel_1/*http://in.travel.yahoo.com/"> Click here!</a></body></html>