hello everyone, I wonder if it is necessary to use “special_bonds” when simulating molecular chains with more than four atoms in LAMMPS (for example, CVFF for organic). If so, how to determine the weight coefficient of each interaction in “special_bonds”? Is it found in the force field file?
It may not be necessary to use the special_bonds command, because there is a default setting.
Yes, the required settings for exclusions need to be part of the force field documentation.