ABSTRACT:A theoretical investigation on the structure, stability, and aromaticity of M-subphthalocyanine (MOSubPc; MAB, Al, and Ga) was performed at the B3LYP/6-31ϩG*//B3LYP/6-31G* level. The comparison between MOSubPc and the corresponding MOphthalocyanine (MOPc) was considered. The geometry optimization of the MOSubPc shows that in the AlOSubPc and GaOSubPc, the steric repulsions among the three azacoupled isoindole moieties increase, as to their macrocycles tend to be far from planarity. The binding energies of ClOM . . . aza-coupled isoindole corrected by the basis set superposition error (BSSE), and the nucleus-independent chemical shift (NICS) values at the ring center, which are a simple and effective local aromaticity probe, were calculated. The results show that AlOSubPc is less stable than both BOSubPc and AlOPc for larger steric repulsion, smaller binding energy, and weaker aromaticity. In the same way, GaOSubPc is less stable than both BOSubPc and GaOPc. In addition, the ring expansion reactivity occurring in BOSubPc was confirmed by the global aromaticity mirrored by the electrophilicity index values. Therefore, the AlOSubPc and GaOSubPc remain unknown, while the corresponding compounds AlOPc and GaOPc are known experimentally.