special_bonds command in combination with the lj/cut/coul/long pair potential

The cut off value 10 is in lj units and not in angstroms ?

Arun

Sorry for the confusion – my last answer was wrong. I am using the following initial settings, so the units are “real” and therefore the default skin cut off should be skin = 2.0 Angstroms (and the cut off value 10 is also in Angstroms).

units real
atom_style full
pair_style lj/cut 10.0
bond_style harmonic
angle_style harmonic
dihedral_style fourier