Dear GMXers,
I'm trying to implementing a flexible water model in "Gordillo, M.C.; Marti,
J. *J. Phys.: Condens. Matter* *2010*, *22*, 284111" using GROMACS.
The intra-molecular water potential has the following form:
V=D{[1-exp(-alpha*r1)]^2 + [1-exp(-alpha*r1)]^2} + k_a /2 *r3^2 +
k_a*r3*(r1+r2) + k_rr *r1*r2
where r1, r2 is the O-H bond distance and r3 is the H-H distance.
The first term can be modeled using Morse potential.
The second term can be model by harmonic potential.
Could anyone give some suggestions about how to implement the remaining two
coupling terms in GROMACS? Or can anyone recommend other packages that
support these function forms?
Thank you.
Best,
Yanbin
--
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
