Thank you to introduce me this potential, I implemented it, just one question :
In this potential you should determine two bonded neighbors ba1,........ . There is a little problem in the last site of every polymer. It just has one bonded neighbor.
May I ask how could I rectify this problem ? Thanks for your potential , its really helpful for me and simplifies my model.
I want to add a virtual site to the end of every tail, I dont know it is really helpful or not.