<html>
<head>
<style>
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
FONT-SIZE: 10pt;
FONT-FAMILY:Tahoma
}
</style>
</head>
<body class='hmmessage'><div style="text-align: left;">Hi,<br><br>There can be many reasons for such problems.<br>I assume you use solvent.<br><br>One option is that you could have made your whole protein<br>a single charge group. grompp in Gromacs 3.3.3 checks for this,<br>but older versions do not.<br></div><br>Berk.<br><br><br><hr id="stopSpelling">&gt; From: affinito@gmail.com<br>&gt; To: gmx-users@gromacs.org<br>&gt; Date: Thu, 8 May 2008 14:54:27 +0200<br>&gt; Subject: [gmx-users] crash nsgrid.c<br>&gt; <br>&gt; <br>&gt; Dear all,<br>&gt; During my simulations of a GFP protein (force field made with amber,  <br>&gt; then converted to gromacs), I usually get this error:<br>&gt; <br>&gt; Fatal error: ci = 8108 should be in 0 .. 7937 [FILE nsgrid.c, LINE 218]<br>&gt; <br>&gt; What this that means?<br>&gt; <br>&gt; If I restart the simulation from the frame just before the crash, the  <br>&gt; run is ok for a (random) time and then it crashes again.<br>&gt; I checked this on different kind of machines and I get always the same  <br>&gt; behavior.<br>&gt; <br>&gt; To make things clearier, I include the input.mdp file:<br>&gt; <br>&gt; title                    =<br>&gt; cpp                      = /lib/cpp<br>&gt; include                  =<br>&gt; define                   =<br>&gt; integrator               = md<br>&gt; dt                       = 0.002<br>&gt; nsteps                   = 500000<br>&gt; nstxout                  = 10000<br>&gt; nstvout                  = 0<br>&gt; nstlog                   = 500<br>&gt; nstenergy                = 500<br>&gt; nstxtcout                = 500<br>&gt; energygrps               = System<br>&gt; pbc                      = xyz<br>&gt; nstlist                  = 15<br>&gt; epsilon_r                = 1.<br>&gt; ns_type                  = grid<br>&gt; coulombtype              = pme<br>&gt; vdwtype                  = Cut-off<br>&gt; rlist                    = 0.8<br>&gt; rcoulomb                 = 0.8<br>&gt; rvdw                     = 0.8<br>&gt; tcoupl                   = Nose-Hoover<br>&gt; tc-grps                  = System<br>&gt; tau_t                    = 2.0<br>&gt; ref_t                    = 300.00<br>&gt; pcoupl                   = Parrinello-Rahman<br>&gt; pcoupltype               = isotropic<br>&gt; tau_p                    = 4.0<br>&gt; compressibility          = 4.5e-5<br>&gt; ref_p                    = 1.0<br>&gt; gen_vel                  = no<br>&gt; gen_seed                 = 173529<br>&gt; constraints              = All-bonds<br>&gt; constraint_algorithm     = lincs<br>&gt; shake_tol                = 0.0001<br>&gt; <br>&gt; Any suggestion is welcome.<br>&gt; <br>&gt; Thanks<br>&gt; <br>&gt; F.A.<br>&gt; _______________________________________________<br>&gt; gmx-users mailing list    gmx-users@gromacs.org<br>&gt; http://www.gromacs.org/mailman/listinfo/gmx-users<br>&gt; Please search the archive at http://www.gromacs.org/search before posting!<br>&gt; Please don't post (un)subscribe requests to the list. Use the <br>&gt; www interface or send it to gmx-users-request@gromacs.org.<br>&gt; Can't post? Read http://www.gromacs.org/mailing_lists/users.php<br><br /><hr />Express yourself instantly with MSN Messenger! <a href='http://clk.atdmt.com/AVE/go/onm00200471ave/direct/01/' target='_new'>MSN Messenger</a></body>
</html>