Re: AMBER: solvent box again

From: David A. Case <case.biomaps.rutgers.edu>
Date: Tue, 16 Dec 2008 15:02:05 -0500

On Tue, Dec 16, 2008, iccy liu wrote:
>
>
> I build a box size of 20Ax20Ax20A with my solvent molecule, and then run
> MD in NVT ensemble to relax it. I did min and heat it up to 300K step
> wise. Later I extracted my final structure, I found out that the shape of
> my solvent box is not a cub anymore.

As Ross said, you probably need to use the "image" command in ptraj.

> My question is now: how to solvate my system with my solvent box? what
> keyword to use? I am very new to Amber. In fact, I still did not figure out
> how to solvate my system with a single molecule as Professor Case mentioned
> in his reply.

The keyword is "solvateBox" or "solvateOct", described in section 3.4.38 of
the Users' Manual. The second argument can be a single molecule or a
pre-equilibrated box. You need to create a LEaP unit for either of these...
if you have done the MD with your solvent alone, you already shold be most of
the way there.

...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
      to majordomo.scripps.edu
Received on Wed Dec 17 2008 - 01:21:22 PST
Custom Search