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.