Molecular magnetic properties of heteroporphyrins: a theoretical analysis‡

B3LYP/6-31G(d) optimization of porphyrin, tetraphenylporphyrin and their 21,23-diheteroatom substituted derivatives with O, S, and Se heteroatoms was performed. A planar macrocycle was found in all cases except 21,23-dioxatetraphenylporphyrin which presents only slight deviations from planarity. A Bader analysis uncovers the presence of S–S and Se–Se interactions in the four corresponding heteroporphyrins, which appreciably distort the original unsubstituted macrocycles. In the minimum energy structures of heterotetraphenylporphyrins the four meso phenyl groups slant alternatively to right or left so that the two pairs of opposite phenyls present a staggered conformation. The π current induced by a perpendicular magnetic field in porphyrin bifurcates across both types of pyrrole subunits but the presence of O, S and Se heteroatoms in 21,23-diheteroporphyrins causes a diminution of the current density through the inner section of the two heterorings and, consequently, the current path goes then through the outer section of these rings. The NICS values at the ring critical points of the heterorings are much larger (in absolute value) than those at the pyrrole ring critical points but appreciably smaller than that at the ring critical point of a pyrrole molecule. In agreement with experimental data the 1H NMR present appreciable downfield shifts for the β H atoms of the heterorings in the 21,23-heterosubstituted molecules.