From: Norman Geist (norman.geist_at_uni-greifswald.de)
Date: Tue Jul 10 2012 - 00:44:57 CDT
I think there's no SHAKE like feature implemented for such usage. It's also
kind of a different simulation method, is mostly called "rigid body".It is,
as it seems, hard to combine with all atom force field, so there are no
implementations I know in no classical md software.
Possibly you could use the extrabonds command to introduce some strong
forced angles/torsions/bonds to your ions to prevent them from moving to
much within, this shouldn't effect interaction with other molecules and
could be worth a try.
Von: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] Im Auftrag
von Evgeny Bulat
Gesendet: Montag, 9. Juli 2012 20:56
Betreff: namd-l: Fixing bonds in dihydrogen phosphate to save computation
Recently, I have explored using dihydrogen phosphate, H2PO4-, for ionizing
my biomolecular MD simulation. The reason is that H2PO4- is more
physiological than, say, Cl-. After generating H2PO4- toppars and ionizing
the system, I performed 10ns of equilibration and noticed that, not
surprisingly, the computation time is significantly slower than when Cl- is
used; in fact, it's just not feasible at the moment.
At this point, I'm looking into somehow freezing most bonds within H2PO4-,
essentially modeling each anion as a rigid structure. I figure this may make
the computation more permissible. I have looked into how NAMD can be used to
freeze bonds between hydrogen and heavy atoms (with SHAKE), and wonder if a
similar method has been developed (or can be developed) for specific bonds
Thanks a bunch for any feedback.
Dept. of Bioengineering, UC Berkeley
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:45 CST