On Thu, May 27, 2004, Tomas Kubar wrote:
>
> I have been trying to simulate the system of two methane molecules
> dissolved in a box of TIP3P water. To build this system, I attempted to
> use the xleap program from the AMBER 8 package. A problem occured when I
> wanted to solvate the two methanes:
>
> > solvatebox CH4 TIP3PBOX 10
>
> The program created a box and filled it with water, but the resulting
> density was 2.2 g/cc. An attempt to equilibrate such a system led to an
> error - excessive number of nonbonded pairs.
>
Works for me, unless there is something unusual about your CH4 unit. I
tried the following:
Welcome to LEaP!
(no leaprc in search path)
> source leaprc.ff94
> solvatebox ACE TIP3PBOX 10
Solute vdw bounding box: 6.053 5.036 4.754
Total bounding box for atom centers: 26.053 25.036 24.754
Solvent unit box: 18.774 18.774 18.774
Total vdw box size: 28.912 28.370 27.815 angstroms.
Volume: 22814.792 A^3
Total mass 8402.468 amu, Density 0.612 g/cc
Added 464 residues.
You could try the above, just in case there is some special about your CH4
unit. Other than that, I don't have any ideas about what might be going on.
You could compare the box size and number of waters you get with numbers like
those shown above; maybe that would provide some clue.
....dac
--
==================================================================
David A. Case | e-mail: case.scripps.edu
Dept. of Molecular Biology, TPC15 | fax: +1-858-784-8896
The Scripps Research Institute | phone: +1-858-784-9768
10550 N. Torrey Pines Rd. | home page:
La Jolla CA 92037 USA | http://www.scripps.edu/case
==================================================================
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Thu May 27 2004 - 16:53:01 PDT
