Metallaaromatics have attracted considerable attention in recent years because they can display properties of both organic and organometallic species. However, it remains unclear whether Clar’s rule could be applied to organometallic chemistry despite its proposal in 1950s. Here, we investigate the relative stabilities of 49 organic and organometallic species by density functional theory (DFT) calculations.

The Clar structure of polybenzenoid hydrocarbons (PBHs) have attracted considerable interest of both theoretical and experimental chemists since it was proposed in the 1950s. However, it remains unclear whether the Clar structure could exist in inorganic PBHs, the boron nitride (BN) analogues where the alternate boron and nitrogen atoms are used to replace the carbon atoms of PBHs. Here, we carry out thorough density functional theory (DFT) calculations to probe the possibility of Clar structures in BN analogues of PBHs.

The degree of p-electron (de)localization and aromaticity of a series of polybenzenoid hydrocarbons (PBHs) has been analyzed through the π-contribution to the electron localization function (ELFπ), calculated at the B3LYP/6-311G(d,p) hybrid density functional theory level. The extent of p-electron delocalization in the various hexagons of a PBH was determined through analysis of the bifurcation values of the ELFp basins (BV(ELFp)), the spans in the bifurcation values in each hexagon (ΔBV(ELFπ)), and the ring-closure bifurcation values of the ELFπ (RCBV(ELFπ)).