Nucleus-independent chemical shifts (NICS) values were calculated at the GIAO-cam-B3LYP/6-311+G(2d,p) level using geometry-optimized structures. Magnetically induced current densities were calculated using the GIMIC program based on the optimized structure. The NMR shielding constants were calculated using gaugeincluding atomic orbitals (GIAO) at cam-B3LYP/6-311+G(2d,p) level. The integration of the current density through a bond (nA T -1 ) was performed using a rectangular plane perpendicular to the plane of the macrocycle. The rectangular plane had dimensions of 10 bohr (5.2 Å) high 4.4 bohr (2.3 Å) wide and spanned the center of the bond of interest. The magnetic field was placed perpendicular to the plane of the macrocycle. Electronic indices of aromaticity (AV1245 and AVmin ) were calculated with the AIMAll and ESI-3D program packages.