No steric clashes, still not minimizing

From: Francesco Pietra (chiendarret_at_gmail.com)
Date: Wed Mar 11 2015 - 06:38:40 CDT

Hallo:
Minimization of the non solvated system under GB conditions runs fine. In
contrast, the TIP3 solvated/ion-added system resists
minimization/velocityQuenching .

No indication of clashes:

lsort -unique -increasing -integer [join [measure contacts 1.4 [atomselect
top "noh"]]]
>Main< (psfgen) 36 %

It is only at 1.5A that the above command reveals a two atoms.

However, at
timestep 0.000001
maximumMove 0.000001

with
rigidBonds water
rigidTolerance 0.000001
exclude scaled1-4
cutoff 12.0
switching on
switchdist 10.0
pairlistdist 13.5 # cutoff +3.5
margin 5
1-4scaling 1.0

error arises:

TCL: Running for 1000 steps
FATAL ERROR: Low global exclusion count! (97807 vs 97818) System unstable
or pairlistdist or cutoff too small.

FATAL ERROR: See http://www.ks.uiuc.edu/Research/namd/bugreport.html
------------- Processor 0 Exiting: Called CmiAbort ------------
Reason: FATAL ERROR: Low global exclusion count! (97807 vs 97818) System
unstable or pairlistdist or cutoff too small.

FATAL ERROR: See http://www.ks.uiuc.edu/Research/namd/bugreport.html

Charm++ fatal error:
FATAL ERROR: Low global exclusion count! (97807 vs 97818) System unstable
or pairlistdist or cutoff too small.

Thanks for suggestions
francesco Pietra

This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:20:57 CST