All geometry optimizations were performed with Gaussian 16 at the B3LYP/6-311G(d,p) level, although a selected species were also calculated with CAM-B3LYP. Aromaticity was evaluated in terms of electronic and magnetic indicators at the same level of theory. The electronic delocalization was evaluated by the electron density of delocalization bond (EDDBH), and magnetic properties were assessed using the nucleus-independent chemical shift (NICS), anisotropy of the induced current density (ACID) plots and gauge-including magnetically induced currents (GIMIC). Regarding EDDBH computations, NBO 3.1 and Multiwfn have been employed, where the former together with Avogadro 1.2 was used to obtaining EDDBH surfaces. NICS-XY scans were performed using the Aroma package and ACID plots were produced using AICD 2.0.0 program.