Dear all I am simulating simple bead-rod polymer chains, so I apply constraints to all the bonds of my system using the p-lincs algorithm. I try to understand how do the constraint forces contribute to the pressure of the system, but I am a bit confused. Can anyone help? Regards, Tommy