[gmx-developers] editing the code of the VDW and electric potentials and forces

[Feenstra, K.A.]

Op 30 aug. 2013 om 08:49 heeft "Asaf Farhi" <asaf.farhi at weizmann.ac.il<mailto:asaf.farhi at weizmann.ac.il>> het volgende geschreven:

Dear GROMACS Developers

My name is Asaf.
I'm trying to implement a novel free energy calculation method and interested to do a very minor code change.
I will really appreciate help on this issue.
What I'm trying to implement in GROMACS is very simple. It's actually only one if statement that is relevant to the VDW and Electric energies and forces.

What is needed is to add for every calculation of energy and force:
If (E_VDW>Ec)
{
E_VDW=Ec;
f_VDW=0;
}
If (E_EL>Ec)
{
E_el=Ec;
f_el=0;
}

where E_VDW is the energy associated with the pair VDW interaction, f_VDW is the force associated with the pair VDW potential and Ec is a constant.
In addition I need to multiply specific VDW and electrical forces by a const (that are associated with certain atoms).

Could anyone advise me how to do this code change?

If I understand what you want to do, i.e. a fixed upper limit to the energies, than you could do that easily, without touching the code, bu using interaction tables. You will find how in the manual.

What you haven't explained is what you want to achieve.

The setup you describe will give an unphysical discontinuity in the forces, which will give rise to some strange artifacts in your simulations for sure. I expect these may be similar, in type and severity, to those associated with the simple interaction cut-off schemes. Using that as a basis for 'a novel free energy calculation method' seems rather unlikely to produce valid results.


Thanks in advance,
Best regards,
Asaf
--
gmx-developers mailing list
gmx-developers at gromacs.org<mailto:gmx-developers at gromacs.org>
http://lists.gromacs.org/mailman/listinfo/gmx-developers
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-developers-request at gromacs.org<mailto:gmx-developers-request at gromacs.org>.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://maillist.sys.kth.se/pipermail/gromacs.org_gmx-developers/attachments/20130830/be3c5302/attachment.html>