The “Anthracene Problem”: Closed-Form Conjugated-Circuit Models of Ring Currents in Linear Polyacenes

Effect of Benzo-Annelation on Triplet State Energies in Polycyclic Conjugated Hydrocarbons

In a series of earlier studies, the effect of benzo-annelation was found to be a useful tool for tuning the aromaticity in polycyclic conjugated compounds to desired level. In this work we studied the (anti)aromaticity of benzo-annelated derivatives of three conjugated hydrocarbons (anthracene, fluoranthene and biphenylene) in their lowest lying singlet (S0) and triplet (T1) states by means of the energy effect (ef), harmonic oscillator model of aromaticity (HOMA), multicentre delocalization indices (MCI), magnetically induced current densities (MICDs) and nucleus independent chemical shifts (NICS). We showed that benzo-annelation is a topology-based effect which can be used to modify the T1 state excitation energies (E(T1)). A quantitative model was established being able to accurately predict the E(T1) based only on the numbers of angularly, linearly and geminally annelated benzene rings. In addition, it was demonstrated that the E(T1) can be directly related to the (anti)aromatic character of the central ring in the studied molecules in their S0 state.

Aromaticity in the Electronic Ground and Lowest Triplet States of Molecules with Fused Thiophene Rings

Analysis of the variations of the off-nucleus isotropic magnetic shielding, σiso(r), around thiophene, thienothiophenes, dithienothiophenes and sulflowers in their electronic ground (S0) and lowest triplet (T1) states reveals that some of the features of aromaticity and bonding in these molecules do not fit in with predictions based on the popular Hückel's and Baird's rules. Despite having 4n π electrons, the S0 states of the sulflowers are shown to be aromatic, due to the local aromaticities of the individual thiophene rings. To reduce its T1 antiaromaticity, the geometry of thiophene changes considerably between S0 and T1: In addition to losing planarity, the carbon-carbon two 'double' and one 'single' bonds in S0 turn into two 'single' and one 'double' bonds in T1. Well-defined Baird-style aromaticity reversals are observed between the S0 and T1 states of only three of the twelve thiophene-based compounds investigated in this work, in contrast, the sulflower with six thiophene rings which is weakly aromatic in S0 becomes more aromatic in T1. The results suggest that the change in aromaticity between the S0 and T1 states in longer chains of fused rings is likely to affect mostly the central ring (or the pair of central rings); rings sufficiently far away from the central ring(s) can retain aromatic character.

Understanding the electronic pi-system of 2D covalent organic frameworks with Wannier functions

We investigate a family of hexagonal 2D covalent organic frameworks (COFs) with phenyl and biphenyl spacer units and different chemical linker species. Chemical trends are elucidated and attributed to microscopic properties of the [Formula: see text]-electron-system spanned by atomic [Formula: see text]-orbitals. We systematically investigate the electronic structure, delocalization of electronic states, effects of disorder, bond torsion, and doping, and correlate these with variable [Formula: see text]-conjugation and nucleus-independent chemical shift (NICS) aromaticity. Molecular orbitals are obtained from maximally localized Wannier functions that have [Formula: see text]- and [Formula: see text]-character, forming distinct [Formula: see text]- and [Formula: see text]-bands for all valence states. The Wannier-orbital description goes beyond simple tight-binding models and enables a detailed understanding of the electronic topology, effective electronic coupling and delocalization. It is shown that a meaningful comparison between COFs with different chemical elements can only be made by examining the entire [Formula: see text]-electron system, while a comparison of individual bands (e.g., bands near the Fermi energy) can be a insufficient to derive general design rules for linker and spacer monomer selection. We further identify delocalized states that are spread across tens or hundreds of pores of the 2D COFs and analyze their robustness against structural and energetic disorders like out-of-plane rotations of molecular fragments, different strength of energetic disorder and energetic shifts due to chemical doping.

Three-Dimensional Fully π-Conjugated Macrocycles: When 3D-Aromatic and When 2D-Aromatic-in-3D?

Several fully π-conjugated macrocycles with puckered or cage-type structures were recently found to exhibit aromatic character according to both experiments and computations. We examine their electronic structures and put them in relation to 3D-aromatic molecules (e.g., closo-boranes) and to 2D-aromatic polycyclic aromatic hydrocarbons. Using qualitative theory combined with quantum chemical calculations, we find that the macrocycles explored hitherto should be described as 2D-aromatic with three-dimensional molecular structures (abbr. 2D-aromatic-in-3D) and not as truly 3D-aromatic. 3D-aromatic molecules have highly symmetric structures (or nearly so), leading to (at least) triply degenerate molecular orbitals, and for tetrahedral or octahedral molecules, an aromatic closed-shell electronic structure with 6n + 2 electrons. Conversely, 2D-aromatic-in-3D structures exhibit aromaticity that results from the fulfillment of Hückel’s 4n + 2 rule for each macrocyclic path, yet their π-electron counts are coincidentally 6n + 2 numbers for macrocycles with three tethers of equal lengths. It is notable that 2D-aromatic-in-3D macrocyclic cages can be aromatic with tethers of different lengths, i.e., with π-electron counts different from 6n + 2, and they are related to naphthalene. Finally, we identify tetrahedral and cubic π-conjugated molecules that fulfill the 6n + 2 rule and exhibit significant electron delocalization. Yet, their properties resemble those of analogous compounds with electron counts that differ from 6n + 2. Thus, despite the fact that these molecules show substantial π-electron delocalization, they cannot be classified as true 3D-aromatics.

Calculated Magnetic Properties of a Recently Synthesized Corannulene-Dibenzobistetracene Hybrid

Comparison is made between two calculated magnetic properties of a recently synthesized 24-ring corannulene-dibenzobistetracene hybrid (Corannulene-DBT): nuclear independent chemical shifts published by Xu et al. and Hückel-London-Pople-McWeeny (HLPM) ring currents and bond currents that are presented here. There is agreement in 18 of the 24 rings in Corannulene-DBT. This is seen as further evidence that the HLPM approach has a surprising ability to reproduce complex patterns of current in large polycyclic hydrocarbons.

Partitioning Hückel–London Currents into Cycle Contributions

Ring-current maps give a direct pictorial representation of molecular aromaticity. They can be computed at levels ranging from empirical to full ab initio and DFT. For benzenoid hydrocarbons, Hückel–London (HL) theory gives a remarkably good qualitative picture of overall current patterns, and a useful basis for their interpretation. This paper describes an implemention of Aihara’s algorithm for computing HL currents for a benzenoid (for example) by partitioning total current into its constituent cycle currents. The Aihara approach can be used as an alternative way of calculating Hückel–London current maps, but more significantly as a tool for analysing other empirical models of induced current based on conjugated circuits. We outline an application where examination of cycle contributions to HL total current led to a simple graph-theoretical approach for cycle currents, which gives a better approximation to the HL currents for Kekulean benzenoids than any of the existing conjugated-circuit models, and unlike these models it also gives predictions of the HL currents in non-Kekulean benzenoids that are of similar quality.

Extension and Quantification of the Fries Rule and Its Connection to Aromaticity: Large-Scale Validation by Wave-Function-Based Resonance Analysis

The Fries rule is a simple, intuitive tool to predict the most dominant Kekulé structures of polycyclic aromatic hydrocarbons (PAHs), which is valuable for understanding the structure, stability, reactivity, and aromaticity of these conjugated compounds. However, it still remains an empirical hypothesis, with limited qualitative applications. Herein, we verify, generalize, and quantify the Fries rule based on the recently developed resonance analysis of the DFT wave functions of over 1500 PAH and fullerene molecules with over a billion Kekulé structures. The extended rules, counting the numbers of electrons within all rings (not just sextets), are able to rank the relative importance of all Kekulé structures for all considered systems. The statistically meaningful quantification also opens a way to evaluate ring aromaticity based on the resonance theory, which generally agrees well with conventional aromaticity descriptors. Furthermore, we propose a purely graph-based aromaticity indicator nicely applicable to PAHs and fullerenes, with no need of any quantum chemistry calculations, so that it can make valuable predictions for molecular properties that are related to local aromaticity.

Detailed Visualization of Aromaticity Using Isotropic Magnetic Shielding

For many years, Clar's aromatic sextet theory has served as a qualitative method for assessing the aromatic character of polycyclic aromatic hydrocarbons. A new approach, based on the calculation of isotropic magnetic shielding (IMS) contour plots, is shown to provide a feature-rich picture of aromaticity that is both quantitative yet still easily interpreted. Chemists are visual creatures who are adept at discerning reactivity and chemical behavior from molecular structures. To quote Roald Hoffmann, "People like pictures. Chemists live off them." Thus, the detailed image analysis we present simultaneously provides quantitative assessment of electronic structure, which is still easy-to-understand through visual inspection, embedded in an aesthetically appealing and intuitive picture that draws the reader in. We provide novel computed IMS contour plots for a representative selection of aromatic molecules. Where Clar's static drawings capture only a partial sketch of the electronic properties of a molecule, IMS contour plots present a detailed, global landscape of a molecule that sums all possible resonance structures. This novel analysis allows us to correct certain drawbacks of Clar's analysis with respect to polycyclic aromatics and quantitatively assess the bonding and electronic structure of acene hydrocarbons.

Ring-Current Maps for Benzenoids: Comparisons, Contradictions, and a Versatile Combinatorial Model

As a key diagnostic property of benzenoids and other polycyclic hydrocarbons, induced ring current has inspired diverse approaches for calculation, modeling, and interpretation. Grid-based methods include the ipsocentric ab initio calculation of current maps, and its surrogate, the pseudo-π model. Graph-based models include a family of conjugated-circuit (CC) models and the molecular-orbital Hückel-London (HL) model. To assess competing claims for physical relevance of derived current maps for benzenoids, a protocol for graph-reduction and comparison was devised. Graph reduction of pseudo-π grid maps highlights their overall similarity to HL maps, but also reveals systematic differences. These are ascribed to unavoidable pseudo-π proximity limitations for benzenoids with short nonbonded distances, and to poor continuity of pseudo-π current for classes of benzenoids with fixed bonds, where single-reference methods can be unreliable. Comparison between graph-based approaches shows that the published CC models all shadow HL maps reasonably well for most benzenoids (as judged by L1-, L2-, and L∞-error norms on scaled bond currents), though all exhibit physically implausible currents for systems with fixed bonds. These comparisons inspire a new combinatorial model (Model W) based on cycle decomposition of current, taking into account the two terms of lowest order that occur in the characteristic polynomial. This improves on all pure-CC models within their range of applicability, giving excellent adherence to HL maps for all Kekulean benzenoids, including those with fixed bonds (halving the rms discrepancy against scaled HL bond currents, from 11% in the best CC model, to 5% for the set of 18 360 Kekulean benzenoids on up to 10 hexagonal rings). Model W also has excellent performance for open-shell systems, where currents cannot be described at all by pure CC models (4% rms discrepancy against scaled HL bond currents for the 20112 non-Kekulean benzenoids on up to 10 hexagonal rings). Consideration of largest and next-to-largest matchings is a useful strategy for modeling and interpretation of currents in Kekulean and non-Kekulean benzenoids (nanographenes).

Interplay of Aromaticity and Antiaromaticity in N-Doped Nanographenes

The aromaticity of three nonplanar, fully conjugated aza-nanographenes built around a pyrrolo[3,2-b]pyrrole core is assessed through the application of two different computational procedures-GIMIC and NICS. We examine the calculated magnetically induced current densities (GIMIC) and nucleus-independent chemical shifts (NICS). The structural differences between these three apparently similar molecules lead to significantly different aromatic properties. GIMIC analysis indicates that the peripheral diatropic ring current of 3.9 nA/T for the studied bowl-shaped diaza-nanographene is the strongest, followed by the double [6]helicene which lacks seven-membered rings, and is practically nonexistent for the double [5]helicene possessing seven-membered rings. The biggest difference however is that in the two not-fully-fused molecules, the central pyrrole rings possess a significant diatropic current of about 4.1 nA/T, whereas there is no such current in the diaza-nanographene. Moreover, the antiaromaticity of the seven-membered rings is increasing while moving from double [5]helicene to diaza-nanographene (from -2.4 to -6.0 nA/T). The induced currents derived from NICS_π,zz-XY-scan analysis for all of the studied systems are in qualitative agreement with the GIMIC results. Subtle differences may originate from σ-electron currents in GIMIC or inaccuracy of NICS_π,zz values due to the nonplanarity of the systems, but the general picture is similar.

The pseudo-π model of the induced magnetic field: fast and accurate visualization of shielding and deshielding cones in planar conjugated hydrocarbons and spherical fullerenes

Hydrogen skeletal models accurately reproduce the π-induced magnetic field of planar PAHs and spherical fullerenes.

Aromaticity of acenes: the model of migrating π-circuits

In this work we extend the concept of migrating Clar's sextets to explain local aromaticity trends in linear acenes predicted by theoretical calculations and experimental data. To assess the link between resonance and reactivity and to rationalize the constant-height AFM image of pentacene we used the electron density of delocalized bonds and other functions of the one-electron density from conceptual density functional theory. The presented results provide evidence for migration of Clar's π-sextets and larger circuits in these systems, and clearly show that the link between the theoretical concept of aromaticity and the real electronic structure entails the separation of intra- and inter-ring resonance effects, which in the case of [n]acenes (n = 3, 4, 5) comes down to solving a system of simple linear equations.

The electron density of delocalized bonds (EDDB) applied for quantifying aromaticity

In this study the recently developed electron density of delocalized bonds (EDDB) is used to define a new measure of aromaticity in molecular rings. The relationships between bond-length alternation, electron delocalization and diatropicity of the induced ring current are investigated for a test set of representative molecular rings by means of correlation and principal component analyses involving the most popular aromaticity descriptors based on structural, electronic, and magnetic criteria. Additionally, a qualitative comparison is made between EDDB and the magnetically induced ring-current density maps from the ipsocentric approach for a series of linear acenes. Special emphasis is given to the comparative study of the description of cyclic delocalization of electrons in a wide range of organic aromatics in terms of the kekulean multicenter index KMCI and the newly proposed EDDB^k index.

The role of the long-range exchange corrections in the description of electron delocalization in aromatic species

In this article, we address the role of the long-range exchange corrections in description of the cyclic delocalization of electrons in aromatic systems at the density functional theory level. A test set of diversified monocyclic and polycyclic aromatics is used in benchmark calculations involving various exchange-correlation functionals. A special emphasis is given to the problem of local aromaticity in acenes, which has been a subject of long-standing debate in the literature. The presented results indicate that the noncorrected exchange-correlation functionals significantly overestimate cyclic delocalization of electrons in heteroaromatics and aromatic systems with fused rings, which in the case of acenes leads to conflicting local aromaticity predictions from different criteria. © 2017 Wiley Periodicals, Inc.

From C(58) to C(62) and back: Stability, structural similarity, and ring current

An increasing number of observations show that non-classical isomers may play an important role in the formation of fullerenes and their exo- and endo-derivatives. A quantum-mechanical study of all classical isomers of C₅₈ , C₆₀ , and C₆₂ , and all non-classical isomers with at most one square or heptagonal face, was carried out. Calculations at the B3LYP/6-31G* level show that the favored isomers of C₅₈ , C₆₀ , and C₆₂ have closely related structures and suggest plausible inter-conversion and growth pathways among low-energy isomers. Similarity of the favored structures is reinforced by comparison of calculated ring currents induced on faces of these polyhedral cages by radial external magnetic fields, implying patterns of magnetic response similar to those of the stable, isolated-pentagon C₆₀ molecule. © 2016 Wiley Periodicals, Inc.

Perimeter ring currents in benzenoids from Pauling bond orders

It is shown that the ring currents in perimeter hexagonal rings of Kekulean benzenoids, as estimated within the Randić conjugated-circuit model, can be calculated directly without tedious pairwise comparison of Kekulé structures or Kekulé counting for cycle-deleted subgraphs. Required are only the Pauling bond orders of perimeter bonds and the number of Kekulé structures of the benzenoid, both readily available from the adjacency matrix of the carbon skeleton. This approach provides easy calculation of complete current maps for benzenoids in which every face has at least one bond on the perimeter (as in the example of cata-condensed benzenoids), and allows qualitative evaluation of the main ring-current contributions to (1)H chemical shifts in general benzenoids. A combined Randić-Pauling model for correlation of ring current and bond length through bond order is derived and shown to be consistent with resilience of current under bond alternation.

Magnetic criteria of aromaticity

This review describes the current state of magnetic criteria of aromaticity. The introduction contains the fundamentals of ring currents in aromatic and antiaromatic systems, followed by a brief description of experimental and computational tools: NMR, diamagnetic susceptibility exaltation, current density analyses (CDA) and nucleus independent chemical shifts (NICS). This is followed by more comprehensive chapters: NMR - focusing on the work of R. Mitchell - NICS and CDA - describing the progress and development of the methods to their current state and presenting some examples of representative work.

Topological ring-currents and bond-currents in the altan-[r,s]-coronenes

π-Electron ring-currents for the altans of the four regular [r,s]-coronenes are calculated by the topological HLPM approach.

Writing with ring currents: selectively hydrogenated polycyclic aromatics as finite models of graphene and graphane

Alternating partial hydrogenation of the interior region of a polycyclic aromatic hydrocarbon gives a finite model system representing systems on the pathway from graphene to the graphane modification of the graphene sheet. Calculations at the DFT and coupled Hartree-Fock levels confirm that sp² cycles of bare carbon centres isolated by selective hydrogenation retain the essentially planar geometry and electron delocalization of the annulene that they mimic. Delocalization is diagnosed by the presence of ring currents, as detected by ipsocentric calculation and visualization of the current density induced in the π system by a perpendicular external magnetic field. These induced 'ring' currents have essentially the same sense, strength and orbital origin as in the free hydrocarbon. Subjected to the important experimental proviso of the need for atomic-scale control of hydrogenation, this finding predicts the possibility of writing single, multiple and concentric diatropic and/or paratropic ring currents on the graphene/graphane sheet. The implication is that pathways for free flow of ballistic current can be modelled in the same way.