On Mon, Aug 08, 2011, Maura Catherine Mooney wrote:
>
> I found out what was causing the 4 crystal waters to be treated as
> 'normal' bonds. Because they are part of the system I want to later
> study with QM, I initially ran gaussian calculations for the RESP calc
> in amber. The 4 waters had names such as O46, O47, O49, O49, H84,
> H85...etc. I looked at the solvent water which are named O and H1 and
> H2. So I changed the names of the these 4 cryst waters to O, H1 and
> H2. When I look at this in Xleap now, they appear as the triangluar
> water models, so it looks like they will be treated in the same manner
> as the solvent water. In saying this, because they are to be considered
> part of my QM system, I need to edit the charges on these waters to
> correspond with my calculated RESP charges. Is this valid to do so?
No, you should not need to edit the charges. Once you run the QM part,
whatever charges were in the MM section will be ignored. *Before* you get
to the QM calculations (e.g. for equilibration) you most likely want to
use a standard water model.
...good luck...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Mon Aug 08 2011 - 07:30:02 PDT