The role of aromaticity in porphyrinoids is a current subject of debate due to the intricate structure of these macrocycles, which can adopt Hückel, Möbius and even figure-eight conformers. One of the main challenges in these large π-conjugated structures is identifying the most conjugated pathway because, among aromaticity descriptors, there are very few that can be applied coherently to this variety of conformers. In this paper, we have investigated the conjugated pathways in nine porphyrinoid compounds using several aromaticity descriptors, including BLA, BOA, FLU and HOMA, as well as the recently introduced AV1245 and AV indices. All the indices agree on the general features of these compounds, such as the fulfillment of Hückel's rule or which compounds should be more or less aromatic from the series. However, our results evince the difficulty of finding the most aromatic pathway in the macrocycle for large porphyrinoids. In fact, only AV is capable of recognizing the annulene pathway as the most aromatic one in the nine studied structures. Finally, we study the effect of the exchange in DFT functionals on the description of the aromaticity of the porphyrinoids. The amount of exact exchange quantitatively changes the picture for most aromaticity descriptors, AV being the only exception that shows the same qualitative results in all cases.

Two of the most popular rules to characterize the aromaticity of molecules are those due to Hückel and Baird, which govern the aromaticity of singlet and triplet states. In this work, we study how these rules fade away as the ring structure increases and an optimal overlap between p orbitals is no longer possible due to geometrical restrictions. To this end, we study the lowest-lying singlet and triplet states of neutral annulenes with an even number of carbon atoms between four and eighteen. First of all, we analyze these rules from the Hückel molecular orbital method and, afterwards, we perform a geometry optimization of the annulenes with several density functional approximations in order to analyze the effect that the distortions from planarity produce on the aromaticity of annulenes. Finally, we analyze the performance of three density functional approximations that employ different percentages of Hartree-Fock exchange (B3LYP, CAM-B3LYP and M06-2X) and Hartree-Fock. Our results reveal that functionals with a low percentage of Hartree-Fock exchange at long ranges suffer from severe delocalization errors that result in wrong geometrical structures and the overestimation of the aromatic character of annulenes.