Aromaticity of giant polycyclic aromatic hydrocarbons with hollow sites: super ring currents in super-rings

We present a systematic theoretical study based on semi-empirical, Hartree-Fock (HF), and density functional theory (DFT) models of a series of polycyclic aromatic hydrocarbons (PAHs) that exhibit hollow sites. In this study we focus particularly on the magnetic criteria of aromaticity, namely (1)H NMR and nucleus-independent chemical shifts (NICS), and on their relationships with other electronic properties. The computed shifts and NICS indices indicate that an external magnetic field induces exceptionally strong ring currents in even-layered PAH doughnuts, in particular in the layer directly adjacent to the central hole of double-layered compounds. These exceptionally strong ring currents also correlate with particularly small HOMO-LUMO gaps and electronic excitation energies and to abnormally high polarizabilities, indicating in turn that these compounds have a more pronounced metallic character. Comparison is made with further depictions of aromaticity in these systems and in [18]-[66]annulene rings by employing topological, structural, and energetic criteria.     

Aromaticity Survival in Hydrofullerenes: The Case of C66H4 with Its π-Aromatic Circuits

The isolated-pentagon rule (IPR) is a determining structural feature that accounts for hollow fullerene stabilization and properties related to C_n (n≥60) cages. The recent characterization of an unprecedented non-IPR hydrofullerene, C_2v C₆₆H₄, bearing two heptagons with adjacent fused-pentagon motifs, largely dismisses this feature. Herein, employing DFT calculations, the ¹³C NMR spectroscopy and aromatic behavior of C_2v C₆₆H₄ are explored. The results show the presence of three π-aromatic circuits at the bottom boat section of C₆₆H₄, indicating the unique features of this hydrofullerene in comparison to those of pristine C₆₀. In addition, under specific orientations of the external field, certain π-aromatic circuits are enabled, resulting in a more aromatic fullerene than that of C₆₀, but lower than that of the spherical aromatic C₆₀⁶⁻ fulleride. Notably, under a field aligned with the saturated carbon atoms, nonaromatic characteristics are exposed. This reveals that spherical-like cages can involve a complex magnetic response that heavily depends on the orientation of the applied field.     

The NICS‐XY‐Scan: Identification of Local and Global Ring Currents in Multi‐Ring Systems

Nucleus‐independent chemical shift (NICS)‐based methods are very popular for the determination of the induced magnetic field under an external magnetic field. These methods are used mostly (but not only) for the determination of the aromaticity and antiaromaticity of molecules and ions, both qualitatively and quantitatively. The ghost atom that serves as the NICS probe senses the induced magnetic field and reports it in the form of an NMR chemical shift. However, the source of the field cannot be determined by NICS. Thus, in a multi‐ring system that may contain more than one induced current circuit (and therefore more than one source of the induced magnetic field) the NICS value may represent the sum of many induced magnetic fields. This may lead to wrong assignments of the aromaticity (and antiaromaticity) of the systems under study. In this paper, we present a NICS‐based method for the determination of local and global ring currents in conjugated multi‐ring systems. The method involves placing the NICS probes along the X axis, and if needed, along the Y axis, at a constant height above the system under study. Following the change in the induced field along these axes allows the identification of global and local induced currents. The best NICS type to use for these scans is NICS_πZZ, but it is shown that at a height of 1.7 Å above the molecular plane, NICS_ZZ provides the same qualitative picture. This method, namely the NICS‐XY‐scan, gives information equivalent to that obtained through current density analysis methods, and in some cases, provides even more details.     

Correlation of delocalization indices and current-density maps in polycyclic aromatic hydrocarbons

Using multicentre delocalization indices, the ring current maps of a large set of polycyclic aromatic hydrocarbons (PAH) are reconstructed and compared with ab initio computations of the same maps in the pseudo-pi version of the ipsocentric approach to magnetic response. The quality of the comparison indicates that both delocalization and ring current approaches capture the same information about the aromatic nature of the PAH. Aromaticity as a global property, requires knowledge of more than single circuits, but the present results suggest no need to introduce a "multidimensional character" for aromaticity.     

On a possible growth mechanism for polycyclic aromatic hydrocarbon di-cations: C(7)H(6)(2+) + C(2)H(2)

The mechanism of the bond-forming reaction between C(7)H(6) (2+) and C(2)H(2) to yield C(9) entities has been investigated by density functional theory calculations with close comparison with experimental data. It is shown that the reaction produces the C(9)H(6) (2+) and C(9)H(7) (2+) di-cations with geometries most probably derived from the indene skeleton. In comparison, the formation of linear structures of di-cations is much more energy-demanding and therefore appears improbable.     

Ring currents as probes of the aromaticity of inorganic monocycles : P5-, As5-, S2N2, S3N3-, S4N3+, S4N42+, S5N5+, S42+ and Se42+

Current-density maps were calculated by the ipsocentric CTOCD-DZ/6-311G** (CTOCD-DZ=continuous transformation of origin of current density-diamagnetic zero) approach for three sets of inorganic monocycles: S(4) (2+), Se(4) (2+), S(2)N(2), P(5) (-) and As(5) (-) with 6 pi electrons; S(3)N(3) (-), S(4)N(3) (+) and S(4)N(4) (2+) with 10 pi electrons; and S(5)N(5) (+) with 14 pi electrons. Ipsocentric orbital analysis was used to partition the currents into contributions from small groups of active electrons and to interpret the contributions in terms of symmetry- and energy-based selection rules. All nine systems were found to support diatropic pi currents, reinforced by sigma circulations in P(5) (-), As(5) (-), S(3)N(3) (-), S(4)N(3) (+), S(4)N(4) (2+) and S(5)N(5) (+), but opposed by them in S(4) (2+), Se(4) (2+) and S(2)N(2). The opposition of pi and sigma effects in the four-membered rings is compatible with height profiles of calculated NICS (nucleus-independent chemical shifts).     

An insight into the local aromaticities of polycyclic aromatic hydrocarbons and fullerenes

In this work we quantify the local aromaticity of six-membered rings in a series of planar and bowl-shaped polycyclic aromatic hydrocarbons (PAHs) and fullerenes. The evaluation of local aromaticity has been carried out through the use of structurally (HOMA) and magnetically (NICS) based measures, as well as by the use of a new electronically based indicator of aromaticity, the para delocalization index (PDI), which is defined as the average of all the Bader delocalization indices between para-related carbon atoms in six-membered rings. The series of PAHs selected includes C(10)H(8), C(12)H(8), C(14)H(8), C(20)H(10), C(26)H(12), and C(30)H(12), with benzene and C(60) taken as references. The change in the local aromaticity of the six-membered rings on going from benzene to C(60) is analyzed. Finally, we also compare the aromaticity of C(60) with that of C(70), open [5,6]- and closed [6,6]-C(60)NH systems, and C(60)F(18).     

Electron delocalization in homoconjugated 7,7-diarylnorbornane systems: a computational and experimental study

A joint computational–experimental study has been carried out to analyze the homoconjugative interactions in 7,7‐diarylnorbornane (DPN) derivatives. The experimentally observed new bands in their UV/Vis have been accurately assigned by means of TD‐DFT calculations. Both experimental data and computations show that aromatic homoconjugation in acyclic systems is an effective mechanism for electron delocalization that resembles the situation described for polyphenylenes and polyenes. The effective homoconjugation length in homoconjugated oligomers is in the range of 6–7 aryl rings. The effect of substituents directly attached to the para carbon atom of the DPN moiety have been also studied. We found that the HOMO→LUMO vertical transitions can indeed be modified by the nature of the aromatic substituents in order to provoke dramatic changes in the electronic properties (i.e., in the absorption spectra) of the studied species.     

Understanding the Reactivity of Planar Polycyclic Aromatic Hydrocarbons: Towards the Graphene Limit

The Diels–Alder reactivity of maleic anhydride towards the bay regions of planar polycyclic aromatic hydrocarbons was explored computationally in the DFT framework. The process becomes more and more exothermic and the associated activation barriers become lower and lower when the size of the system increases. This enhanced reactivity follows an exponential behavior that reaches its maximum for systems having 18–20 benzenoid rings in their structures. This peculiar behavior was analyzed in detail by using the activation strain model of reactivity in combination with energy decomposition analysis. The influence of the change in the aromaticity of the polycyclic compound during the process on the respective activation barriers was also studied.     

Impact of C=C/B−N Replacement on the Diels–Alder Reactivity of Curved Polycyclic Aromatic Hydrocarbons

The influence of the replacement of C=C bonds by isoelectronic B−N moieties on the reactivity of π-curved polycyclic aromatic hydrocarbons has been computationally explored by means of density functional theory calculations. To this end, we selected the Diels–Alder cycloaddition reactions of the parent corannulene and its BN-doped counterparts with either cyclopentadiene or maleic anhydride. In addition, the analogous reactions involving larger buckybowls, such as BN-hemifullerene, BN-circumtrindene, and BN-fullerene, have been also considered. It has been found that whereas corannulene behaves as a dienophile, its BN counterpart better acts as a diene. In contrast, the larger BN-curved systems cannot be used as dienes in Diels–Alder reactions, but undergo facile (i.e., low barrier) cycloaddition reactions with cyclopentadiene. The observed trends in reactivity, which cannot be directly explained by using typical frontier molecular orbital arguments, are quantitatively described in detail by means of state-of-the-art computational methods, namely the activation strain model of reactivity combined with the energy decomposition analysis method. The results of our calculations highlight the crucial role of the curvature of the system on the reactivity and its influence on the strength of the orbital interactions between the deformed reactants during their transformations.