I found that the weighting factors set by special_bonds ARE used by the pair style Tersoff, which compute many-body interactions.
For example, for a system that uses Tersoff for a structure with lists of bonds, angles, dihedrals (for debuging purpose, the bonded interactions are set to zero by settin K = 0.0), seting "special_bonds lj 0.0 0.0 0.5 coul 0.0 0.0 0.5" will apply the factor 0.0 to the Tersoff interactions; If setting ""special_bonds lj 1.0 1.0 1.0 coul 1.0 1.0 1.0", the Tersoff interactions are calculated as if for the same structure but without list of bonds, angles, and dihedrals.
Could you please check for that? Thank you!
I'm using lammps-11Aug17.