Re: [AMBER] Ptraj for diffusion coefficient of water

From: Ganesh Kamath <gkamath9173.gmail.com>
Date: Sun, 23 Jan 2011 13:01:36 -0500

Dear Dr. Case,

  Thank you for the suggestion of just doing NVE for the TIP3P water model.
I now get D=5.49*10-5 cm2/sec
which is comparable to literature.
I would also like to thank Lennart Nilsson for his useful comments.

Regards,
ganesh kamath
On Sat, Jan 22, 2011 at 5:39 PM, case <case.biomaps.rutgers.edu> wrote:

> On Sat, Jan 22, 2011, Ganesh Kamath wrote:
> >
> > I am using NVT ensemble for 901 molecules of water, and langevin
> > Temperature control at 298 K.
>
> If you want to compute diffusion constants, your production run (after NPT
> equilibration) should be NVE, that is no thermostat or barostat. The
> langevin
> thermostat could seriously affect the diffusion constant (and all other
> time-dependent properties).
>
> Otherwise, I don't see any problems in what you are doing.
>
> ...regards...dac
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Jan 23 2011 - 10:30:04 PST
Custom Search