I am working on a system of particles which are bonded. I have 5000 bonds from one bonded interaction style. During my simulation, I need to remove a group of these bonds. How could I do that?
I have another question: Assume that you have a harmonic potential:
inter 10 harmonic $K $R
I want to change the R of a special bond during my simulation. If I define: inter 10 harmonic $K $R_new