I recently made an RHDO of water with A=36, and have posted my script to the archive.

In brief, I create a sphere with the same volume (and therefore number of water molecules) by cutting a spherical octant (1/8th) from a box and rotating 3 times to make a sphere. This can be used in general to make a sphere of radius R from a cube with an edge greater than R.

Since the volumes are equivalent (by design) and RHDO is nearly spherical, the sphere can be converted to RHDO by simply imposing the crystal lattice and minimizng. I then equilibrated with a few 100 ps of dynamics, before using it as a solvent bath for some simulations.

By using the volume formula for RHDO in crystl.doc and that for the volume of a sphere, you can establish the sphere radius needed for an RHDO with any value of A.

*Edited by rmv (04/03/04 05:52 PM)*
_________________________

Rick Venable

computational chemist