AMBER Archive (2007)

Subject: Re: AMBER: perturbed atoms are not bonded on using TIP5P water

From: David A. Case (case_at_scripps.edu)
Date: Mon Jan 15 2007 - 20:23:17 CST


On Tue, Jan 16, 2007, sugino tatsuro wrote:
>
> > edit model ( and select a water molecule and set perturbation by GUI )

Although this would seem to make sense, I'm not sure LEaP works that way.

Try something like this (untested!)

source leaprc.ff99SB
loadamberparams frcmod.tip5p
MTP = copy TP5
edit MTP (and make the perturbation changes; check with "desc" that
            the MTP unit has the proper bonds.)
saveoff MTP mtp.off (in case you need it later, or possibly to modify it
                      by hand).

solvatebox MTP TP5 20.0 0.8
saveAmberParm .....
quit

This procedure helps to make sure that the special (perturbed) water uses
parameters different from the regular tip5p water. Just choosing one water
visually to perturb might be getting everything mixed up.

...hope this helps...dac

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu