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

Substituent effects and electron delocalization in five-membered N-heterocycles

Five-membered N-heterocycles are principal constituents of many compounds of vital importance in various fields of chemistry, biochemistry or pharmaceutical chemistry. For this reason, unequivocal identification of structural factors determining electron donating/withdrawing properties of specific groups attached to the heterocyclic moiety becomes an utmost need together with elucidation of the substitution-induced changes in cyclic and noncyclic electron delocalization. Thus, quantum-chemical calculations were performed for pyrrole, imidazole, pyrazole, 1,2,3- and 1,2,4-triazole, and their C-substituted mono-derivatives (X = NO2, CN, Br, Cl, F, SH, OH, NH2). The obtained dataset contains information on substituent properties (cSAR - charge of the substituent active region method), delocalization (EDDB - electron density of delocalized bonds) and geometry. It follows that the positions of endocyclic N atoms relative to the substituent influence in the most profound manner its properties. N atoms in ortho positions significantly boost the electron-donation and weaken the electron-withdrawal by induction. Another factor is the resonance charge transfer from the substituents to N atoms, and then inductive interactions with further (non-ortho) N atoms. While substituent constants correctly describe the changes of their properties (including those attached to the heterocycles), a testimony to Hammett's genius, quantum chemical models must be used to quantify the exact properties. In most heterocycles, electron-donating substituents hinder the cyclic delocalization, except 4-pyrazole. The applied recent EDDB method allows to study this phenomenon in detail. It follows that changes in aromaticity originate from the π-electronic effects of substituents on the ring bonds, changing the localization and delocalization of particular bonds in a correlated manner.

Revisiting Pentazole: An Investigation into the Intriguing Molecule Exhibiting Dual Organic and Inorganic Characteristics

Pentazole (cyclo-HN5) is a unique heterocycle categorized as both an organic and inorganic compound. However, attempts to synthesize and characterize cyclo-HN5 have been unsuccessful thus far. In this study, we synthesized a cyclo-HN5 solution and investigated the spectra, structure, aromaticity, acidity, and stability of cyclo-HN5. The lone pair of electrons on the protonated N atom of cyclo-HN5 participates in π-electron delocalization, forming two N═N bonds. Further investigations suggest that cyclo-HN5 exhibits significantly decreased π aromaticity and slightly lower σ aromaticity than cyclo-N5-. Experimental results suggest that pure cyclo-HN5 is unstable at ambient temperatures and pressures, but it can be isolated at high pressures or stabilized in solution by abundant hydrogen bonds. The pKa of cyclo-HN5 was determined as 1.63 (H2O, 25 °C) via potentiometric titration, indicating that cyclo-HN5 is a medium-strong acid. This study reveals the fundamental structure and properties of cyclo-HN5, thereby providing important data for advancing cyclo-HN5 chemistry.

Aromaticity in Semi-Condensed Figure-Eight Molecules

Electron delocalization and aromaticity was comparatively evaluated in recently synthesized figure-eight molecules made of two condensed U-shaped polycyclic aromatic hydrocarbon moieties connected either by two single bonds or by two para-phenylene groups. The selected examples include molecules that incorporate eight-membered and sixteen-membered rings, as well as a doubly [5]helicene-bridged (1,4)cyclophane. We probe whether some electron delocalization could occur through the stereogenic single bonds in these molecules: Is aromaticity purely (semi-)local, or possibly also global in these molecules? It was concluded that the situation can go from a purely (semi-)local character when the dihedral angle at the connecting single bonds is large, such as in biphenyl, to a predominantly (semi-)local character with a minor global contribution when the dihedral angle is small, such as in the para-phenylene connectors of the [5] helicene-bridged cyclophane.

Formation and reactivity of a unique M⋯C-H interaction stabilized by carborane cages

Broadening carborane applications has consistently been the goal of chemists in this field. Herein, compared to alkyl or aryl groups, a carborane cage demonstrates an advantage in stabilizing a unique bonding interaction: M⋯C-H interaction. Experimental results and theoretical calculations have revealed the characteristic of this two-center, two-electron bonding interaction, in which the carbon atom in the arene ring provides two electrons to the metal center. The reduced aromaticity of the benzene moiety, long distance between the metal and carbon atom in arene, and the upfield shift of the signal of M⋯C-H in the nuclear magnetic resonance spectrum distinguished this interaction from metal⋯C π interaction and metal-C(H) σ bonds. Control experiments demonstrate the unique electronic effects of carborane in stabilizing the M⋯C-H bonding interaction in organometallic chemistry. Furthermore, the M⋯C-H interaction can convert into C-H bond metallization under acidic conditions or via treatment with t-butyl isocyanide. These findings deepen our understanding regarding the interactions between metal centers and carbon atoms and provide new opportunities for the use of carboranes.

Insight into the Stability of Pentazolyl Derivatives based on Covalent Bond

Pentazole is regarded as a unique inorganic molecule that possess organic heterocyclic structure. Therefore, the research on pentazolyl derivatives represents a cutting-edge direction in both contemporary inorganic chemistry and heterocyclic chemistry. Moreover, their synthesis is regarded as the most significant research topic in the field of energetic materials due to the great potential of pentazolyl derivatives to breakthrough the energy bottleneck of CHNO-based energetic materials. However, synthesizing pentazolyl derivatives is challenging. To provide a theoretical support for the synthesis, we conducted theoretical studies on six single-ring pentazolyl derivatives with different functional groups. The results suggest that derivatization reduces the bond strength and weakens the aromaticity of the pentazolate ring. Further analysis showed that derivatization mainly affects the π aromaticity of the pentazolate ring, and ultimately causing poor stability of the pentazolyl derivatives. Among the six derivatives investigated in this study, fluoro pentazole (cyclo-N5-F) and hydroxyl pentazole (cyclo-N5-OH) possess good aromaticity, which is similar to the reported cyclo-N5-NCHN(CH3)2. Further calculations show that the kinetic stability of cyclo-N5-OH is higher than that of cyclo-N5-F. These results collectively indicate that cyclo-N5-OH is a promising candidate for synthesizing single-ring pentazolyl derivatives.

The n,π* States of Heteroaromatics: When are They the Lowest Excited States and in What Way Can They Be Aromatic or Antiaromatic?

Heteroaromatic molecules are found in areas ranging from biochemistry to photovoltaics. We analyze the n,π* excited states of 6π-electron heteroaromatics with in-plane lone pairs (nσ, herein n) and use qualitative theory and quantum chemical computations, starting at Mandado's 2n + 1 rule for aromaticity of separate spins. After excitation of an electron from n to π*, a (4n + 2)π-electron species has 2n + 2 πα-electrons and 2n + 1 πβ-electrons (or vice versa) and becomes πα-antiaromatic and πβ-aromatic. Yet, the antiaromatic πα- and aromatic πβ-components seldom cancel, leading to residuals with aromatic or antiaromatic character. We explore vertically excited triplet n,π* states (3n,π*), which are most readily analyzed, but also singlet n,π* states (1n,π*), and explain which compounds have n,π* states with aromatic residuals as their lowest excited states (e.g., pyrazine and the phenyl anion). If the πβ-electron population becomes more (less) uniformly distributed upon excitation, the system will have an (anti)aromatic residual. Among isomers, the one that has the most aromatic residual in 3n,π* is often of the lowest energy in this state. Five-membered ring heteroaromatics with one or two N, O, and/or S atoms never have n,π* states as their first excited states (T1 and S1), while this is nearly always the case for six-membered ring heteroaromatics with electropositive heteroatoms and/or highly symmetric (D2h) diheteroaromatics. For the complete compound set, there is a modest correlation between the (anti)aromatic character of the n,π* state and the energy gap between the lowest n,π* and π,π* states (R2 = 0.42), while it is stronger for monosubstituted pyrazines (R2 = 0.84).

Helicene Aromaticity Deviates from the Clar Rule-On the Electronic Dissimilarity of Large Isomeric Fibonacenes

This combined experimental and theoretical study illustrates the profound consequences of non-planarity on the electronic properties of polycyclic arenes. Three isomeric [10]fibonacene tetraesters were synthesized through a robust and regiocontrolled Perkin/Mallory approach: a nearly planar [10]phenacene derivative, a moderately twisted [10]semicircle derivative, and a 3D non-planar [10]helicene derivative. The photophysical properties of the 3D [10]helicene isomer were found to be dramatically different from the comparable ones of the [10]phenacene and [10]semicircle isomers. The aromatic properties of the [10]phenacene and [10]semicircle isomers conform well with their predictive Kekulé and Clar analyses, but the [10]helicene isomer deviates from these general topological rules, which appears to be a general phenomenon for [n]fibonacenes with n≥9.

[K2(Bi@Pd12@Bi20)]4-: An Endohedral Inorganic Fullerene with Spherical Aromaticity

Inorganic fullerene clusters have attracted widespread attention due to their highly symmetrical geometric structures and intrinsic electronic properties. However, cage-like clusters composed of heavy metal elements with high symmetry are rarely reported, and their synthesis is also highly challenging. In this study, we present the synthesis of a [K2(Bi@Pd12@Bi20)]4- cluster that incorporates a {Bi20} cage with pseudo-Ih symmetry, making it the largest main group metal cluster compound composed of the bismuth element to date. Magnetic characterization and theoretical calculations suggest that the spin state of the overall cluster is a quartet. Quantum chemical calculations reveal that the [Bi20]3- cluster has a similar electronic configuration to C606- and the [Bi@Pd12@Bi20]6- cluster exhibits a unique open-shell aromatic character.

3,4-Dimethylenecyclobutene: A Building Block for Design of Macrocycles with Excited State Aromatic Low-Lying High-Spin States

3,4-Dimethylenecyclobutene (DMCB) is an unusual isomer of benzene. Motivated by recent synthetic progress to substituted derivatives of this scaffold, we carried out a theoretical and computational analysis with a particular focus on the extent of (anti)aromatic character in the lowest excited states of different multiplicities. We found that the parent DMCB is non-aromatic in its singlet ground state (S0), lowest triplet state (T1), and lowest singlet excited state (S1), while it is aromatic in its lowest quintet state (Q1) as this state is represented by a triplet multiplicity cyclobutadiene (CBD) ring and two uncoupled same-spin methylene radicals. Interestingly, the Q1 state, despite having four unpaired electrons, is placed merely 4.8 eV above S0, and there is a corresponding singlet tetraradical 0.16 eV above. The DMCB is potentially a highly useful structural motif for the design of larger molecular entities with interesting optoelectronic properties. Here, we designed macrocycles composed of fused DMCB units, and according to our computations, two of these have low-lying nonet states (i. e., octaradical states) at energies merely 2.40 and 0.37 eV above their S0 states as a result of local Hückel- and Baird-aromatic character of individual 6π- and 4π-electron monocycles.

N2 cleavage by silylene and formation of H2Si(μ-N)2SiH2

Fixation and functionalisation of N2 by main-group elements has remained scarce. Herein, we report a fixation and cleavage of the N ≡ N triple bond achieved in a dinitrogen (N2) matrix by the reaction of hydrogen and laser-ablated silicon atoms. The four-membered heterocycle H2Si(μ-N)2SiH2, the H2SiNN(H2) and HNSiNH complexes are characterized by infrared spectroscopy in conjunction with quantum-chemical calculations. The synergistic interaction of the two SiH2 moieties with N2 results in the formation of final product H2Si(μ-N)2SiH2, and theoretical calculations reveal the donation of electron density of Si to π* antibonding orbitals and the removal of electron density from the π bonding orbitals of N2, leading to cleave the non-polar and strong NN triple bond.

A Janus carbaporphyrin pseudo-dimer

Carbaporphyrin dimers, investigated for their distinctive electronic structures and exceptional properties, have predominantly consisted of systems containing identical subunits. This study addresses the associated knowledge gap by focusing on asymmetric carbaporphyrin dimers with Janus-like characteristics. The synthesis of a Janus-type carbaporphyrin pseudo-dimer 5 is presented. It displays antiaromatic characteristics on the fused side and nonaromatic behavior on the unfused side. A newly synthesized tetraphenylene (TPE) linked bis-dibenzihomoporphyrin 8 and a previously reported dibenzo[g,p]chrysene (DBC) linked bis-dicarbacorrole 9 were prepared as controls. Comprehensive analyses, including 1H NMR spectral studies, single crystal X-ray diffraction analyses, and DFT calculations, validate the mixed character of 5. A further feature of the Janus pseudo-dimer 5 is that it may be transformed into a heterometallic complex, with one side coordinating a Cu(III) center and the other stabilizing a BODIPY complex. This disparate regiochemical reactivity underscores the potential of carbaporphyrin dimers as versatile frameworks, with electronic features and site-specific coordination chemistry controlled through asymmetry. These findings position carbaporphyrin dimers as promising candidates for advances in electronic structure studies, coordination chemistry, materials science, and beyond.

Central Ring Puckering Enhances the Stokes Shift of Xanthene Dyes

Small-molecule dyes are generally designed based on well-understood electronic effects. However, steric hindrance can promote excited-state geometric relaxation, increasing the difference between the positions of absorption and emission bands (the Stokes shift). Accordingly, we hypothesized that sterically induced central ring puckering in xanthene dyes could be used to systematically increase their Stokes shift. Through a combined experimental/quantum-chemical approach, we screened a group of (9-acylimino)-pyronin dyes with a perturbed central ring geometry. Our results showed that an atom with sp3 hybridization in position 10 of (9-acylimino)-pyronins induces central ring puckering and facilitates excited-state geometric relaxation, thereby markedly enhancing their Stokes shifts (by up to ~2000 cm-1). Thus, we prepared fluorescent (9-acylimino)-pyronin pH sensors, which showed a Stokes shift disparity between acid and base forms of up to ~8700 cm-1. Moreover, the concept of ring puckering-enhanced Stokes shift can be applied to a wide range of xanthene analogues found in the literature. Therefore, central ring puckering may be reliably used as a strategy for enhancing Stokes shifts in the rational design of dyes.

In Silico Analysis of the Aromaticity of Some Carbo-Metallabenzenes and Carbo-Dimetallabenzenes (Carbo-mers Proposed from Metallabenzenes)

In the current work, we introduce a novel class of molecules termed carbo-metallabenzenes, and their aromaticity has been comprehensively analyzed. The molecules were strategically designed with the insertion of acetylene (C≡C or C2) units in some selected metallabenzenes. Furthermore, if a second metallic unit is inserted (replacing a sp2 carbon) in the carbo-metallabenzenes rings, a new family of carbo-mers is generated, and this second group has been named as carbo-dimetallabenzenes. The primary objective of this work is to ascertain, through various methodologies, whether these newly proposed molecules retain the aromatic characteristics observed in carbo-benzene. The methodologies employed for bond analysis and aromaticity exploration include the analysis of the molecular orbitals, energy decomposition analysis, electron density of delocalized bonds, magnetically induced current density, and the induced magnetic field (Bind). This study sheds light on that the insertion of the metallic centers reduces the electronic delocalization and their aromaticity is, in some cases, comparable with the electronic delocalization of the inorganic iminobora-borazine and also provides valuable insights into their electronic structure through a multifaceted analysis.

Understanding Aromaticity in [5]Helicene-Bridged Cyclophanes: A Comprehensive Study

This study explores the aromaticity of doubly [5]helicene-bridged (1,4)cyclophane and triply [5]helicene-bridged (1,3,5)cyclophane via calculations of the magnetic response and of electronic aromaticity indices. The primary objective is to assess the π-electron delocalization to determine whether they sustain global ring currents associated with π aromaticity. The molecules show local ring currents in the presence of an external magnetic field. The ring currents flow diatropically in the stacked six-membered rings and in the helicene arms. However, these π currents are not interconnected due to the discontinuity of the π delocalization at the C-C single bonds connecting the central six-membered rings to the helicene arms. Electronic indices suggest that the helicene-arm systems have significantly smaller electron delocalization than benzene. The reduction in the delocalization does not compromise their ability to exhibit ring currents in the presence of an external magnetic field. The analysis provides further evidence that the magnetic criteria yield a different degree of aromaticity for the helicene arms than obtained in the calculation of the electronic aromaticity indices. However, both approaches confirm that the studied molecules are not globally aromatic.

Craig-Hückel Hybrid Aromatic Metalla-dehydro[11]annulenes Constructed by a Formal [10+1] Cycloaddition Reaction

Aromatic metalla-annulenes are important aromatic compounds, research into which has been mainly concentrated on metal-benzenes and their lower homologues. Reports on their superior homologs are rare, and this has greatly limited the systematic study of their properties. In this work, a series of osma-dehydro[11]annulenes with good air and thermal stability were prepared in high yields through a simple [10+1] strategy, by incorporating a metal fragment into conjugated ten-carbon chains in a one-pot reaction. They are the first monometallic aromatic metalla-[n]annulenes with the ring size larger than 6, and their Craig-Hückel hybrid aromaticity is supported by various physical and computational parameters. Besides, these complexes show versatile reactivities, not only giving further evidence for their aromaticity, but also demonstrating their physical and chemical properties can easily be regulated. This work enriches the metalla-aromatic chemistry, and provides a new avenue for the synthesis of large metalla-annulenes with different ring sizes.

%VBur index and steric maps: from predictive catalysis to machine learning

Steric indices are parameters used in chemistry to describe the spatial arrangement of atoms or groups of atoms in molecules. They are important in determining the reactivity, stability, and physical properties of chemical compounds. One commonly used steric index is the steric hindrance, which refers to the obstruction or hindrance of movement in a molecule caused by bulky substituents or functional groups. Steric hindrance can affect the reactivity of a molecule by altering the accessibility of its reactive sites and influencing the geometry of its transition states. Notably, the Tolman cone angle and %VBur are prominent among these indices. Actually, steric effects can also be described using the concept of steric bulk, which refers to the space occupied by a molecule or functional group. Steric bulk can affect the solubility, melting point, boiling point, and viscosity of a substance. Even though electronic indices are more widely used, they have certain drawbacks that might shift preferences towards others. They present a higher computational cost, and often, the weight of electronics in correlation with chemical properties, e.g. binding energies, falls short in comparison to %VBur. However, it is worth noting that this may be because the steric index inherently captures part of the electronic content. Overall, steric indices play an important role in understanding the behaviour of chemical compounds and can be used to predict their reactivity, stability, and physical properties. Predictive chemistry is an approach to chemical research that uses computational methods to anticipate the properties and behaviour of these compounds and reactions, facilitating the design of new compounds and reactivities. Within this domain, predictive catalysis specifically targets the prediction of the performance and behaviour of catalysts. Ultimately, the goal is to identify new catalysts with optimal properties, leading to chemical processes that are both more efficient and sustainable. In this framework, %VBur can be a key metric for deepening our understanding of catalysis, emphasizing predictive catalysis and sustainability. Those latter concepts are needed to direct our efforts toward identifying the optimal catalyst for any reaction, minimizing waste, and reducing experimental efforts while maximizing the efficacy of the computational methods.

Unraveling the Push-Pull Effect in Acenes, Polyenes and Polyynes

Substituent effects (SEs) are fundamental for predicting molecular reactivity, while polyene, polyyne and acene derivatives are precursors to compounds with diverse applications. Computations were performed for Y-R-X systems, where reaction sites Y=NO2 and O- , substituents X=NO2 , CN, Cl, H, OH, NH2 , and spacers R=polyene, polyyne (n=1-5, 10 repeating units) and acene (up to tetracene). The cSAR (charge of the substituent active region) approach allowed to present, for the first time, quantitative relations describing the spacer's electron-donating and withdrawing properties as a function of n and the spacer type. The electronic properties of the X substituents depend on the type of spacer, its length and the Y group, which is an example of the reverse SE. To describe how the SE between Y and X weakens with n, two approaches were compared: cSAR and SESE (SE stabilization energy). The EDDB (electron density of delocalized bonds) characterize changes in electron delocalization in spacers due to the SE. A new approach - EDDB differential maps - allow to extract the effect of X substitution on the electron delocalization. The charges at spacer's C atoms correlate with cSAR; changes in the slopes confirm the charge transfer by resonance.

Metal Influence on Cyaphide-Azide 1,3-Dipolar Cycloaddition Reactions: Aromaticity and Activation Strain

The factors governing 1,3-dipolar cycloaddition reactions involving C≡P-containing compounds are computationally explored in detail using quantum chemical tools. To this end, the parent process involving tBuN3 and tBuCP is analyzed and compared to the analogous reaction involving organometallic cyaphide complexes (metal=Au, Pt, Ge, Mg), in order to understand the role of the metal fragment in such transformations. It is found that while the metal fragment does not significantly influence the aromaticity of the corresponding concerted transition states or the regioselectivity of the transformation, it may modify the reactivity of the cyaphide complexes (i. e. Ge and Mg cyaphide complexes are comparatively more reactive). The computed reactivity trends and the factors behind the regioselectivity of the cycloaddition reaction are quantitatively analyzed with the help of the activation strain model in combination with the energy decomposition analysis method.

Theoretical investigation of tautomerism of 2- and 4-pyridones: origin, substituent and solvent effects

Computational investigation at the BHandHLYP/6-311+G(d,p) level of theory of the gas-phase tautomerism of 2- and 4-pyridones confirmed the slight prevalence of lactim in the case of the former, but its dominance in the case of the latter, as shown previously. Examination of aromaticity by using HOMA, EDDB, NBOdel, NICS and AICD led to the conclusion that tautomerization of 4-pyridone results in greater aromaticity gain. It is also driven by the Pauli repulsion relief, which was revealed by the tautomerization energy decomposition analysis. By contrast, in the case of 2-pyridone, lactim is favoured by orbital and electrostatic interactions and disfavoured by the Pauli repulsion. Aromaticity gain in this case is smaller. The position of the tautomeric equilibrium can be modulated by substituent inductive effects (Cl and F), inductive and resonance effects (NH2 and NO2), hydrogen bonding (NO2), and medium polarity, the increase of which increases lactam population.

Redox Activity of IrIII Complexes with Multidentate Ligands Based on Dipyrido-Annulated N-Heterocyclic Carbenes: Access to High Valent and High Spin State with Carbon Donors

Synthetic strategies to access high-valent iridium complexes usually require use of π donating ligands bearing electronegative atoms (e. g. amide or oxide) or σ donating electropositive atoms (e. g. boryl or hydride). Besides the η5 -(methyl)cyclopentadienyl derivatives, high-valent η1 carbon-ligated iridium complexes are challenging to synthesize. To meet this challenge, this work reports the oxidation behavior of an all-carbon-ligated anionic bis(CCC-pincer) IrIII complex. Being both σ and π donating, the diaryl dipyrido-annulated N-heterocyclic carbene (dpa-NHC) IrIII complex allowed a stepwise 4e- oxidation sequence. The first 2e- oxidation led to an oxidative coupling of two adjacent aryl groups, resulting in formation of a cationic chiral IrIII complex bearing a CCCC-tetradentate ligand. A further 2e- oxidation allowed isolation of a high-valent tricationic complex with a triplet ground state. These results close a synthetic gap for carbon-ligated iridium complexes and demonstrate the electronic tuning potential of organic π ligands for unusual electronic properties.

Substituent Effect versus Aromaticity─A Curious Case of Fulvene Derivatives

A computational study on amino- and nitro-substituted penta- and heptafulvenes reveals the interplay between the aromaticity and the substituent effect (SE). Ring substitution alone has little influence on the aromaticity, but in combination with an exo substituent of opposite properties, it substantially enhances the cyclic π-electron delocalization. Despite the SE being stronger for β substitution, only γ substitution leads to higher aromaticity. An explanation is provided by the electron density of delocalized bonds (EDDB) method, which proves to be a valuable tool in analyzing both cyclic delocalization and the SE.

Single─Not Double─3D-Aromaticity in an Oxidized Closo Icosahedral Dodecaiodo-Dodecaborate Cluster

3D-aromatic molecules with (distorted) tetrahedral, octahedral, or spherical structures are much less common than typical 2D-aromatic species or even 2D-aromatic-in-3D systems. Closo boranes, [BnHn]2- (5 ≤ n ≤ 14) and carboranes are examples of compounds that are singly 3D-aromatic, and we now explore if there are species that are doubly 3D-aromatic. The most widely known example of a species with double 2D-aromaticity is the hexaiodobenzene dication, [C6I6]2+. This species shows π-aromaticity in the benzene ring and σ-aromaticity in the outer ring formed by the iodine substituents. Inspired by the hexaiodobenzene dication example, in this work, we explore the potential for double 3D-aromaticity in [B12I12]0/2+. Our results based on magnetic and electronic descriptors of aromaticity together with 11B{1H} NMR experimental spectra of boron-iodinated o-carboranes suggest that these two oxidized forms of a closo icosahedral dodecaiodo-dodecaborate cluster, [B12I12] and [B12I12]2+, behave as doubly 3D-aromatic compounds. However, an evaluation of the energetic contribution of the potential double 3D-aromaticity through homodesmotic reactions shows that delocalization in the I12 shell, in contrast to the 10σ-electron I62+ ring in the hexaiodobenzene dication, does not contribute to any stabilization of the system. Therefore, the [B12I12]0/2+ species cannot be considered as doubly 3D-aromatic.

Excited-State (Anti)Aromaticity Explains Why Azulene Disobeys Kasha's Rule

Fluorescence exclusively occurs from the lowest excited state of a given multiplicity according to Kasha's rule. However, this rule is not obeyed by a handful of anti-Kasha fluorophores whose underlying mechanism is still understood merely on a phenomenological basis. This lack of understanding prevents the rational design and property-tuning of anti-Kasha fluorophores. Here, we propose a model explaining the photophysical properties of an archetypal anti-Kasha fluorophore, azulene, based on its ground- and excited-state (anti)aromaticity. We derived our model from a detailed analysis of the electronic structure of the ground singlet, first excited triplet, and quintet states and of the first and second excited singlet states using the perturbational molecular orbital theory and quantum-chemical aromaticity indices. Our model reveals that the anti-Kasha properties of azulene and its derivatives result from (i) the contrasting (anti)aromaticity of its first and second singlet excited states (S1 and S2, respectively) and (ii) an easily accessible antiaromaticity relief pathway of the S1 state. This explanation of the fundamental cause of anti-Kasha behavior may pave the way for new classes of anti-Kasha fluorophores and materials with long-lived, high-energy excited states.

π-Aromaticity Dominating in a Saturated Ring: Neutral Aromatic Silicon Analogues of Cyclobutane-1,3-diyls

The synthesis, structures, and reactivity of the first neutral 2π-aromatic Si4 rings [LSiSiAr(X)]2 (3: X = Br; 4: X = Cl; L = PhC(NtBu)2, Ar = 2,4,6-Me3C6H2) were described. Compounds 3 and 4 were obtained by 1,3-halogenation of tetrasilacyclobutadiene (LSiSiAr)2 (2), which was prepared by the reductive cross-coupling of trisilane (ArSiCl2)2SiHAr with two equiv of chlorosilylene LSiCl. The reaction of 3 with two equiv of PhLi yielded the corresponding substitution Si4 ring [LSiSiAr(Ph)]2 (5). Single-crystal X-ray diffraction analysis of 3 disclosed that it adopts both puckered (3a) and planar (3b) structures in the solid state, whereas 4 and 5 exhibit only a puckered structure. DFT calculations suggested that the puckered 3a features almost the same electronic structure with fully delocalized 2π planar 3b. The dominant 2π-aromaticity of 3 in a σ-frame has been demonstrated by DFT calculations, providing the first example of aromatics featuring both planar and puckered structures.

Isolable Tetragold(0) Clusters with Polarity-Tunable exo-Au-Au Bond via Intramolecular σ-Aromatization

Intramolecular π-aromatization is a trait of many organic compounds that enhances the stability of their structures and polarizes related C-C π bonds. In contrast, rare study is focused on this phenomenon in metal clusters. Many existing homometallic clusters exhibit aromaticity, often characterized by nonpolar metal-metal bonds and a high degree of symmetry. However, synthesizing low-symmetric homometallic clusters with high-polar metal-metal bonds is challenging due to their limited thermodynamic stability. Herein, we report a facile strategy for the synthesis of [Au(μ2 -ER2 )]3 -AuPMe3 (E=Ge, Sn; R2 =1,1,4,4-tetrakis(trimethylsilyl)butane-1,4-diyl) clusters and reveal a novel stabilization mode, intramolecular σ-aromatization. Our electronic structure analyses show that these low-symmetric clusters possess a ten-electron σ-aromatic system, which is achieved via intramolecular σ-aromatization. Moreover, the strength of σ-aromaticity gives rise to a polarity-tunable exo-Au-Au bond.

Actinide-Actinide Bonding: Electron Delocalisation and σ-Aromaticity in the Tri-Thorium Cluster [{Th(η8 -C8 H8 )(μ-Cl)2 }3 K2 ]

The tri-thorium cluster [{Th(η8 -C8 H8 )(μ3 -Cl)2 }3 {K(THF)2 }2 ]∞ (Nature 2021, 598, 72-75) was reported to feature intriguing σ-aromatic bonding between the thorium atoms, a mode of metal-metal bonding unique in the actinide series. However, the presence of this bonding motif has since been challenged by others. Here, we computationally explore electron delocalisation in a molecular cluster fragment of [{Th(η8 -C8 H8 )(μ3 -Cl)2 }3 {K(THF)2 }2 ]∞ and examine its responses to an applied magnetic field using a variety of methods. We also discuss the importance of the choice of basis set for the Th atoms and issues regarding locating QTAIM bond critical points. When taken together, the computed data consistently suggest the presence of delocalised Th-Th bonding and Th3 σ-aromaticity.

Overcrowded Triply Fused Carbo[7]helicene

This paper presents the synthesis and comprehensive analysis of a highly contorted and doubly negatively curved multihelicene compound, composed of three carbo[7]helicene units fused within a central six-membered ring. The synthesis of this compound involved a [2 + 2 + 2] cycloaddition reaction of 13,14-picyne, employing a Ni(0) catalyst, which exhibited superior performance compared to conventional Pd(0) catalysts. The evaluation of aromaticity in this triple carbo[7]helicene, utilizing magnetic and electronic criteria, led to noteworthy insights challenging the limitations of Clar's model of aromaticity.

Benzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole) and Its Bromo Derivatives: Molecular Structure and Reactivity

Benzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole) (isoBBT) is a new electron-withdrawing building block that can be used to obtain potentially interesting compounds for the synthesis of OLEDs and organic solar cells components. The electronic structure and delocalization in benzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole), 4-bromobenzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole), and 4,8-dibromobenzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole) were studied using X-ray diffraction analysis and ab initio calculations by EDDB and GIMIC methods and were compared to the corresponding properties of benzo[1,2-c:4,5-c']bis[1,2,5]thiadiazole (BBT). Calculations at a high level of theory showed that the electron affinity, which determines electron deficiency, of isoBBT was significantly smaller than that of BBT (1.09 vs. 1.90 eV). Incorporation of bromine atoms improves the electrical deficiency of bromobenzo-bis-thiadiazoles nearly without affecting aromaticity, which increases the reactivity of these compounds in aromatic nucleophilic substitution reactions and, on the other hand, does not reduce the ability to undergo cross-coupling reactions. 4-Bromobenzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole) is an attractive object for the synthesis of monosubstituted isoBBT compounds. The goal to find conditions for the selective substitution of hydrogen or bromine atoms at position 4 in order to obtain compounds containing a (het)aryl group in this position and to use the remaining unsubstituted hydrogen or bromine atoms to obtain unsymmetrically substituted isoBBT derivatives, potentially interesting compounds for organic photovoltaic components, was not set before. Nucleophilic aromatic and cross-coupling reactions, along with palladium-catalyzed C-H direct arylation reactions for 4-bromobenzo[1,2-d:4,5-d']bis([1,2,3]thiadiazole), were studied and selective conditions for the synthesis of monoarylated derivatives were found. The observed features of the structure and reactivity of isoBBT derivatives may be useful for building organic semiconductor-based devices.

Rupturing aromaticity by periphery overcrowding

The balance between strain relief and aromatic stabilization dictates the form and function of non-planar π-aromatics. Overcrowded systems are known to undergo geometric deformations, but the energetically favourable π-electron delocalization of their aromatic ring(s) is typically preserved. In this study we incremented the strain energy of an aromatic system beyond its aromatic stabilization energy, causing it to rearrange and its aromaticity to be ruptured. We noted that increasing the steric bulk around the periphery of π-extended tropylium rings leads them to deviate from planarity to form contorted conformations in which aromatic stabilization and strain are close in energy. Under increasing strain, the aromatic π-electron delocalization of the system is broken, leading to the formation of a non-aromatic, bicyclic analogue referred to as 'Dewar tropylium'. The aromatic and non-aromatic isomers have been found to exist in rapid equilibrium with one another. This investigation demarcates the extent of steric deformation tolerated by an aromatic carbocycle and thus provides direct experimental insights into the fundamental nature of aromaticity.

Helical Trilayer Nanographenes with Tunable Interlayer Overlaps

The synthesis of well-defined nanocarbon multilayers, beyond the bilayer structure, is still a challenging goal. Herein, two trilayer nanographenes were synthesized by covalently linking nanographene layers through helicene bridges. The structural characterization of the trilayer nanographenes revealed a compact trilayer-stacked architecture. The introduction of a furan ring into the helicene linker modulates the interlayer overlap and π-conjugation of the trilayer nanographenes, enabling the tuning of the interlayer coupling, as demonstrated by optical, electrochemical, and theoretical analyses. Both synthesized trilayer nanographenes are rigid chiral nanocarbons and show a chirality transfer from the helicene moiety to the stacked nanographene layers. These helical trilayer nanographenes reported here represent the covalently linked multilayer nanographenes rather than bilayer ones, showing the tunable multilayer stacking structure.

Switch from local to global aromatic character in Möbius carbon nanobelts upon dioxidation. Evaluation of magnetic behavior in neutral and charged species

Here we show that recent Möbius carbon nanobelts (MCNBs) can be switched from a local to a global aromatic behavior upon dioxidation. Hence, large aromatic structures can be achieved by the choice of the charge states, giving rise to shielding cone characteristics extended within the overall structure at the nanoscale regime, pushing the limit of aromatic circuits to 198 π-electrons.

Aromaticity controls the excited-state properties of host-guest complexes of nanohoops

π-Conjugated organic molecules have exciting applications as materials for batteries, solar cells, light-emitting diodes, etc. Among these systems, antiaromatic compounds are of particular interest because of their smaller HOMO-LUMO energy gap compared to aromatic compounds. A small HOMO-LUMO gap is expected to facilitate charge transfer in the systems. Here we report the ground and excited-state properties of two model nanohoops that are nitrogen-doped analogs of recently synthesized [4]cyclodibenzopentalenes - tetramers of benzene-fused aromatic 1,4-dihydropyrrolo[3,2-b]pyrrole ([4]DHPP) and antiaromatic pyrrolo[3,2-b]pyrrole ([4]PP). Their complexes with C60 fullerene show different behavior upon photoexcitation, depending on the degree of aromaticity. [4]DHPP acts as an electron donor, whereas [4]PP is a stronger electron acceptor than C60. The ultrafast charge separation combined with the slow charge recombination that we found for [4]PP⊃C60 indicates a long lifetime of the charge transfer state.

Isolation of an Annulated 1,4-Distibabenzene Diradicaloid

The first 1,4-distibabenzene-1,4-diide compound [(ADC)Sb]₂ (5) based on an anionic dicarbene (ADC) (ADC=PhC{N(Dipp)C}₂ , Dipp=2,6-iPr₂ C₆ H₃ ) is reported as a bordeaux-red solid. Compound 5, featuring a central six-membered C₄ Sb₂ ring with formally Sb^I atoms may be regarded as a base-stabilized cyclic bis-stibinidene in which each of the Sb atoms bears two lone-pairs of electrons. 5 undergoes 2 e-oxidation with Ph₃ C[B(C₆ F₅ )₄ ] to afford [(ADC)Sb]₂ [B(C₆ F₅ )₄ ]₂ (6) as a brick-red solid. Each of the Sb atoms of 6 has an unpaired electron and a lone-pair. The broken-symmetry open-shell singlet diradical solution for (6)²⁺ is calculated to be 2.13 kcal mol^-1 more stable than the closed-shell singlet. The diradical character of (6)²⁺ according to SS-CASSCF (state-specific complete active space self-consistent field) and UHF (unrestricted Hartree-Fock) methods amounts to 36 % and 39 %, respectively. Treatments of 6 with (PhE)₂ yield [(ADC)Sb(EPh)]₂ [B(C₆ F₅ )₄ ]₂ (7-E) (E=S or Se). Reaction of 5 with (cod)Mo(CO)₄ affords [(ADC)Sb]₂ Mo(CO)₄ (8).

Direct Conversion of N2 and O2 to Nitric Oxide at Room Temperature Initiated by Double Aromaticity in the Y2BO+ Cation

The conversion of dinitrogen to more useful and reactive molecules has been the focus of intense research by chemists. In contrast to reductive N₂ fixation, direct oxidation of N₂ by O₂ to nitric oxide under mild conditions via a thermochemical process is extremely challenging. Herein, we report the first example of N₂ and O₂ activation and coupling under thermochemical conditions through the remarkable ability of Y₂BO⁺ to react with one N₂ and two O₂ molecules. Detailed mechanistic studies using mass spectrometry and quantum chemical calculations revealed that the N₂ activation by Y₂BO⁺ is facilitated by the double aromatic character of the Y₂BON₂⁺ intermediate. Subsequent oxidation with O₂ releases NO in a dearomatization process driven by the formation of stronger Y-O bonds over the Y-N bonds. Our findings represent the first example of N₂ and O₂ activation and coupling under thermochemical conditions at room temperature, providing a novel strategy for small-molecule activation.

Superatomic icosahedral-CnB12-n (n = 0, 1, 2) Stuffed mononuclear and binuclear borafullerene and borospherene nanoclusters with spherical aromaticity

Boron and boron-based nanoclusters exhibit unique structural and bonding patterns in chemistry. Extensive density functional theory calculations performed in this work predict the mononuclear walnut-like Ci C50B54 (1) (C2B10@C48B44), C1 C50B54 (2) (CB11@C49B43), and S10 C50B54 (3) (B12@C50B42) which contain one icosahedral-CnB12-n core (n = 0, 1, 2) at the center following the Wade’s skeletal electron counting rules and the approximately electron sufficient binuclear peanut-like Cs C88B78 (4) ((C2B10)2@C84B58), Cs C88B78 (5) ((CB11)2@C86B56), Cs C88B78 (6) ((B12)2@C88B54), Cs B180 (7) ((B12)2@B156), Cs B182 (8) ((B12)2@B158), and Cs B184 (9) ((B12)2@B160) which encapsulate two interconnected CnB12-n icosahedrons inside. These novel core–shell borafullerene and borospherene nanoclusters appear to be the most stable species in thermodynamics in the corresponding cluster size ranges reported to date. Detailed bonding analyses indicate that the icosahedral B122−, CB11−, and C2B10 cores in these core–shell structures possess the superatomic electronic configuration of 1S21P61D101F8, rendering spherical aromaticity and extra stability to the systems. Such superatomic icosahedral-CnB12-n stuffed borafullerenes and borospherenes with spherical aromaticity may serve as embryos to form bulk boron allotropes and their carbon-boron binary counterparts in bottom-up approaches.

Aromaticity Concepts Derived from Experiments

Aromaticity, a very important term in organic chemistry, has never been defined unambiguously. Various ways to describe it come from different phenomena that have been experimentally observed. The most important examples related to some theoretical concepts are presented here.

Predicting dinitrogen activation by borenium and borinium cations

The activation of thermodynamically stable and kinetically inert dinitrogen (N₂) has been a great challenge due to the significant strength of the triple bond. Recently, in an experimental study on N₂ activation by boron species, a highly reactive two-coordinated borylene broke through the limitations of traditional strategies of N₂ activation by metal species. Still, studies on metal-free N₂ activation remain underdeveloped. Here, we systematically investigate a frustrated Lewis pair (FLP) combining carbene and borenium (or borinium) cations to screen potential candidates for N₂ activation via density functional theory calculations. As a result, we found that two FLPs (closed form FLP, borenium and open form FLP, borinium) are able to activate N₂ in a thermodynamically and kinetically favorable manner, with a low energy barrier of 9.6 and 7.3 kcal mol^-1, respectively. Furthermore, aromaticity was found to play an important role in the stabilization of the products, supported by nucleus-independent chemical shift (NICS), anisotropy of the current-induced density (ACID) and electron density of delocalized bonds (EDDB) analysis. Our findings provide an alternative approach for metal-free N₂ activation, highlighting the importance of FLP chemistry and aromaticity in N₂ activation.

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.

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.

Energetic and Geometric Characteristics of Substituents, Part 3: The Case of NO2 and NH2 Groups in Their Mono-Substituted Derivatives of Six-Membered Heterocycles

Substituted heterocyclic arenes play important roles in biochemistry, catalysis, and in the design of functional materials. Exemplary six-membered heteroaromatic molecules, that differ from benzene by inclusion of one heteroatom, are pyridine, phosphorine, arsabenzene, and borabenzene. This theoretical study concerns the influence of the heteroatom present in these molecules on the properties of substituents of two types: electron-donating (ED) NH2 group and electron-accepting (EA) NO2 group, attached at the 2-, 3-, or 4-position. The effect is evaluated by the energy of interaction (Erel) between the substituent and the substituted system and electronic properties of the substituents described by the charge of the substituent active region (cSAR) index. In addition, several geometric descriptors of the substituent and heteroaromatic ring, as well as changes in the aromaticity, are considered. The latter are assessed using the Electron Density of Delocalized Bonds (EDDBs) property of delocalized π electrons. The obtained results show that the electronegativity (EN) of the heteroatom has a profound effect on the EA/ED properties of the substituents. This effect is also reflected in the geometry of studied molecules. The Erel parameter indicates that the relative stability of the molecules is highly related to the electronic interactions between the substituent and the heteroarene. This especially applies to the enhancement or weakening of π-resonance due to the EN of the heteroatom. Additionally, in the 2-heteroarene derivatives, specific through-space ortho interactions contribute to the heteroatom effects.

Revival of Hückel Aromatic (Poly)benzenoid Subunits in Triplet State Polycyclic Aromatic Hydrocarbons by Silicon Substitution

By employing density functional theory (DFT) calculations we show that mono- and disilicon substitution in polycyclic aromatic hydrocarbons, having two to four benzene units, quenches their triplet state antiaromaticity by creating Hückel aromatic (poly)benzenoid subunit(s) and weakly antiaromatic, or almost nonaromatic silacycle. Therefore, such systems are predicted to be globally aromatic in both the ground state and the first excited triplet state. Putting the silicon atom(s) into various positions of a hydrocarbon provides an opportunity to tune the singlet-triplet energy gaps. They depend on the global aromaticity degree which, in turn, depends on the type of aromatic carbocyclic subunit(s) and the extent of their aromaticity. On the basis of the set of studied compounds, some preliminary rules on how to regulate the extent of global, semiglobal and local aromaticity are proposed. The results of this work extend the importance of Hückel aromaticity concept to excited triplet states which are usually characterized by the Baird type of (anti)aromaticity.

Metalla-phenalene complexes: synthesis, structure and aromaticity

Metallaaromatics show a diversity of aromaticity. In this work, we report the synthesis and structural characterization of the first rhena-phenalene complexes. In addition to the Hückel aromaticity and σ-aromaticity, pseudo π anti-aromaticity is observed. DFT computations show that this anti-aromaticity (paramagnetic properties) is induced by the fused aromatic naphthyl ring.

A Conjugated Figure-of-Eight Oligoparaphenylene Nanohoop with Adaptive Cavities Derived from Cyclooctatetrathiophene Core

A fully conjugated figure-of-eight nanohoop is presented with facile synthesis. The molecule's lemniscular skeleton features the combination of two strained oligoparaphenylene loops and a flexible cyclooctatetrathiophene core. Its rigid yet guest-adaptive cavities enable the formation of the peanut-like 1:2 host-guest complexes with C₆₀ or C₇₀ , which have been confirmed by X-ray crystallography and characterized in solution. Further computational studies suggest notable geometric variations and non-covalent interactions of the cavities upon binding with different fullerenes, as well as overall conjugation comparable to cycloparaphenylenes.

Valence non-Lewis density as an approach to describe and measure aromaticity of organic and inorganic molecules

Based on the linear combination of atomic orbital-molecular orbital by the natural bond orbitals (NBO) theory, the attractive donor-acceptor superposition interaction between filled (Lewis-type) and vacant (non-Lewis-type) orbitals provide a general mechanism for quantal energy lowering. This interaction has a direct impact on the quantity of the second-order stabilization energy. Therefore, the valence non-Lewis density (VNLD) index, the electron density of unoccupied valence nonbonding and antibonding orbitals, is introduced as an approach to describe and measure aromaticity. This index is based on the frontier orbital concept. To investigate the validity of the proposed aromaticity index, we selected several test sets of organic and inorganic molecules such as different ring sizes in cyclic and heterocyclic hydrocarbons, and all-metal and semimetal clusters, and compared our findings with previous aromaticity analysis. According to the results, VNLD values are well correlated and anticipated the order of aromaticity with the formerly introduced criteria. Furthermore, VNLD reveals that the rings with π-sextet electrons localized in a ring are more aromatic than the other rings, thus, it is in line with Clar's rule. Our proposed aromaticity index has advantages such as, easy to obtain from NBO analysis, and does not require reference molecules which made it more applicable for realizing the aromaticity order in many organic and inorganic compounds.