Hi Yao,
As far as I know, current Amber doesn't have this kind of facility (?)
If you are only interested with the electrostatic QM/MM interaction
energy, one simple way you could try is to modify the prmtop file by
zero the charge of the atoms within the MM part that you are
interested in. Then calculating the QM/MM interaction energy with the
original prmtop and with the modified prmtop file separately. The
difference between two interaction energy will be the contribution
(approximately) of the zero charged MM part.
Or, if you don't mind to tweak the source code a little bit. You can
first save the converged density matrix with original prmtop file.
Then do the second calculation with the saved density matrix after
zero the charges. This way, from my personal experience, you may get
better result.
Best,
Qiantao
On Fri, Oct 7, 2011 at 3:27 PM, Yao Yao <yao0o.ymail.com> wrote:
> Hi,
>
> I am running a QMMMMD simulation. For some reason, I wanna check the interaction energy between different parts
>
> of my system, like solutes and solvent, or QM part and MM part. Is there any way to do this in Amber?
>
> Thanks,
>
> Yao
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Oct 07 2011 - 13:30:02 PDT
