[gmx-users] G_ENERGY COMMAND

asaffarhi at post.tau.ac.il asaffarhi at post.tau.ac.il
Fri Jan 13 10:59:28 CET 2017

Dear Gromacs Users,

We would like to introduce a recently published article on calculation  
of relative free energies.


The article includes a decoupling analysis in which the partition  
functions of the transformed molecules are decomposed into two  
partition functions. It enables to keep dihedral angle terms and  
coupled bond angle terms at the separation point on, independently of  
the function type, also when there is a bond junction at the  
separation point between the shared and unshared submolecules. The  
decoupling analysis is demonstrated in a relative solvation free  
energy calculation of p-Cl and p-CH3 and the calculation is shown to  
be more efficient since less terms are removed.
In addition, it is proved analytically that when capping the  
non-bonded potentials (soft core technique) the integrated function  
(in TI) is monotonic. It is also shown mathematically that when  
capping the non-bonded potentials and the two molecules are simulated  
separately, the integrated function can be non-steep. It is also  
explained how when the systems have rugged energy landscape they can  
be equilibrated in the same sampling dimension.

Other important references there:
Ref. 34 (33), where it has been suggested to simulate the two  
molecules separately and to cap the potential with accessible energies  
of ~5kcal/mol (page 28, denoted by E_cutoff).
Ref. 35 where it has been suggested to remove terms only of the atoms  
of the unshared submolecules (v1 Eq. (21), v6 Fig. 7) and to cap  
potential terms with accessible energies of 7kcal/mol (v1 Eq. (24)).
Ref. 38 which is about an exact calculation of molecular free energies  
in a general environment. The calculations include the free energies  
of bond stretching, bond angle, dihedral angle, bond junctions, and  
complex structures. These free energies are also relevant for  
restraints in binding.

Thank you for your attention,
Best regards,

More information about the gromacs.org_gmx-users mailing list