[gmx-users] widom TPI insertion

[deepak bapat]

Dear Gromacs Users

I have carried out TPI inserion of small LJ particle in SPC/E water.
NPT ensemble was used. 1 bar, 300 K.
3 nm box was used for simulation without test particle.
LJ parameters of the TPI particle were 9e-02 and 6.51e-02.
Successful insertion probability is around 18%.

tpi.xvg file reports 7 datasets.
The first data set is -kT log(<Vexp(-bU)>/<V>).
avg value of this series is around 1.15 kJ/mol.nm3 +/- 3e-4. (obtained from
block averaging)
How should I convert this value to excess chemical potential in kJ/mol
units considering volume effects?

Regards
-- 
Deepak U. Bapat