I wanted to make sure I understand things correctly: in my simulations of non-interacting particles, I only use harmonic bonds, which at some point, need to be broken. I use the following commands to obtain this:
fix 1 all bond/break 10 1 2
bond_coeff 1 10.0 1
special_bonds lj/coul 0 1 1
delete_bonds all bond 0 remove
I am slightly confused as to whether this is the correct way of doing this because I do not want LJ and Coulomb interactions to have an impact on my harmonic potential. Could you please let me know whether this is the case here?
Thank you very much in advance for your contribution.
With best wishes,