On Thu, Jul 21, 2022, Amit, Amit (amit.amit--- via AMBER wrote:
>Please help me for : How to avoid SHAKE on a specific QM atom in QM/MM
>simulations?
>
>
>I tried "noshakemask=.12" flag but than sander asked for "qmshake=0". I
>don't want to avoid SHAKE for all of my QM atoms, instead I want to avoid
>SHAKE only for a specific H atom, that is being transferred from atom A to
>B. Aim is to capture the energy barrier for hydride transfer, and for that
>I don't want to restrain this particular bond involve this H atom.
The manual implies, and the code agrees, that shaking hydrogens in the qm
region is an all-or-nothing proposition. You can't shake most of them, but
leave one out.
You could play with the code starting at line 782 of set.F90, removing the
test that qmmm_nml%qmshake == 0 (lines 798 and 806) and see if bad things
happen. Try a small system first, and critically examine what is going on.
Or, be conservative: if have a hydrogen that is being transfered, you might
want to have a 1 fs time step anyway, since the forces near the trasnfer
point might be changing very quickly.
...good luck...dac
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Aug 04 2022 - 13:35:19 PDT
