This work revealed the spherical aromaticity of some inorganic E4 cages and their protonated E4H+ ions (E=N, P, As, Sb, and Bi). For this purpose, we employed several evaluations like (0D-1D) nucleus independent chemical shift (NICS), multidimensional (2D-3D) off-nucleus isotropic shielding σiso(r), and natural bond orbital (NBO) analysis. The magnetic calculations involved gauge-including atomic orbitals (GIAO) with two density functionals B3LYP and WB97XD, and basis sets of Jorge-ATZP, 6-311+G(d,p), and Lanl2DZp. The Jorge-ATZP basis set showed the best consistency. Our findings disclosed non-classical aromatic characters in the above molecules, which decreased from N to Bi cages. Also, the results showed more aromaticity in E4 than E4H+. The NBO analysis attributed the aromaticity in the above molecules to the residual density of the overlapping σ-bonding orbitals. So, the aromaticity in these molecules is unlike the classical aromaticity that is associated with electron delocalization. Scanning 1D σiso(r) variation along E−E bonds indicated a lowering in the shielding trend from N to Bi cages. The 3D results showed a similar decrease in the relative volumetric diffusion of the magnetic activity, whereas the volumetric ratio of V1ppm/V2ppm is almost constant for all the E4 cages.
