> Can anyone suggest me a good way to equilibrate a box of Lennard-Jones
> particles at 300K and 1atm starting with certain lattice conformation? I
> tried minimization and a long NVT MD. But as long as I switched to NPT,
> the system just blew up after some picoseconds. Thanks.
Maybe the lattice conformation is not stable under NPT, i.e. that at 300K
it wants to melt. You have to be careful to have a good balance with the
LJ particles. I would suggest either increasing the well depth or running
at higher pressure with a long relaxation time. You could also try
running the lattice in vacuum to see if it still falls apart (my guess is
that it will). It seems unlikely to be an artifact of the pressure
coupling.
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu Sep 02 2004 - 01:15:11 PDT
