From: Aron Broom (broomsday_at_gmail.com)
Date: Tue Mar 06 2012 - 08:44:57 CST
What if for the water molecules in your system you calculate from the
center of mass of the waters rather than just the atomic coordinates? I
don't know what the real answer is, but that would shrink your box ever so
slightly.
On Tue, Mar 6, 2012 at 2:48 AM, mjyang <mjyang_at_hku.hk> wrote:
> Dear NAMD users,
>
> I performed a simulation by NAMD for a protein system solvated by a
> cubic water box. Now I extracted the snapshots from NAMD trajectory for
> postpone analysis, in which the protein should be centered in the box. The
> box size was calculated by (max_x-min_x, max_y-min_y, max_z-min_z), where
> max_x,y,z and min_x, y, z are the maximum and minimum coordinate values of
> the extracted frame. However, the box size calculated in this way is a
> little larger than the real one, resulting in gaps in the recentered box.
> It seems that the box size written out by NAMD in the *.xsc file is better
> than what I calculated. But the NAMD output size (in *.xsc file) only
> corresponds to the last frame of a segment of trajectory file and I need
> the size for each frame in the traj. Could someone please tell me how NAMD
> calculate the box along MD simulations?
>
>
> Many thanks.
>
>
> Mingjun
>
>
-- Aron Broom M.Sc PhD Student Department of Chemistry University of Waterloo
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:17 CST