Thank you Mark What do you think about making a rerun on the trajectories generated previously with PME but this time using coulombtype: cut-off? Could you suggest a cut off value? Thanks again Williams