<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>Did you check how your box dimensions evolve?<br><br>For a homogeneous system (as your small test system)<br>one should not use semi or anisotropic pressure coupling,<br>since the box shape can distort unlimitedly causing<br>the system to crash.<br><br>Berk.<br><br></div><br><hr id="stopSpelling">&gt; Date: Tue, 3 Jun 2008 17:00:14 +0200<br>&gt; From: till@bph.rub.de<br>&gt; To: gmx-users@gromacs.org<br>&gt; Subject: [gmx-users] neighborsearching in semiisotropic systems<br>&gt; <br>&gt; Dear users,<br>&gt; <br>&gt; During protein simulations occurred the range checking error of the <br>&gt; neighborsearching. I checked the comments on this mistake in the user <br>&gt; list and then checked my system on any hints of an exploding system but <br>&gt; it seems to be stable. In order to be sure I cut down the system to a <br>&gt; small peptide consisting of only ten aminoacids and used the ffG43a1 <br>&gt; forcefield.<br>&gt; <br>&gt; <br>&gt; I performed 10ns simulations of this small system in a cubic <br>&gt; spc-waterbox with different pressure parameters and tried single and <br>&gt; double precision. I got the following results:<br>&gt; <br>&gt; <br>&gt; For<br>&gt; <br>&gt; pcoupltype = semiisotropic<br>&gt; <br>&gt; tau_p = 0.5 0.5<br>&gt; <br>&gt; compressibility = 4.5E-5 4.5E-5<br>&gt; <br>&gt; ref_p = 1 1<br>&gt; <br>&gt; <br>&gt; with double precision occurred the mistake:<br>&gt; <br>&gt; <br>&gt; -------------------------------------------------------<br>&gt; <br>&gt; Program mdrun_d.mpi, VERSION 3.3.1<br>&gt; <br>&gt; Source code file: nsgrid.c, line: 226<br>&gt; <br>&gt; <br>&gt; Range checking error:<br>&gt; <br>&gt; Explanation: During neighborsearching, we assign each particle to a grid<br>&gt; <br>&gt; based on its coordinates. If your system contains collisions or parameter<br>&gt; <br>&gt; errors that give particles very high velocities you might end up with some<br>&gt; <br>&gt; coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot<br>&gt; <br>&gt; put these on a grid, so this is usually where we detect those errors.<br>&gt; <br>&gt; Make sure your system is properly energy-minimized and that the potential<br>&gt; <br>&gt; energy seems reasonable before trying again.<br>&gt; <br>&gt; <br>&gt; Variable ci has value 284. It should have been within [ 0 .. 200 ]<br>&gt; <br>&gt; Please report this to the mailing list (gmx-users@gromacs.org <br>&gt; &lt;mailto:gmx-users@gromacs.org&gt;)<br>&gt; <br>&gt; <br>&gt; Then, I changed tau_p to 5.0 and the mistake did not occurce anymore <br>&gt; during the 10 ns. With single precision the mistake never occurres. I <br>&gt; also tested double precision with pcoupltype = isotropic and this also <br>&gt; works.<br>&gt; <br>&gt; Next, I wanted to reproduce the mistake so I took the same tpr file of <br>&gt; the crashed simulation with tau_p = 0.5 and started the simulation <br>&gt; again. The mistake occurred again but not at the same timestep. The <br>&gt; first run crashed after step 1558630 and the second run after step <br>&gt; 2597500 and also the ci variable is not exactly the same (Variable ci <br>&gt; has value 283. It should have been within [ 0 .. 200 ]).<br>&gt; <br>&gt; <br>&gt; I think there might be a problem between the neighborsearching and the <br>&gt; semiisotropic pressurcoupling together with double precision. I am quite <br>&gt; sure that the peptide is not the problem. Maybe there is a problem with <br>&gt; a water molecule. I would be glad if someone can give me an explanation <br>&gt; why this mistake is not reproducable at the same timestep and why it <br>&gt; occures with tau_p=0.5 and not with tau_p=5.0.<br>&gt; <br>&gt; <br>&gt; Thanks<br>&gt; <br>&gt; <br>&gt; Till<br>&gt; <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>