VMD-L Mailing List

From: David Sáez (davidsaezsan_at_gmail.com)
Date: Wed Jan 29 2020 - 07:57:40 CST

Hi James and Giuseppe,

Thanks for your suggestions. I used the charges provided by CGenFF (much
less polarized than those obtained in the calculation) and the profile
improved a lot:

https://www.dropbox.com/s/ezw6ulylgpa4sxp/Selecci%C3%B3n_025.png?dl=0

The MM potential in the picture has an RMSE of 0.875. Is that a reasonable
value, or should I keep searching a profile more similar to the QME?

PS: Regarding the charge optimization I noticed that in different
iterations I obtained quite different values for the atomic charges,
without signs of convergence. What is the target function in the process?
Reproducing the electrostatic potential like in ESP?

Regards.

On Mon, Jan 27, 2020 at 10:52 PM Gumbart, JC <gumbart_at_physics.gatech.edu>
wrote:

> Hi David,
>
> You should look carefully at the dihedral scan that gives the anomalously
> low values. There is probably some interaction that is too favorable,
> possibly due to unrealistic charges. If you suspect that, you can modify
> the charges to see if it makes a difference.
>
> Best,
> JC
>
> On Jan 25, 2020, at 12:07 AM, David Sáez <davidsaezsan_at_gmail.com> wrote:
>
> Hello everyone,
>
> I am working on the parametrization of a small molecule using the superb
> tool called FFToolkit and Gaussian 16. Up to now, it has worked great. Now,
> in the dihedral section I am facing problems. When trying to find
> parameters to reproduce the QM potential, I get this plot:
>
> https://www.dropbox.com/s/w0h6t27vs0z05ru/Selecc_024.png?dl=0
>
> Is it normal for the preliminary MM potential to be so high in comparison
> to the QM potential? Should I take a look on the other elements of
> interaction (charges, angles, bonds) because they could be raising the
> energy too much?
>
> Thanks for your time.
>
> David.
>
>
>