Dear ESPResSo users,
I want to simulate a polymer with nearly constant length. The polymer experiences some large forces so I need almost rigid bonds to simulate it. The first idea is using a harmonic bond with large stiffness, and the second is using a combination of
Weeks-Chandler-Andersen potential and FENE with large stiffness. But when I was reading the user guide documentation I founded there is a rigid bond potential in Espresso. My question is about 3 parameters which this potential need. Can you tell me a little bit about this potential, parameters, and efficiency of that? Is this potential faster than those two initial ideas? Because the speed of simulation is important in my work. It would be my pleasure if you have any other suggestion?