Re: AMBER: CHCl3 solvent box

From: David A. Case <case.biomaps.rutgers.edu>
Date: Mon, 5 Jan 2009 08:48:44 -0500

On Mon, Jan 05, 2009, WJ Ding wrote:

> I built a solvent box of chloroform with tleap. When I visualized the pdb file
> with VMD, it is found that the bond lengths of C-Cl even in the same molecule
> are different (1.72-1.83 angstrom, I just cheked three molecules), which will
> cause the C-Cl bond to be cleaved during the MD calculation. I do not think tt
> is reasonable. Is there anyone who ever encountered this problem? What's
> wrong with my solvent box?

Bonds can never be cleaved in MD simulations -- why not minimize or run
dynamics and see what happens.

You didn't say how to created the solvent box, but my guess is that you
used the CHCL3BOX unit. This file indeed has unequal C-Cl bond lengths
-- it could be a snapshot from an MD simulation, although the magnitude
of the deviations from ideality seems large.

...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 Tue Jan 06 2009 - 15:22:13 PST
Custom Search