To gain an insight into the thermodynamics associated with the removal of molecule from the crystal surface, a free energy landscape for each surface/solvent pair is calculated from a number of independent WTmetaD simulations using Mean Force Integration. By implementing a general approach rooted in Thermodynamic Integration, MFI allows us to consistently use sampling obtained in independent biased simulations to estimate a joint FES, as discussed in ref. The FES describes the free energy associated with a detachment event and the basins recovered in CV space correspond to high probability configuration visited in the pathway towards detachment. The FES are converged up to the barrier associated with detachment from the surface.