Previous Thread
Next Thread
Print Thread
Units of the parameters in CHARMM36 forecefiled
#37384 03/15/19 04:55 PM
Joined: Mar 2019
Posts: 11
M
Mohsen Offline OP
Forum Member
OP Offline
Forum Member
M
Joined: Mar 2019
Posts: 11
Dear all,

Hi!

Where can I find units of the parameters in the CHARMM36? For instance, what are the units of the sigma and epsilon in the ffnonbonded.itp file.

Thank you
Mohsen.

Re: Units of the parameters in CHARMM36 forecefiled
Mohsen #37386 03/15/19 05:50 PM
Joined: Sep 2003
Posts: 8,485
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,485
Most of the units are in the "usage" doc file on the charmm.org web site. Energies are typically in kcal/mol, and distances are in ┼nstroms.

The actual CHARMM force field parameter files have names such as par_all36_prot.prm, so it seems you are using a derived file in some other package. The distributed files tabulate Rmin/2 instead of sigma for the nonbond terms.


Rick Venable
computational chemist

Re: Units of the parameters in CHARMM36 forecefiled
Mohsen #37388 03/15/19 06:17 PM
Joined: Mar 2019
Posts: 11
M
Mohsen Offline OP
Forum Member
OP Offline
Forum Member
M
Joined: Mar 2019
Posts: 11
Thanks for your kind answer

I have downloaded CHARMM36 from this webpage:

charmm36-nov2018.ff.tgz in

http://mackerell.umaryland.edu/charmm_ff.shtml

So do you mean that this force field is a derived file?
So finally what are the units of the parameters of sigma and epsilon? Also, do you think that tabulating Rmin/2 instead of sigma is a correct model?

Re: Units of the parameters in CHARMM36 forecefiled
Mohsen #37390 03/15/19 06:51 PM
Joined: Sep 2003
Posts: 8,485
rmv Offline
Forum Member
Offline
Forum Member
Joined: Sep 2003
Posts: 8,485
That is a conversion of the distributed CHARMM files into GROMACS format, so presumably the change of units and computing sigma from Rmin/2 has already been performed. The units should conform to those used by GROMACS.

Sigma and Rmin/2 can be readily interconverted, but early CHARMM developers chose Rmin/2 to facilitate computing the atom pair values from the simple sum of the self terms (arithmetic combining rule).


Rick Venable
computational chemist


Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 5.6.33-0+deb8u1 Page Time: 0.011s Queries: 22 (0.005s) Memory: 0.9122 MB (Peak: 0.9941 MB) Data Comp: Off Server Time: 2020-08-14 13:57:44 UTC
Valid HTML 5 and Valid CSS