influenced by N-doping. As in the case of pristine one, the overall surface is predominantly negative for the NB7 model. The obtained magnitudes of surface minima reveal that the surface reactivities of the NC6 and NC8 doped BC3NTs are larger than those of NB7 doped one.This trend means that, surface reactivity of the NC6 and NC8 models of the zigzag BC3NT detect more changes with respect to the NB7 due to the nitrogen impurity.