The iso-chemical shielding surface (ICSS) analysis visualizes componential magnetic shielding tensor, which is essentially equivalent to extending the NICS analysis to three-dimensional space, providing more comprehensive information about molecular aromaticity. The isosurfaces of z-component of the ICSS, referred to as ICSSzz, as well as the color-filled contour map of ICSSzz in the slice plane parallel and perpendicular to the plane of benzene, carborazine, and borazine are presented in Figure carborazine and borazine are aromatic like benzene. Upon careful observation, it was found that benzene and carborazine have similar continuous peripheral deshielding isosurfaces at the same isovalue, while the deshielding region of the former is wider, indicating that although carborazine also shows strong aromaticity, its degree is slightly weaker than that of benzene. Prominently, the deshielding isosurfaces of borazine exhibit discontinuity near the H(B) atoms, which intuitively demonstrates that its aromaticity is significantly weaker than that of the carborazine. It is noteworthy that the borazine has a deshielding region wrapped in the shielding isosurface at its ring center, which is apparently due to the insufficient π electron conjugation over the molecular ring. The NICSzz curves in Figure indeed show that over a considerable distance near its center, the NICSzz value of borazine is positive, while the NICSzz of benzene and carborazine are negative at any position. The NICSzz values at the ring centers, referred to as NICS(0)zz, are calculated to be of -17.0, -9.61, and 11.9 ppm for benzene, carborazine, and borazine, respectively. Therefore, the NICS(0)zz erroneously determines borazine as an antiaromatic species. From Figure , it can be seen that NICS(0)zz is largely contributed by essentially localized  electrons of the systems, demonstrating that NICS(1)zz is much more reliable than NICS(0)zz in determining  aromaticity.