[AMBER] EEL become ****** in minimization

From: Suguru ASAI <suguruasai.gmail.com>
Date: Thu, 28 Jul 2016 23:04:56 +0900

Dear AMBER users,

# Summary

EEL value in mdout become ******* in minimization and a program failed.
I would like to solve this problem but don't know where to check or how
to fix.


# Detail

I've tried to use TIP4P-D water model (1) which is derived from
TIP4P/2005 water model (2).
So I've created frcmod and off files of both water models and performed
water only simulations.
However, EEL values in a simulation of TIP4P-D became ********** within
100 steps like:

    NSTEP ENERGY RMS GMAX NAME NUMBER
       1 4.7263E+08 2.1830E+08 6.9776E+09 O 569

  BOND = 0.1578 ANGLE = 0.0000 DIHED = 0.0000
  VDWAALS = ************* EEL = 737.3763 HBOND = 0.0000
  1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT =
0.0000


    NSTEP ENERGY RMS GMAX NAME NUMBER
     100 -1.1110E+08 1.1744E+12 4.3130E+13 H1 570

  BOND = 970.3621 ANGLE = 0.0000 DIHED = 0.0000
  VDWAALS = 2936.1737 EEL = ************* HBOND = 0.0000
  1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT =
1016.3076
  EAMBER = *************


    NSTEP ENERGY RMS GMAX NAME NUMBER
     200 -6.1725E+08 3.6244E+13 1.3291E+15 H1 570

  BOND = 970.3620 ANGLE = 0.0000 DIHED = 0.0000
  VDWAALS = 2936.1736 EEL = ************* HBOND = 0.0000
  1-4 VDW = 0.0000 1-4 EEL = 0.0000 RESTRAINT =
1016.3076
  EAMBER = *************

And sander/pmemd showed SEGV and stopped.
It did not happen with TIP4P/2005 and the differences between TIP4P-D
and TIP4P/2005 are

- ACOEFF (TIP4PD: 9.04E+05, TIP4P2005: 7.31E+05; Found in prmtop file)
- BCOEFF (TIP4PD: 8.99E+02, TIP4P2005: 7.36E+02; Found in prmtop file)
- HW charge (TIP4PD: 0.58, TIP4P2005: 0.5564; Written in OFF file)

So I'm guessing relatively large ACOEFF/BCOEFF of TIP4P-D would cause
the issue but I don't really know if it is or way to fix.
So I would be grateful if someone give me a direction to solve this.

In case if you need more details, I uploaded parameter files (frcmod,
off, PDBs) I used on

TIP4PD: https://www.dropbox.com/s/18e2297nowsp9dy/TIP4PD.zip?dl=0
TIP4P2005: https://www.dropbox.com/s/lg3g5oe7tb5i5gq/TIP4P2005.zip?dl=0


Thanks for reading.


Asai


1. Piana S, Donchev AG, Robustelli P, Shaw DE (2015) Water Dispersion
Interactions Strongly Influence Simulated Structural Properties of
Disordered Protein States. J Phys Chem B 119: 5113–5123.
2. Abascal JLF, Vega C (2005) A general purpose model for the condensed
phases of water: TIP4P/2005. J Chem Phys 123: 234505.


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Jul 28 2016 - 07:30:03 PDT
Custom Search