octahedral volume

From: Peter Varnai <Peter.Varnai_at_ibpc.fr>
Date: Fri 17 Jan 2003 11:26:58 +0100 (CET)

Dear AMBER list,

I would like to know the volume of a truncated octahedral cell used in a
simulation. The topology file gives the angle (109.471) and some length
parameter. Can someone tell me how this length parameter is related to
the truncated octahedron?
When constructing the system, leap reports a volume. This has
been changed between version 6 and 7 according to leap/tools.c:
         dVolume = dMaxX * dMaxY * dMaxZ;
                 dVolume *= 0.5;
                 dVolume *= 0.7698004;
 /* = sqrt(1 - 3*cost^2 - 2*cost^3), where cost = -1/3 = cos(109.471) */
Can someone explain where the ver7 formula comes from?
Thank you very much for your help.
Peter Varnai
Received on Fri Jan 17 2003 - 02:26:58 PST
Custom Search