hi, all

I want to reproduce a part work of a refernce paper which models CO2 in a slit containing two parallel SIO2 structure. The free energy difference for the interface CO2 molecule between the distances Z and Z0 are calculated by integrating the average force acting on the CO2 normal to the surface with coordinate Z. how can i obtain it or any other method to calculate free energy in MS?