[gmx-users] System neutralization in Thermodynamic Integration

[Leandro Bortot]

Dear All,

     In order to calculate the binding free energy of a ligand to a protein
I am using Thermodynamic Integration, i.e. making the ligand "disappear".

     Because this ligand is charged (-1), at lambdas different than zero
(when the ligand is already fading) the total system charge increases until
it is +1 when lambda=1.

     I am using 1 nm cutoff for the non-bonded interactions and PME. Should
I make something to account for this total system charge? Won't it induce
artifacts in my calculation?


thanks in advance,
Leandro Bortot