<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7233.28">
<TITLE>TFE-water minimization</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>

<P><FONT SIZE=2>Hello,<BR>
<BR>
I have created my own TFE.itp file for a mixed solvent simulation containing TFE and TIP4P water.&nbsp; To test my TFE file I put one TFE molecule in a box of solvent and and tried to minimize it using the following parameters:<BR>
<BR>
<BR>
define&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; =&nbsp; -DPOSRES<BR>
cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; /lib/cpp<BR>
;constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; all-bonds<BR>
integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; steep<BR>
nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 2000<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; =&nbsp; 0.05<BR>
emstep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.01<BR>
nstcgsteep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1000<BR>
nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1<BR>
ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; grid<BR>
;deltagrid&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 2<BR>
rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.4<BR>
rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.4<BR>
rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1.4<BR>
Tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<BR>
Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<BR>
gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; no<BR>
lincs_iter&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; =&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>
<BR>
<BR>
Unfortunately my run ends in error with the following message:<BR>
<BR>
------------------------------------------------------<BR>
Program mdrun, VERSION 3.3.1<BR>
Source code file: nsgrid.c, line: 226<BR>
<BR>
Range checking error:<BR>
Explanation: During neighborsearching, we assign each particle to a grid<BR>
based on its coordinates. If your system contains collisions or parameter<BR>
errors that give particles very high velocities you might end up with some<BR>
coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot<BR>
put these on a grid, so this is usually where we detect those errors.<BR>
Make sure your system is properly energy-minimized and that the potential<BR>
energy seems reasonable before trying again.<BR>
<BR>
Variable ci has value -2147483648. It should have been within [ 0 .. 512 ]<BR>
Please report this to the mailing list (gmx-users@gromacs.org)<BR>
-------------------------------------------------------<BR>
<BR>
<BR>
Can you tell me what I can do to solve this problem?&nbsp;<BR>
<BR>
Thank youin advance,<BR>
<BR>
Michael Owen</FONT>
</P>

</BODY>
</HTML>