Date: Tue, 12 Aug 2003 13:32:36 -0600 (Mountain Daylight Time)
From: "Thomas E. Cheatham, III"
Subject: Re: AMBER: question to calculate diffusion coefficient under periodic


> I was wondering if anyone is willing to share how to get the diffusion
> coefficient of particles when the periodic boundary conditions are
> included.

Following up on what David Case just mentioned, it is also possible to
unwrap a periodically imaged trajectory (if the time between frames is
short enough to not allow large scale translation of a molecule, i.e. >
1/2 box size between frames). In other words, if the simulation was run
with IWRAP=1, code can detect that a particle was imaged and unimage it
for the diffusion analysis. This is currently done in ptraj (ONLY for
orthorhombic i.e. 90x90x90 unit cells) within the diffusion command. To
make this more general, ptraj would have to modified to create an unimage
command or the diffusion command modified to allow for general triclinic
unit cells.

Thomas E. Cheatham, III (Assistant Professor) College of Pharmacy