On 2012-10-15 16:53, Peter C. Lai wrote:
Isn't Mark's script outdated for this purpose?
I just uploaded a new script that seems to work:
http://www.gromacs.org/@api/deki/files/185/=charmm2gromacs-pvm.py
Please check the comment on
http://www.gromacs.org/Downloads/User_contributions/Other_software
charmm forcefields specify epsilon and sigmas so you only need to convert
them:
gromacs(sigma) = (charmm(Rmin/2)/10) * (2/(2^(1/6)))
gromacs(epsilon) = charmm(eps) * 4.184
For 1-4 pair interactions,
gromacs(sigma1,4 i,j) = (charmm(Rmin/2_1,4_i) + charmm(Rmin/2_1,4_j))/(2^1/6)
gromacs(eps1,4 i,j) = sqrt(charmm(eps_1,4_i) * charmm(eps_1,4_j)) * 4.184
On 2012-10-15 07:27:51AM -0700, spin wrote:
Hello, everyone.
I used the Mark's script's to convert the CGenff (version 2b7 ) parameter
file to Gromacs .itp files. In the ffcharmmnb.itp, the script gave the c6
and c12, while the charmm's ffnonbond.itp showed epsilon and sigma.I do not
understand the relation between the c6/c12 and epsilon and sigma, and I
have a poor Perl skill. Can someone give me a solution? In addition, the
script makes all atoms' charge zero in the file, which is not the case in
the ffnonbond.itp. Why is it?
Thank you!
Qing Liu
--
View this message in context:
http://gromacs.5086.n6.nabble.com/The-problem-of-converting-CGenff-parameters-to-those-of-Gromacs-tp5002042.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
--
David van der Spoel, Ph.D., Professor of Biology
Dept. of Cell & Molec. Biol., Uppsala University.
Box 596, 75124 Uppsala, Sweden. Phone: +46184714205.
sp...@xray.bmc.uu.se http://folding.bmc.uu.se
--
gmx-users mailing list gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists