Current-Density Calculations on Zn-Porphyrin₄₀ Nanorings

Two porphyrinoid nanorings have been studied computationally. They were built by linking 40 Zn-porphyrin units with butadiyne bridges. The molecular structures belonging to the D₄₀ point group were fully optimized with the Turbomole program at the density functional theory (DFT) level using the B3LYP functional and the def2-SVP basis sets. The aromatic character was studied at the DFT level by calculating the magnetically induced current-density (MICD) susceptibility using the GIMIC program. The neutral molecules are globally non-aromatic with aromatic Zn-porphyrin units. Charged nanorings could not be studied because almost degenerate frontier orbitals led to vanishing optical gaps for the cations. Since DFT calculations of the MICD are computationally expensive, we also calculated the MICD using three pseudo-π models. Appropriate pseudo-π models were constructed by removing the outer hydrogen atoms and replacing all carbon and nitrogen atoms with hydrogen atoms. The central Zn atom was either replaced with a beryllium atom or with two inner hydrogen atoms. Calculations with the computationally inexpensive pseudo-π models yielded qualitatively the same magnetic response as obtained in the all-electron calculations.