<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">

<html>

<head>

  <meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">

  <title></title>

</head>

<body bgcolor="#ffffff" text="#000000">

<pre><small><i> /&gt;/ 

</i>&gt;<i> /&gt;/ Range checking error:

</i>&gt;<i> /&gt;/ Explanation: During neighborsearching, we assign each particle to a grid

</i>&gt;<i> /&gt;/ based on its coordinates. If your system contains collisions or parameter

</i>&gt;<i> /&gt;/ errors that give particles very high velocities you might end up with some

</i>&gt;<i> /&gt;/ coordinates being +-Infinity or NaN (not-a-number). Obviously, we cannot

&gt; /&gt;/ put these on a grid, so this is usually where we detect those errors.

</i>&gt;<i> /&gt;/ Make sure your system is properly energy-minimized and that the potential

&gt; /&gt;/ energy seems reasonable before trying again.

&gt; /Please follow the advice in the above two sentences. Check your structure.

</i>&gt;<i> 

</i>&gt;<i> I did the energy minimization by increasing the number of steps using steepest descent, nothing happened. Then I tried cg using -DFLEXIBLE

</i>&gt;<i> and steep using -DFLEXIBLE in my mdp file and I was able to compute using mdrun but I got the following in my log file

</i>&gt;<i> 

</i>&gt;<i> Initiating Steepest Descents

</i>&gt;<i> Center of mass motion removal mode is Linear

</i>&gt;<i> We have the following groups for center of mass motion removal:

</i>&gt;<i>   0:  rest, initial mass: 24594

</i>&gt;<i> Started Steepest Descents on node 0 Mon Mar 13 17:07:04 2006

</i>&gt;<i> Removing pbc first time

</i>&gt;<i> Done rmpbc

</i>&gt;<i> Steepest Descents:

</i>&gt;<i>    Tolerance (Fmax)   =  1.00000e+01

</i>&gt;<i>    Number of steps    =        10000

</i>&gt;<i> Grid: 10 x 10 x 20 cells

&gt; Configuring nonbonded kernels...

&gt; Testing AMD 3DNow support... not present.

</i>&gt;<i> Testing ia32 SSE support... present.

</i>&gt;<i> 

</i>&gt;<i> 

</i>&gt;<i> Then with this when I do position restrained md, using -DPOSRES in mdp file for grompp and md integrator and then mdrun for position restrained

</i>&gt;<i> I get the following

</i>&gt;<i> 

</i>&gt;<i> Step 0, time 0 (ps)  LINCS WARNING

</i>&gt;<i> relative constraint deviation after LINCS:

</i>&gt;<i> max 1391229833052160.000000 (between atoms 655 and 657) rms 40665803128832.000000

</i>&gt;<i> bonds that rotated more than 30 degrees:

</i>

<big>turn off the constraints

</big>

&gt;</small></pre>

<br>

Hi David<br>

<br>

I did turn off the constraints and also made the starting structure

unconstarined in the mdp file. The mdrun for position restrained gives

Segmentation fault and the log file is shown below.<br>

Could you please suggest what could I do. I attached a pdb file in the

earlier mail which has to be reviewed by the moderator of gmx list.<br>

<br>

Thanks<br>

Gaurav<br>

<br>

<br>

-------------------------------------------<br>

Input Parameters:<br>

&nbsp;&nbsp; integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = md<br>

&nbsp;&nbsp; nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; init_step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Grid<br>

&nbsp;&nbsp; nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>

&nbsp;&nbsp; ndelta&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<br>

&nbsp;&nbsp; bDomDecomp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>

&nbsp;&nbsp; decomp_dir&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; comm_mode&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Linear<br>

&nbsp;&nbsp; nstcheckpoint&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10000<br>

&nbsp;&nbsp; nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstvout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 250<br>

&nbsp;&nbsp; nstfout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 250<br>

&nbsp;&nbsp; nstxtcout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 250<br>

&nbsp;&nbsp; init_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; delta_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.002<br>

&nbsp;&nbsp; xtcprec&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; nkx&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 42<br>

&nbsp;&nbsp; nky&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 42<br>

&nbsp;&nbsp; nkz&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 84<br>

&nbsp;&nbsp; pme_order&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4<br>

&nbsp;&nbsp; ewald_rtol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1e-05<br>

&nbsp;&nbsp; ewald_geometry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; epsilon_surface&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; optimize_fft&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>

&nbsp;&nbsp; ePBC&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = xyz<br>

&nbsp;&nbsp; bUncStart&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = TRUE<br>

&nbsp;&nbsp; bShakeSOR&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>

&nbsp;&nbsp; etc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Berendsen<br>

&nbsp;&nbsp; epc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Berendsen<br>

&nbsp;&nbsp; epctype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Isotropic<br>

&nbsp;&nbsp; tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20<br>

&nbsp;&nbsp; ref_p (3x3):<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 0]={ 1.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 1.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ref_p[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 1.00000e+00}<br>

&nbsp;&nbsp; compress (3x3):<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 0]={ 4.50000e-05,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 4.50000e-05,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; compress[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 4.50000e-05}<br>

&nbsp;&nbsp; andersen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 815131<br>

&nbsp;&nbsp; rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = PME<br>

&nbsp;&nbsp; rcoulomb_switch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; vdwtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Cut-off<br>

&nbsp;&nbsp; rvdw_switch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; epsilon_r&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; epsilon_rf&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; tabext&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; gb_algorithm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Still<br>

&nbsp;&nbsp; nstgbradii&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; rgbradii&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 2<br>

&nbsp;&nbsp; gb_saltconc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; implicit_solvent&nbsp;&nbsp;&nbsp;&nbsp; = No<br>

&nbsp;&nbsp; DispCorr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = No<br>

&nbsp;&nbsp; fudgeQQ&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.5<br>

&nbsp;&nbsp; free_energy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = no<br>

&nbsp;&nbsp; init_lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; sc_alpha&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; sc_power&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; sc_sigma&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.3<br>

&nbsp;&nbsp; delta_lambda&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; disre_weighting&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Conservative<br>

&nbsp;&nbsp; disre_mixed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>

&nbsp;&nbsp; dr_fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; dr_tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstdisreout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100<br>

&nbsp;&nbsp; orires_fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; orires_tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstorireout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100<br>

&nbsp;&nbsp; dihre-fc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; dihre-tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstdihreout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 100<br>

&nbsp;&nbsp; em_stepsize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0.01<br>

&nbsp;&nbsp; em_tol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>

&nbsp;&nbsp; niter&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 20<br>

&nbsp;&nbsp; fc_stepsize&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; nstcgsteep&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000<br>

&nbsp;&nbsp; nbfgscorr&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 10<br>

&nbsp;&nbsp; ConstAlg&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Lincs<br>

&nbsp;&nbsp; shake_tol&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1e-04<br>

&nbsp;&nbsp; lincs_order&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 4<br>

&nbsp;&nbsp; lincs_warnangle&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 30<br>

&nbsp;&nbsp; lincs_iter&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

&nbsp;&nbsp; bd_fric&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; ld_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1993<br>

&nbsp;&nbsp; cos_accel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; deform (3x3):<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 0]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 1]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; deform[&nbsp;&nbsp;&nbsp; 2]={ 0.00000e+00,&nbsp; 0.00000e+00,&nbsp; 0.00000e+00}<br>

&nbsp;&nbsp; userint1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userint2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userint3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userint4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userreal1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userreal2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userreal3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; userreal4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

grpopts:<br>

&nbsp;&nbsp; nrdf:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 51456<br>

&nbsp;&nbsp; ref_t:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 300<br>

&nbsp;&nbsp; tau_t:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>

anneal:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; No<br>

ann_npoints:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>

&nbsp;&nbsp; acc:&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>

&nbsp;&nbsp; nfreeze:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; N<br>

&nbsp;&nbsp; energygrp_flags[&nbsp; 0]: 0<br>

&nbsp;&nbsp; efield-x:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; efield-xt:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; efield-y:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; efield-yt:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; efield-z:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; efield-zt:<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; n = 0<br>

&nbsp;&nbsp; bQMMM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = FALSE<br>

&nbsp;&nbsp; QMconstraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; QMMMscheme&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

&nbsp;&nbsp; scalefactor&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1<br>

qm_opts:<br>

&nbsp;&nbsp; ngQM&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 0<br>

CPU=&nbsp; 0, lastcg= 8223, targetcg= 4112, myshift=&nbsp;&nbsp;&nbsp; 0<br>

nsb-&gt;shift =&nbsp;&nbsp; 1, nsb-&gt;bshift=&nbsp; 0<br>

Neighbor Search Blocks<br>

nsb-&gt;nodeid:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>

nsb-&gt;nnodes:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>

nsb-&gt;cgtotal:&nbsp; 8224<br>

nsb-&gt;natoms:&nbsp; 24594<br>

nsb-&gt;shift:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1<br>

nsb-&gt;bshift:&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0<br>

Nodeid&nbsp;&nbsp; index&nbsp; homenr&nbsp; cgload&nbsp; workload<br>

&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 24594&nbsp;&nbsp;&nbsp; 8224&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8224<br>

<br>

Max number of graph edges per atom is 4<br>

Table routines are used for coulomb: TRUE<br>

Table routines are used for vdw:&nbsp;&nbsp;&nbsp;&nbsp; FALSE<br>

Using a Gaussian width (1/beta) of 0.320163 nm for Ewald<br>

Cut-off's:&nbsp;&nbsp; NS: 1&nbsp;&nbsp; Coulomb: 1&nbsp;&nbsp; LJ: 1<br>

System total charge: -0.000<br>

Generated table with 1000 data points for Ewald.<br>

Tabscale = 500 points/nm<br>

Generated table with 1000 data points for LJ6.<br>

Tabscale = 500 points/nm<br>

Generated table with 1000 data points for LJ12.<br>

Tabscale = 500 points/nm<br>

Generated table with 500 data points for 1-4 COUL.<br>

Tabscale = 500 points/nm<br>

Generated table with 500 data points for 1-4 LJ6.<br>

Tabscale = 500 points/nm<br>

Generated table with 500 data points for 1-4 LJ12.<br>

Tabscale = 500 points/nm<br>

<br>

Enabling SPC water optimization for 7441 molecules.<br>

<br>

Will do PME sum in reciprocal space.<br>

<br>

++++ PLEASE READ AND CITE THE FOLLOWING REFERENCE ++++<br>

U. Essman, L. Perela, M. L. Berkowitz, T. Darden, H. Lee and L. G.

Pedersen <br>

A smooth particle mesh Ewald method<br>

J. Chem. Phys. 103 (1995) pp. 8577-8592<br>

-------- -------- --- Thank You --- -------- --------<br>

<br>

Center of mass motion removal mode is Linear<br>

We have the following groups for center of mass motion removal:<br>

&nbsp; 0:&nbsp; rest, initial mass: 150847<br>

There are: 24594 Atoms<br>

Removing pbc first time<br>

Done rmpbc<br>

Started mdrun on node 0 Tue Mar 14 05:52:54 2006<br>

Initial temperature: 430.169 K<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Step&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Time&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Lambda<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00000<br>

<br>

Grid: 10 x 10 x 20 cells<br>

Configuring nonbonded kernels...<br>

Testing AMD 3DNow support... not present.<br>

Testing ia32 SSE support... present.<br>

<br>

<br>

There are 24594 atoms in your xtc output selection

</body>

</html>