[gmx-users] Adding 3 body force field term into Gromacs
lindahl at csb.stanford.edu
Wed Jun 2 19:50:29 CEST 2004
This has been covered before on the list, so you might get some more
info by using the search function at http://www.gromacs.org .
Adding an interaction isn't the problem, but rather how to keep track
of your interacting atoms. If you have a fixed list of indices that
interact it is quite easy to implement it as a new bonded interaction.
If not, the cost of finding all combinations of three atoms close in
space quickly becomes prohibitive, not to mention that you might need
several gigabytes of memory just for the neighborlists...
On Jun 2, 2004, at 4:38 PM, Tandia, Adama wrote:
> Dear ALL:
> Has anyone tried to add a tree body term in Gromacs?
> Any functional like Stillinger-Weber (SW) or Axilrod-Teller (AT)?
> SW: e*a*exp(b*(1/(r12/c-d)+1/(r13/c-d)))*(cos(theta)-cos(theta_0))**2
> where e,a,b,c,d,theta_0 are fixed
> Where Z_123 is fixed
> Thanks a lot,
> Adama Tandia
> Modeling & Simulation
> Corning INC
> Corning, NY 14831 USA
> Tel: 607 248 1036
> Fax: 607 974 3405
> gmx-users mailing list
> gmx-users at gromacs.org
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-users-request at gromacs.org.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Size: 2791 bytes
Desc: not available
More information about the gromacs.org_gmx-users