Hello,
I am trying to include constraints in my forcefield document which I need to run a mesoscale simulation. Let's say I want the bond length of A-B to be fixed at the equilibrium bond length, then I would have to constraint the bond force constant somehow. I tried to enter a very high value for the spring constant but the geometry optimization was terminated abruptly and it did not converge. Hence, the DPD run also failed as there was an unphysical integration step. Does anyone know how to work this out?
Regards,
Ricky