Hi,
> I have been trying to simulate a DNA chain and Na+ counterions (added using
> leap) by way of sander and the
> implicit GB approx.
> My problem is that on starting md to equilibrate the system the Na+ ions
> quickly fly away from the DNA chain.
> I have tried a couple of different protocols for equilibriation and
> minimisation to no avail.
> Is there some obvious reason why this is happening that I am overlooking?
...yes. It is rather easy to show that the (thermodynamic) average
distance between an ion and the solute is infinite. So what you see has
not only been observed by many (me included) but is also not unexpected.
One can try various things to change that (temporarily by somehow causing
stronger attachment between solute and ions) of "forever" by enclosing
both e.g. inside an elastic sphere. I guess both tricks won't lead to
new insights about molecules, but can be fun to implement :-)
good luck
Andreas
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Jul 23 2006 - 06:07:11 PDT
