Multiwfn forum / Esp grid points

Re: Esp grid points

dummy@example.com (julieborah66) — Sun, 14 Oct 2018 03:53:50 +0000

Okay Sir, thankyou Sir for updating

Re: Esp grid points

dummy@example.com (sobereva) — Fri, 12 Oct 2018 19:13:18 +0000

julieborah66 wrote:

Thank you Sir for being so prompt in replying to queries. I do need an option where I can calculate the RMSE error of a specific atom in a given molecule, if I upload the chg file. It will be so useful if such an option will be available in Multiwfn

Today I have updated Multiwfn 3.6(dev) on the official site, the function you need has been supported. Please check Section 4.7.8 of the latest version of manual on how to use this function.

Re: Esp grid points

dummy@example.com (sobereva) — Wed, 10 Oct 2018 20:39:37 +0000

julieborah66 wrote:

Thank you Sir for being so prompt in replying to queries. I do need an option where I can calculate the RMSE error of a specific atom in a given molecule, if I upload the chg file. It will be so useful if such an option will be available in Multiwfn

I will update Multiwfn soon, new version will include this option. Once update is finished, I will inform you here.

Re: Esp grid points

dummy@example.com (julieborah66) — Wed, 10 Oct 2018 16:16:21 +0000

Thank you Sir for being so prompt in replying to queries. I do need an option where I can calculate the RMSE error of a specific atom in a given molecule, if I upload the chg file. It will be so useful if such an option will be available in Multiwfn

Re: Esp grid points

dummy@example.com (sobereva) — Wed, 10 Oct 2018 13:38:48 +0000

julieborah66 wrote:

Hello Sir,
Apologies to disturb you once again, is there any way out in Multiwfn to calculate the RMSE error around a specific atom?

There is no such an option. If you really need it, I will consider add such a function in the MK module

Re: Esp grid points

dummy@example.com (sobereva) — Wed, 10 Oct 2018 13:32:36 +0000

julieborah66 wrote:

Sir,
I uploaded a file containing CM5 charges. The phenol.wfn file was taken from the example section and then the CM5 charges were calculated. I uploaded the .chg file in the CHELP module I got RMSE error 0. I also did it with MK scheme. I got RMSE. Why is that so? Is there something wrong in my procedure. I am attaching a screenshot of the results
https://i.postimg.cc/VJ82WYVz/SCREENSHOT.png

You didn't use the latest version, the CHELPG and MK interfaces of the latest version of 3.6(dev) on the Multiwfn website look like this

------------ Calculation of CHELPG charges ------------
-3 Toggle using atomic charges in external file, current: No
-2 Toggle loading additional fitting centers from external file, current: No
-1 Toggle using fitting points in external file, current: No
0 Return
1 Start calculation!
2 Set grid spacing, current: 0.567 Bohr ( 0.300 Angstrom)
3 Set box extension, current: 5.291 Bohr ( 2.800 Angstrom)
5 Choose ESP type, current: Nuclear + Electronic
6 Toggle if exporting fitting points with ESP after the task, current: No

You should use a file containing wavefunction information (e.g. wfn, fch) as input file, and then enter this module. Before calculation via option 1, you should first select option -3 to load atomic charges from specified .chg file.

Re: Esp grid points

dummy@example.com (julieborah66) — Wed, 10 Oct 2018 12:08:57 +0000

Sir,

I uploaded a file containing CM5 charges. The phenol.wfn file was taken from the example section and then the CM5 charges were calculated. I uploaded the .chg file in the CHELP module I got RMSE error 0. I also did it with MK scheme. I got RMSE. Why is that so? Is there something wrong in my procedure. I am attaching a screenshot of the results

Re: Esp grid points

dummy@example.com (julieborah66) — Wed, 10 Oct 2018 11:43:50 +0000

Hello Sir,

Apologies to disturb you once again, is there any way out in Multiwfn to calculate the RMSE error around a specific atom?

Re: Esp grid points

dummy@example.com (sobereva) — Mon, 08 Oct 2018 08:56:54 +0000

Please check end of Section 4.7.8 of the latest version of Multiwfn manual, in which I exemplified how to evaluate reproducibility of ESP (measured by RMSE and RRMSE) of given atomic charges. Note that only currently latest version of Multiwfn has this capacity.

Re: Esp grid points

dummy@example.com (julieborah66) — Mon, 08 Oct 2018 08:18:02 +0000

Hello Sir,
Thanks for your reply, Actually, I am referring to a paper " Charge Anisotropy: Where Atomic Multipoles Matter Most: Christian Kramer, Alexander Spinn, and Klaus R. Liedl, Journal of Chemical Theory and Computation, 2014, 10, 4488?4496" in the page 4490 of statistical analysis they mention a line " Each grid point is assigned to the nearest atom for the analysis " (grid points at a distance of 1.66?2.2 times the vdW radius)

How can I assign in my case if I want to calculate the RMSE of the ESP between the quantum and the CM5 charges. I want not using any fitting CHELP or RESP. It's just that I want to calculate to error metrics.

Re: Esp grid points

dummy@example.com (sobereva) — Fri, 28 Sep 2018 00:39:30 +0000

Dear julieborah66,

There is no any correspondence between ESP fitting points and atoms. In principle, ESP value at all points affect charge of all atoms. If you are confused, please check Section 3.9.10 for introduction of CHELPG type of ESP fitting method.

Best regards,

Tian

Esp grid points

dummy@example.com (julieborah66) — Thu, 27 Sep 2018 14:04:23 +0000

Sir,
I generated esp (in grids) for my molecule of atomic charges and got it into a text file espfit.txt?
Is there a way to know which atom is assigned which grid points, since the text file contains x,y,z coordinate of the grid points and fourth column has esp but how can I know the atom numbers?