[gmx-users] about "LINCS WARNING"
Kun-Lin Yang
klyang at che.wisc.edu
Tue Jan 14 04:35:13 CET 2003
Hi there:
I happened to simulate a similar problem with decanol and water interface, and I don't have any problems with mdrun. I guess you have to check your initial configuration (to see if there is any overlap or unwanted contact) rather than the parameters in em.mdp. The parameters there come with GROMACS are just fine. BTW, I noticed that you erased a line, define = -DFLEX_SPC, which may cause trouble too.
Kun-Lin, 1/13
Dear Sir or Madam:
I am trying to simulate the interface between alkyls & water, the decane (C10 , declared in ffgmx.itp) and SOL (OPLS) are used.
But I've met some problems:
1. mdrun always throw some error message, then stoped
2. how to enlarge "energy minimize " step length? I'd modified the "emstep ", but the mdrun shinked it rapidly.
Could you please give me some advice?
Thank you very much in advance!
rayn
Step 0, time 0 (ps) LINCS WARNING
relative constraint deviation after LINCS:
max 40.729252 (between atoms 1854 and 1855) rms 1.751049
bonds that rotated more than 30 degrees:
atom 1 atom 2 angle previous, current, constraint length
391 392 46.6 0.1530 0.1549 0.1530
392 393 45.9 0.1530 0.1517 0.1530
738 739 35.7 0.1530 0.1851 0.1530
739 740 90.0 0.1530 0.1464 0.1530
1245 1246 35.2 0.1530 0.1874 0.1530
.....
My em.mdp :
constraints = none
integrator = steep
dt = 0.0001
nsteps = 5000
nstlist = 0
ns_type = grid
rlist = 1.0
rcoulomb = 1.0
rvdw = 1.0
emtol = 10.0
emstep = 0.5
nstcgsteep = 100
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-users/attachments/20030113/d9904b9e/attachment.html>
More information about the gromacs.org_gmx-users
mailing list