[gmx-users] DrudeFF topology transfer from CHARMM to GROMACS
tomasz.k.piskorz at gmail.com
Thu Jun 1 14:54:22 CEST 2017
Thank you, Justin. That solves my problem.
Tomasz K. Piskorz
On 31 May 2017 at 13:51, Justin Lemkul <jalemkul at vt.edu> wrote:
> On 5/31/17 7:06 AM, Tomasz Piskorz wrote:
>> Dear all,
>> I would like to simulate cyclohexane using Drude FF in GROMACS. I've seen
>> that the topology is already available in CHARMM-DrudeFF and I would like
>> to transfer it to GROMACS. However, I don't know how to calculate the
>> charge of Drude particle, which is required by topology in GROMACS. Does
>> anyone know how to do it?
>> Atomic polarizability can be described by equation:
>> alpha = (q_D)^2/k,
>> where k = 418400 kJ/(mol*nm^2) and alpha is given in CHARMM topology (i.g.
>> for O in glycine alpha=0.000651 nm^3). However, when I try to calculate
>> charge I always get a wrong result. Is this question which I should use or
>> there is something more which I don't take under consideration?
> Yes, there is an implicit conversion factor of Coulomb's constant.
> qD = sqrt((alpha/kD)/CCELEC)
> In CHARMM units, this is
> qD = sqrt((0.651/1000)/332.0716)
> I'll add the model compounds to the next release of the GROMACS-formatted
> Drude force field. I have other issues to sort out before doing that,
> Justin A. Lemkul, Ph.D.
> Ruth L. Kirschstein NRSA Postdoctoral Fellow
> Department of Pharmaceutical Sciences
> School of Pharmacy
> Health Sciences Facility II, Room 629
> University of Maryland, Baltimore
> 20 Penn St.
> Baltimore, MD 21201
> jalemkul at outerbanks.umaryland.edu | (410) 706-7441
> Gromacs Users mailing list
> * Please search the archive at http://www.gromacs.org/Support
> /Mailing_Lists/GMX-Users_List before posting!
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or
> send a mail to gmx-users-request at gromacs.org.
More information about the gromacs.org_gmx-users