I don’t know how to use the special_bonds command? If I do Langevin simulation of polymer with coarse-gained bead spring model and FENE bond potential. How can I give the special_bonds numeric coefficients ? What is the difference of special_bonds 0.0 1.0 0.0, special_bonds 0.0 1.0 1.0, and special_bonds 1.0 1.0 1.0? And, how can I choose them? Could some experts give me any suggestion?
Thanks in advance!