Dear amber users,
I have an equilibrated periodic system (truncated octahedron box) and I
want to do a very small modification of the system using Leap. In the
meantime I want to retain all info regarding periodicity (box shape and
dimensions). What I do is: make a pdb file from the top file and the
last restart file (after performing imaging in ptraj), load the new pdb
into leap, modify what I need to modify, and then save a new top and crd
file. Now, the problem is that I want to retain the box dimensions and
shape of the system into the new top and crd files. For that, if the box
was cubic, I always used setbox command in Leap to fix my box dimensions
to the desired values and to have box info in the new top file .... Now,
the problem is that setBox sets automatically a cubic box even if I give
it the angle values of 109..., it still sets them to 90.
Does anyone have any idea how to overcome this problem?
Thanks
vlad
P.S. The modification is replacing one -OH group with a -H atom in my
system.
--
Dr. Vlad Cojocaru
EML Research gGmbH
Molecular and Cellular Modeling Group
Schloss-Wolfsbrunnenweg 33
69118 Heidelberg, Germany
Phone: +49-6221-533266
Fax: +49-6221-533298
e-mail: Vlad.Cojocaru.eml-r.villa-bosch.de
http://projects.villa-bosch.de/mcm/people/cojocaru/
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Wed Sep 13 2006 - 06:07:09 PDT