If you have done SMD using NAMD, the tutorial itself has a tcl script for calculating PMF. I presume using that would be staight forward. Is there any specific reason why you like to go for bio3d and R?
i want to perform SMD on namd for trimer protein. on which two chain should be stretched out in opposite direction on x axis. can you please help me regarding this?