[gmx-users] Error in equilibration

s lal badshah shahbiochemist at yahoo.com
Tue Mar 25 05:56:54 CET 2008


Dear Justin & Mark,
Hi! 
I made changes in my pr.mdp file but still the errors came.Please guide me.
title           = pr.mdp
cpp             = /usr/bin/cpp
define          = -DPOSRES
constraints     = all-bonds
integrator      = md
dt              = 0.002
nsteps          = 100000
nstcomm         = 1
nstxout         = 250
nstvout         = 1000
nstlog          = 10
nstenergy       = 10
nstlist         = 10
ns_type         = grid
rlist           = 1.0
coulombtype     = PME
rcoulomb        = 1.0
vdwtype         = cut-off
rvdw            = 1.4
fourierspacing  = 0.12
fourier_nx      = 0
fourier_ny      = 0
fourier_nz      = 0
pme_order       = 4
ewald_rtol      = 1e5
optimize_fft    = yes
;Berendsen temperature coupling is on
Tcoupl          = berendsen
tau_t           = 0.1       0.1 
tc-grps          = protein  non-protein
ref_t           = 298       298 
;Pressure coupling is on
Pcoupl          = berendsen
tau_p           = 20
compressibility = 4.5e-5
ref_p           = 1.0
;Generate velocities is on at 298 K
gen_vel         = yes
gen_temp        = 298
gen_seed        = 173529



Step 0, time 0 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
max 224.227219 (between atoms 1932 and 1934) rms 4.182889
bonds that rotated more than 30 degrees:
 atom 1 atom 2  angle  previous, current, constraint length
    198    200   49.7    0.1449   0.1615      0.1449
    200    201   90.1    0.1090   1.2644      0.1090
    200    202   46.8    0.1529   0.1690      0.1529
    200    211   41.2    0.1522   0.1686      0.1522
    264    265   70.9    0.1229   0.1349      0.1229
    264    266   47.2    0.1335   0.1443      0.1335
    266    267   67.0    0.1010   0.1148      0.1010
----------------------------------------------------------------------------------
step 0Warning: 1-4 interaction between 196 and 201 at distance 1.104 which is larger than the 1-4 table size 1.000 nm
These are ignored for the rest of the simulation
This usually means your system is exploding,
if not, you should increase table-extension in your mdp file

Step 1  Warning: pressure scaling more than 1%, mu: -6.71376e+08 -6.71376e+08 -6.71376e+08

Step 1, time 0.002 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
max 10623.531250 (between atoms 1932 and 1933) rms 207.528702
bonds that rotated more than 30 degrees:
 atom 1 atom 2  angle  previous, current, constraint length
    170    172   63.2    0.1335   0.1515      0.1335
    172    173   41.5    0.1449   0.1423      0.1449
    172    181   87.7    0.1449   0.2689      0.1449
    173    174   38.5    0.1090   0.1402      0.1090
    173    184   33.3    0.1522   0.1855      0.1522
    175    178   64.7    0.1529   0.2016      0.1529
----------------------------------------------------------------------------------------
   2565   2576   59.4    0.1522   0.1316      0.1522
   2567   2568   38.4    0.1090   0.0818      0.1090
-------------------------------------------------------
Program mdrun, VERSION 3.3.1
Source code file: ns.c, line: 265

Fatal error:
Box was shifted at least 10 times. Please see log-file.
-------------------------------------------------------

Regards,
Lal badshah

 Send instant messages to your online friends http://uk.messenger.yahoo.com 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20080325/c12ecf3a/attachment.html>


More information about the gromacs.org_gmx-users mailing list