Why is Benzene Unique? Screening Magnetic Properties of C6H6Isomers

On the aromaticity of actinide compounds

The chemistry of actinides has flourished since the late 2010s with the synthesis of new actinide complexes and clusters. On the theoretical side, a range of tools is available for the characterization of these heavy element-containing compounds, but discrepancies in the assessment of aromaticity using different tools have led to controversies. In this Perspective, we examine the origin of controversies relating to the aromaticity of metallic compounds, with a focus on actinides. The aromaticity of actinides is important, not because these molecules are numerous or have a special role in catalysis or reactivity, but because this topic pushes theories of aromaticity to their limits. Owing to its reference independence, the magnetic criterion of aromaticity has been the most popular choice for the characterization of the aromaticity of metallic compounds, including actinide compounds. Through examination of several case studies, we show why this criterion might be misleading for metallic species and explain how findings relating to actinide compounds could reshape theories of aromaticity, not just for actinides but perhaps also for well-known hydrocarbons.

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.

Is azulene's local aromaticity and relative stability driven by the Glidewell-Lloyd rule?

The local aromaticity of azulene has been studied to understand their electronic properties. For this purpose, we have used the magnetic criterion through magnetically induced current density maps, ring current strengths, NICSzz(1), and the bifurcation value of three-dimensional surfaces of NICSzz. On the other hand, the delocalization criterion was used by calculating the MCI and ELFπ. The results show that the five-membered ring (5-MR) is more aromatic than the seven-membered ring (7-MR) and more aromatic than the free C5H5- ring. The opposite case is seen for the seven-membered ring, which is less aromatic than the free C7H7+. The local aromatic rings in azulene are formed due to an intramolecular electron transfer from the 7-MR to the 5-MR. In addition, the proposed resonance structures that allow explaining the properties of azulene, such as the dipole moment or the relative stability (in comparison to other isomers), show a preference for the formation of 5-MRs; for this reason, it is possible to conclude that the aromaticity and relative stability of azulene is driven by the Glidewell-Lloyd rule.

A reinvestigation of the boron cluster B15+/0/-: a benchmark of density functionals and consideration of aromaticity models

This study presents a thorough reinvestigation of the B15+/0/- isomers, first employing coupled-cluster theory CCSD(T) calculations to validate the performance of different DFT functionals. The B15+ cation has two planar lowest-lying isomers, while the first 3D isomer is less stable than the global minimum by ∼10 kcal mol-1. The PBE functional, within this benchmark survey, has proved to be reliable in predicting relative energies for boron isomers. Other functionals such as the TPSSh, PBE0 and HSE06 result in good energy ordering of isomers but warrant reconsideration when distinguishing between 2D and 3D forms. Caution is needed for structures having high spin contamination, as it may lead to significant errors. The anomalously lower stability of the B15- anion with respect to its neighbours, in terms of electron detachment energy, was explained through a competition between both rectangle and disk models for its geometry. This elucidates its stability with 12 electrons in rectangle model and instability with 10 electrons in disk-shaped structure, emphasizing the value of employing such geometric models. The proximity of the σ* LUMO to the π HOMO also contributes to the weakening of the B15- stability.

Update for Isomerization Stabilization Energies: The Fulvenization Approach

An alternative approach for calculating aromatic stabilization energies is proposed based on transforming an (anti)aromatic ring into a fulvene isomer. This fulvenization process gives a value of 34.05 kcal·mol-1 for benzene in the singlet state and a value of -17.85 kcal·mol-1 in the triplet state. Additionally, it is possible to use experimental values (as long as they exist) for the calculation as the gas-phase formation enthalpies of benzene and fulvene, whose difference is 33.72 kcal·mol-1. On the other hand, this same approach has been evaluated on several six-membered rings, including those persubstituted, biradicals, azines, and inorganic analogues, giving results in agreement with those reported in the literature using different criteria. Additionally, it is possible to differentiate the aromaticity of the rings in polycyclic aromatic hydrocarbons according to Clar's rules. Assigning the (anti)aromatic character in various nonbenzenoid rings (neutral and charged), except for five- and seven-membered rings, is also possible. The construction of the fulvene isomers in PAHs is set such that nonaromaticity-related effects are not considered. The results show that the fulvenization approach is an effective and efficient approach that can serve as an alternative or complement to existing tools.

On the aromaticity and stability of benzynes in the ground and lowest-lying triplet excited states

In this work, we have revisited the aromaticity of benzyne isomers at the unrestricted density functional theory level (UDFT) using the energetic, magnetic, and delocalization criteria. In addition, this last criterion has also been analyzed employing complete active space (CASSCF) calculations. The results show conservation of aromaticity in these monocycles. Additionally it is observed that this trend is maintained in polycyclic aromatic hydrocarbon derivatives such as biradical didehydrophenanthrenes. Do these results imply a violation of Baird's rule? The answer is No, because this conservation in aromaticity is due to the loss of hydrogen atoms affects only the electronic σ skeleton and exerts a minor influence on the π cloud. Additionally, we have analyzed the relative stability of benzyne isomers and their relationship with experimental ΔES-T values. According to the literature, the stability of the benzynes in the singlet state is due to an effective interaction between the electrons of the biradical centers; however, this effect is completely reversed in the triplet state, which explains why the para isomer has the lowest ΔES-T gap.

Aromaticity: Quo Vadis

Aromaticity is one of the most deeply rooted concepts in chemistry. But why, if two-thirds of existing compounds can be classified as aromatic, is there no consensus on what aromaticity is? σ-, π-, δ-, spherical, Möbius, or all-metal aromaticity… why are so many attributes needed to specify a property? Is aromaticity a dubious concept? This perspective aims to reflect where the aromaticity community is and where it is going.

Aromaticity in P8 allotropes and (CH)8 analogues: significance of their 40 valence electrons?

The currently unknown phosphorus allotrope P₈ is of interest since its 40 total valence electrons is a "magic number" corresponding to a filled 1S²1P⁶1D¹⁰1S²1F¹⁴2P⁶ shell such as found in the relatively stable main group element clusters Al₁₃^- and Ge₉^4-. However, P₈ still remains as an elusive structure not realized experimentally. The lowest energy P₈ structure by a margin of ∼9 kcal mol^-1 is shown by density functional theory to be a cuneane analogue with no PP double bonds and two each of P₅, P₄, and P₃ rings. Higher energy P₈ structures are polycyclic systems having at most a single PP double bond. These P₈ systems are not "carbon copies" of the corresponding (CH)₈ hydrocarbons with exactly one hydrogen atom bonded to each carbon atom. Thus the lowest energy (CH)₈ structure is cyclooctatetraene with four CC bonds followed by benzocyclobutene with three CC bonds. The cuneane (CH)₈ structure is a relatively high energy isomer lying ∼36 kcal mol^-1 above cyclooctatetraene. The cubane P₈ and (CH)₈ structures are even higher energy structures, lying ∼37 and ∼74 kcal mol^-1 in energy above the corresponding global minima. Our results demonstrate differences in medium sized aggregates of elemental phosphorus and isolobal hydrocarbon species.

α-Rhenabenzofuran with nonaromatic T0 and aromatic S1 states

The first α-rhenabenzofuran complexes 2-5 are obtained from one-pot reactions of ReCl₃(PMePh₂)₃ with o-ethynyl phenols. With a delocalized structure, these paramagnetic compounds are nonaromatic at the ground state, which is a triplet. But they exhibit an aromatic singlet excited state, revealed by DFT studies.

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.

Fluorene analogues of xanthenes - low molecular weight near-infrared dyes

Fluorene-based analogues of fluorescein, rhodol, and rhodamine exhibit absorption and fluorescence beyond 800-900 nm in water, 300-400 nm red-shifted compared to the original oxygen-bridged xanthene dyes, giving potential access to low molecular weight fluorescent markers for the second biological window (NIR-II, ca. 1000-1350 nm).

Synthesis, structure and aromaticity of metallapyridinium complexes

The first rhena-analogues of pyridinium (cyclopropametalla-2-isoquinolinium complexes) are obtained from o-ethynyl benzonitriles. Structural analysis and DFT calculations confirm their aromatic nature. Compared to rhenapyrylium, rhenapyridinium has a slightly stronger Hückel π-aromaticity, while a chlorine substituent on the rhenapyridinium ring decreases its aromaticity, which is revealed by NICS, EDDB, MCI and ΔBV(ELF_π) analysis.

A Novel Automated Screening Method for Combinatorially Generated Small Molecules

A main challenge in the enumeration of small-molecule chemical spaces for drug design is to quickly and accurately differentiate between possible and impossible molecules. Current approaches for screening enumerated molecules (e.g., 2D heuristics and 3D force fields) have not been able to achieve a balance between accuracy and speed. We have developed a new automated approach for fast and high-quality screening of small molecules, with the following steps: (1) for each molecule in the set, an ensemble of 2D descriptors as feature encoding is computed; (2) on a random small subset, classification (feasible/infeasible) targets via a 3D-based approach are generated; (3) a classification dataset with the computed features and targets is formed and a machine learning model for predicting the 3D approach's decisions is trained; and (4) the trained model is used to screen the remainder of the enumerated set. Our approach is ≈8× (7.96× to 8.84×) faster than screening via 3D simulations without significantly sacrificing accuracy; while compared to 2D-based pruning rules, this approach is more accurate, with better coverage of known feasible molecules. Once the topological features and 3D conformer evaluation methods are established, the process can be fully automated, without any additional chemistry expertise.

Bonding and Aromaticity in Electron-Rich Boron and Aluminum Clusters

In this work bonding and aromaticity of triply bonded atoms of group 13 elements (M≡M, M = B and Al) in recently characterized B₂Al₃^-, Na₃Al₂^-, and Na₄Al₂ are studied. Here, I show that although molecular orbital-based analyses characterize triple bonds, the electropositive nature of group 13 elements gives these bonds unique characteristics. The bond orders derived from the delocalization index, topology of the electron density, and local characteristics of (3, -1) critical points, as defined within the context of quantum theory of atoms in molecules, do not conform with those of ordinary triple bonds. In Na₃Al₂^- and Na₄Al₂ clusters non-nuclear attractors form between the electropositive Al atoms acting like pseudo atoms. The bond between boron atoms in B₂Al₃^- is more similar to an ordinary triple covalent bond benefiting from the exchange-correlation component of the interatomic interaction energy as defined via interacting quantum atom theory. However, extreme electrostatic repulsion between negatively charged boron atoms attenuates this bond. Finally, current density analysis suggests that B₂Al₃^- is a magnetic aromatic system, nearly 50% more aromatic compared to benzene.

The symmetry principle of antiaromaticity

Current definitions of aromaticity are purely phenomenological and relate symmetry, reactive stability and the occurrence of molecular diamagnetic response currents. The antithetical concept of antiaromaticity provides a connection between the contrary properties: structural instability or distortion out of higher symmetry, a small HOMO-LUMO gap, and paramagnetic response currents. We reveal the symmetry principle that is underlying antiaromaticity by showing an intimate and strict symmetry induced relation between these properties. The principle is mathematically rigorous and can be formulated like: First order (and related) Jahn-Teller distorted molecules out of non-cubic and non-icosahedral point groups are prone to induced paramagnetism in magnetic fields parallel to the main axis of symmetry. We show by the exemplary cases of cyclobutadiene, cyclcooctatetraene, pentalene and manganese trifluoride how this principle works and discuss this new perspective on antiaromaticity.

Relativity or aromaticity? A first-principles perspective of chemical shifts in osmabenzene and osmapentalene derivatives

The topology of the magnetically induced current density in osmabenzene suggests that the molecule is a novel type of Craig–Möbius aromatic system.

Local aromaticity mapping in the vicinity of planar and nonplanar molecules

We report on nucleus-independent magnetic shielding (NICS) scans over the centers of six- and five-membered rings in selected metal phthalocyanines (MPc) and fullerene C60 for more accurate characterization of local aromaticity in these compounds. Detailed tests were conducted on model aromatic molecules including benzene, pyrrole, indole, isoindole, and carbazole and subsequently applied to H₂ Pc, ZnPc, Al(OH)Pc, and CuPc. Similar behavior of three selected magnetic probes, Bq, ³ He, and ⁷ Li⁺ , approaching perpendicularly the ring centers, was observed. For better visualization of shielding zone over the centers of aromatic rings, we introduced a simple mathematical procedure: the first and second derivatives of scan curves with respect to magnetic probe position enabled their additional examination. It allowed an easier localization of curve minimum and discrimination between areas in space varying by the magnetic field magnitude and to illustrate local aromaticity of two different kinds of rings in MPc with better resolution. Our results supported earlier reports on very low aromaticity indexes of pyrrole ring incorporated into MPc and significant aromaticity of the central macrocycle. No direct dependence between harmonic oscillator model of aromaticity and NICS was observed. Instead, a correlation between position of scan curve minimum and its magnitude were observed. In addition, the NICS values and ³ He chemical shifts in the middle of neutral C60 and C60^6- anion agreed well with the reported experimental NMR values for He@C60 and He@C60^6- .

Delocalization energy retrieved from the current density tensor

The anisotropy of the magnetically induced current density tensor can be computed by its original formulation, known in the literature as ACID, or by its revised definition, AACID, which takes into account the asymmetric nature of the tensor. In polycyclic aromatic hydrocarbons, the excess of the integrated value of AACID over that computed for isolated ethylenes correlate with the Hückel delocalization energy.

Magnetic Diversity in Heteroisocorroles: Aromatic Pathways in 10-Heteroatom-Substituted Isocorroles

A recent study on magnetically induced currents in 10-isocorrole derivatives indicated that both the free-base and metal-complexed forms of the unsubstituted macrocycle are homoaromatic. Furthermore, depending on the substituents at the 10-position, the aromatic character was found to swing between substantially homoaromatic to substantially antihomoaromatic. Heteroisocorroles, in which the saturated 10-position has been replaced by a heteroatom-containing group X, are predicted to exhibit even more dramatic variations in aromatic character, ranging from strongly aromatic (X = O, NH, PH, and S) to strongly antiaromatic (X = BH and CO). Interestingly, the experimentally studied X = SiMe2 case does not appear to sustain a significant global ring current.

The electronic structure and stability of germanium tubes Ge30H12 and Ge33H12

To explore the possibility of large cylindrical Ge clusters, we investigated the molecular and electronic structure of the germanium tube Ge₃₀H₁₂, composed of six parallel, planar hexagons using DFT calculations, and Ge₃₃H₁₂ generated from the insertion of Ge atoms at the center of three inner hexagons of Ge₃₀H₁₂.