[gmx-users] changing form of nonbonded force fields

[Andre Ghizoni]

Dear gromacs users,

I would like to know how and if it is possible to change the form of the 
lennard jones potential in gromacs.  What I want to do is shift the LJ 
function some R radius to the right for one species, i.e. f(r-R).  In effect 
my potential is f(r)=C1/(r-R)^12-C2/(r-R)^6, where r is the true distance 
between particles.  If it is not possible to do that, I would appreciate any 
suggestions on approximating this function with other functions that gromacs 
uses.

Thank you,
Andre

_________________________________________________________________
Stay in touch with old friends and meet new ones with Windows Live Spaces 
http://clk.atdmt.com/MSN/go/msnnkwsp0070000001msn/direct/01/?href=http://spaces.live.com/spacesapi.aspx?wx_action=create&wx_url=/friends.aspx&mkt=en-us