We optimize Tersoff potential to simulate the cBN better under nanoidentation.
Dislocations slip more easily along <110> and <112> directions on the {111} plane.
Shuffle-set dislocation slip along <112> direction on {111} plane first.
A tetrahedron structure is found in the initial stage of the indentation.