- Contemporary messages sorted: [ by date ] [ by thread ] [ by subject ] [ by author ] [ by messages with attachments ]

From: Lucas Bandeira <bandeiralucas97.gmail.com>

Date: Thu, 20 Dec 2018 08:21:04 -0200

Hello Saeed

Thank you very much for your help.

Yes, I'm using a NPT ensemble to calculate self-diffusion of water. I'm

glad to know that thermostat can cause this kind of problem. I'll try

your suggestion.

Thank you very much.

Lucas

On Wed, Dec 19, 2018 at 6:52 PM Saeed Izadi <izadi.vt.edu> wrote:

*> Hi Lucas,
*

*>
*

*> What type of ensemble are you using for the calculation? .. note that
*

*> unlike many other properties of water which are calculated in NPT ensemble,
*

*> self-diffusion coefficient is safer to be estimated with no thermostat
*

*> (i.e., NVE ensemble). The use of a thermostat affects how the water
*

*> molecules move, usually leading to under estimation of self-diffusion
*

*> value. To make sure the system’s temperature does not deviate from the
*

*> target temperature, it is suggested to run a large number of short (NVE)
*

*> simulations starting from well- equilibrated evenly spaced NPT-sampled
*

*> coordinates and velocities as initial conditions. The Einstein relation can
*

*> then be used to infer the diffusion constant from the time dependence of
*

*> the mean-squared displacement, and then final value of the self-diffusion
*

*> coefficient can be obtained by averaging over all the NVE trajectories. You
*

*> can also take a look at the following paper for more information:
*

*> https://onlinelibrary.wiley.com/doi/full/10.1002/wcms.1347
*

*>
*

*> best,
*

*> Saeed
*

*>
*

*> On Wed, Dec 19, 2018 at 9:14 AM Lucas Bandeira <bandeiralucas97.gmail.com>
*

*> wrote:
*

*>
*

*> > Hello,
*

*> >
*

*> > I'm trying to calculate the self coefficient diffusion of water using the
*

*> > FF TIP3P. I know that to TIP3P the value of the diffusion coefficient of
*

*> > water is different than the obtained experimentally, 5.19 * 10⁻⁹ m²*s⁻¹
*

*> > instead of 2.3 * 10⁻⁹ m² * s⁻¹, but when I try to calculated the
*

*> diffusion
*

*> > coefficient I get the value 4.21 * 10⁻⁹ m² * s⁻¹. I'm using the formula
*

*> on
*

*> > AMBER Manual:
*

*> >
*

*> > [image: image.png]
*

*> >
*

*> > Using n = 3. The MSD value I compute using the out_r.xmgr generated when
*

*> I
*

*> > use cpptraj.
*

*> > My input looks like:
*

*> >
*

*> > diffusion :WAT.O 0.2 out
*

*> >
*

*> > I'm using Amber14 so I can't calculate diffusion coefficient directly.
*

*> >
*

*> > Thank you in advance,
*

*> >
*

*> > Lucas
*

*> > _______________________________________________
*

*> > 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
*

*>
*

_______________________________________________

AMBER mailing list

AMBER.ambermd.org

http://lists.ambermd.org/mailman/listinfo/amber

Received on Thu Dec 20 2018 - 02:30:01 PST

Date: Thu, 20 Dec 2018 08:21:04 -0200

Hello Saeed

Thank you very much for your help.

Yes, I'm using a NPT ensemble to calculate self-diffusion of water. I'm

glad to know that thermostat can cause this kind of problem. I'll try

your suggestion.

Thank you very much.

Lucas

On Wed, Dec 19, 2018 at 6:52 PM Saeed Izadi <izadi.vt.edu> wrote:

_______________________________________________

AMBER mailing list

AMBER.ambermd.org

http://lists.ambermd.org/mailman/listinfo/amber

Received on Thu Dec 20 2018 - 02:30:01 PST

Custom Search