Direct estimate of conjugation and aromaticity in cyclic compounds with the EDA method

Aromatic stabilization energies in excited states at the multiconfigurational level: assessment in archetypal organic rings

In this study, the excited state (anti)aromaticity of archetypal rings: benzene, cyclobutadiene, and cyclooctatetraene, was investigated using the energetic criterion by calculating aromatic stabilization energies. Calculations were performed at the multiconfigurational level, including dynamic correlation effect corrections using the N-electron valence state perturbation theory (NEVPT2) method. Results were compared with previously published data based on the magnetic and delocalization criteria. Aromaticity was assessed for the ground state, singlet excited states (S1, S2, and S3), and triplet excited states (T1, T2, T3, and T4). (Anti)aromaticity assignments using the energetic criterion demonstrate both agreement and discrepancies with the other criteria, particularly for higher energy electronic states demonstrating the complexity of aromaticity assignment beyond the ground state. Finally, an approximate equation is proposed for the calculation of aromatic stabilization energies in excited states using experimental data such as formation enthalpies and well-resolved absorption spectra.

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.

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.

Annulated Azuleno[2,1,8-ija]azulenes: Synthesis and Properties

Non-alternant non-benzenoid hydrocarbons exhibit very different optical and electronic properties than their well-studied benzenoid analogues. However, preparing such structures with extended conjugation length, remains challenging. Herein, we report the synthesis and properties of azuleno[2,1,8-ija]azulene derivatives using a two-step sequence involving a four-fold aldol condensation between aromatic dialdehydes and readily available tetrahydropentalene-2,5-(1H,3H)-dione. Molecules with band gap values ranging from 1.69 to 2.14 eV and molar extinction coefficients (ϵ) of nearly 3×10⁵  M^-1  cm^-1 have been prepared. These annulene-like structures exhibit significant diatropic ring currents (aromatic), as supported by ¹ H NMR spectroscopy and DFT calculations. Field-effect transistors (OFETs) using azuleno[2,1,8-ija]azulene derivatives as semiconductors exhibit charge mobility values of up to 0.05 cm²  V^-1  s^-1 .

Benzene and Borazine, so Different, yet so Similar: Insight from Experimental Charge Density Analysis

Although benzene and borazine are isoelectronic and isostructural, they have very different electronic structures, mainly due to the polar nature of the B-N bond. Herein, we present an experimental study of the charge density distribution obtained from the multipole model formalism and Hirshfeld atom refinement (HAR) based on high-resolution X-ray diffraction data of borazine B₃N₃H₆ (1) and B,B',B″-trichloroborazine (2) crystals. These data are compared to those obtained from HAR for benzene (4) and 1,3,5-trichlorobenzene (5) and further compared with values obtained from density functional theory calculations in the gas phase, where N,N',N″-trichloroborazine (3) was also included. The results confirm that, unlike benzene, borazines are only weakly aromatic with an island-like electronic delocalization within the B₃N₃ ring involving only the nitrogen atoms. Furthermore, delocalization indices and interacting quantum atom energy for bonded and non-bonded atoms were found to be highly suitable indicators capable of describing the origin of the discrepancies observed when the degree of aromaticity in 2 and 3 is evaluated using common aromaticity indices. Additionally, analysis of intermolecular interactions in the crystals brings further evidence of a weakly aromatic character of the borazines as it reveals surprising similarities between the crystal packing of borazine and benzene and also between B,B',B″-trichloroborazine and 1,3,5-trichlorobenzene.

The hidden aromaticity in borazine

The aromaticity of borazine (B3N3H6), also known as the inorganic benzene, is a controversial issue since this compound has several characteristics that could qualify it as an aromatic compound. However, recent studies using magnetic criteria indicate that this compound should be considered as a non-aromatic system. This assignment is mainly due to diatropic currents in the nitrogen atoms without observation of ring currents. The present work shows by means of the magnetic criteria that borazine has a ring current hidden by the local contributions of degenerate orbitals π1 and π2. Additionally, the study of borazine's first triplet state antiaromaticity using the magnetic and energetic criteria by means of isomerization stabilization energies (ISEs) together with Baird's and Hückel's rules suggests that borazine is best described as an (weakly) aromatic system.

Aromaticity in catalysis: metal ligand cooperation via ligand dearomatization and rearomatization

Unlike the conventional model of transition metal catalysis, ligands in metal-ligand cooperative (or bifunctional) catalysis are involved in the substrate activations. Such processes have offered unique mechanistic understandings and led to new concepts for the catalyst design. In particular, unprecedented activities were discovered when the ligand could undergo dearomatization-rearomatization reactions during the catalytic cycle. Aromatization can provide an extra driving force to thermodynamics; consequently, it brings a new perspective to ligand platform design for catalysis. While numerous applications were demonstrated, the influences of changing ligand aromatic properties were often overlooked. In this article, representative ligand systems will be highlighted and a comparison between the Milstein and the Huang pincer systems will be discussed to provide theoretical and conceptual insights.

Searching for double σ- and π-aromaticity in borazine derivatives

Inspired by the double-aromatic (σ and π) C6H3 +, C6I6 2+, and C6(SePh)6 2+ ring-shaped compounds, herein we theoretically study their borazine derivative analogues. The systems studied are the cation and dications with formulas B3N3H3 +, B3N3Br6 2+, B3N3I6 2+, B3N3(SeH)6 2+, and B3N3(TeH)6 2+. Our DFT calculations indicate that the ring-shaped planar structures of B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ are more stable in the singlet state, while those of B3N3Br6 2+ and B3N3(SeH)6 2+ prefer the triplet state. Besides, exploration of the potential energy surface shows that the ring-shaped structure is the putative global minimum only for B3N3I6 2+. According to chemical bonding analysis, B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ have σ and π delocalized bonds. The number of delocalized σ/π electrons is 2/6 for the first, and 10/6 for the second and third, similar to what their carbon analogs exhibit. Finally, the analysis of the magnetically induced current density allows B3N3H3 +, B3N3I6 2+, and B3N3(TeH)6 2+ to be classified as strongly σ aromatic, and poorly π aromatic compounds.

Stacking interactions of borazine: important stacking at large horizontal displacements and dihydrogen bonding governed by electrostatic potentials of borazine

Potential energy surfaces of borazine-benzene and borazine-borazine stacking interactions were studied by performing DFT, CCSD(T)/CBS and SAPT calculations. The strongest borazine-benzene stacking was found in a parallel-displaced geometry, with a CCSD(T)/CBS interaction energy of -3.46 kcal mol-1. The strongest borazine-borazine stacking has a sandwich geometry, with a CCSD(T)/CBS interaction energy of -3.57 kcal mol-1. The study showed that borazine forms significant stacking interactions at large horizontal displacements (over 4.5 Å), with energies of -2.20 kcal mol-1 for the borazine-benzene and -1.96 kcal mol-1 for the borazine-borazine system. The strength of interactions and their geometrical preferences can be rationalized by observing the electrostatic potentials of borazine and benzene, which is in agreement with SAPT analysis showing that electrostatics is the most important energy component for borazine stacking. All the interactions found in crystal structures of borazine and related compounds were identified either as potential curve minima or the geometries obtained from their optimizations. We also report a new dihydrogen bonding dimer with a CCSD(T)/CBS interaction energy of -2.37 kcal mol-1, which is encountered in the borazine crystal structures and enables the formation of additional simultaneous interactions that contribute to the overall stability of the crystals.

Recent Advances in Metallaaromatic Chemistry

Metallaaromatics can be broadly defined as aromatic compounds in which one of the ring atoms is a transition metal. The metallabenzenes are one important class of these compounds that has undergone extensive study recently. Closely related species such as fused-ring metallabenzenes, heterometallabenzenes, π-coordinated metallabenzenes and metallabenzynes have also attracted considerable attention. Although many metallaaromatics can be considered as metalla-analogues of classic organic aromatic compounds, this is not always the case. Recent seminal studies have shown that metallapentalenes and metallapentalynes, which are metalla-analogues of the anti-aromatic compounds pentalene and pentalyne, are in fact aromatic and highly stable. Very unusual spiro-metallaaromatic compounds have also recently been isolated. In this concepts article, key features of all these intriguing metallaaromatic compounds are discussed with reference to the structural, spectroscopic, reactivity and theoretical studies that have been undertaken. These compounds continue to generate much interest, not only because of the contributions they make to fundamental chemical understanding, but also because of the promise of possible practical applications.

Theoretical study of azido gauche effect and its origin

The strength and origin of the azido gauche effect were studied by ab initio calculations and compared with the well-known fluorine gauche effect.

Testing the effectiveness of the isoelectronic substitution principle through the transformation of aromatic osmathiophene derivatives into their inorganic analogues

The isoelectronic substitution principle has been employed in the [OsCl₂(SC₃H₃)(PH₃)₂]⁺ complex, generating a new set of complexes, where C has been replaced by B⁻ and/or N⁺.

Planar vs. three-dimensional X62−, X2Y42−, and X3Y32− (X, Y = B, Al, Ga) metal clusters: an analysis of their relative energies through the turn-upside-down approach

Despite the fact that B and Al belong to the same group 13 elements, the B₆²⁻ cluster prefers the planar D_2h geometry, whereas Al₆²⁻ favours the O_h structure, which is caused by orbital interactions.

Dyotropic rearrangement of bridgehead substituents in closed dithienylethenes; conjugated verses non-conjugated analogues

Type I dyotropic rearrangement reactions of halogen and methyl substituents at the bridgehead position of diarylethenes and dihydroarylethenes have been studied through density functional theory at B3LYP/6-31+G(d) level. The calculations have been performed to explore the dyotropic rearrangement as a possible factor for the elusive nature of halogenated dithienylethenes (closed). The dyotropic rearrangement process in closed dithienylethenes is then compared with the dihydro analogues. Moreover, the effect of hetero atom and conjugation is also explored through quantum mechanical calculations.

Triplet State Aromaticity: NICS Criterion, Hyperconjugation, and Charge Effects

Aromaticity, one of the most important concepts in organic chemistry, has attracted considerable interest from both experimentalists and theoreticians. It remains unclear which NICS index is best to evaluate the triplet-state aromaticity. Here, we carry out thorough density functional theory (DFT) calculations to examine this issue. Our results indicate that among the various computationally available NICS indices, NICS(1)zz is the best for the triplet state. The correlations can be improved from 0.840 to 0.938 when only neutral species are considered, demonstrating the significant effect of the charge on the triplet-state aromaticity. In addition, calculations suggest that five-membered cyclic species with "hyperconjugative" aromaticity (and antiaromaticity) in the S0 state will become antiaromatic (and aromatic) in the T1 state, indicating an important role of hyperconjugation. Finally, a moderate correlation (r(2) =0.708) is identified between the NICS(1)zz values and spin distributions.

Study of Electron Delocalization in 1,2-, 1,3-, and 1,4-Azaborines Based on the Canonical Molecular Orbital Contributions to the Induced Magnetic Field and Polyelectron Population Analysis

The electron delocalization in 1,2-azaborine, 1,3-azaborine, and 1,4-azaborine is studied using canonical molecular orbital contributions to the induced magnetic field (CMO-IMF) method and polyelectron population analysis (PEPA). Contour maps of the out-of-plane component of the induced magnetic field (Bz(ind)) of the π system show that the three azaborines, in contrast with borazine, sustain much of benzene's π-aromatic character. Among them, 1,3-azaborine exhibits the strongest π delocalization, while 1,4-azaborine is the weakest. Contour maps of Bz(ind) for individual π orbitals reveal that the differentiation of the magnetic response among the three isomers originates from the π-HOMO orbitals, whose magnetic response is governed by rotational allowed transitions to unoccupied orbitals. The low symmetry of azaborines enables a paratropic response from HOMO to unoccupied orbitals excitations, with their magnitude depending on the shape of interacting orbitals. 1,3-Azaborine presents negligible paratropic contributions to Bz(ind) from HOMO to unoccupied orbitals transitions, where 1,2- and 1,4-azaborine present substantial paratropic contributions, which lead to reduced diatropic response. Natural bond orbital (NBO) analysis employing PEPA shows that only the 1,3-azaborine contains π-electron fully delocalized resonance structures.

Combined activation strain model and energy decomposition analysis methods: a new way to understand pericyclic reactions

The recently introduced activation strain model (ASM) has allowed us to gain more insight into the intimacies of different fundamental processes in chemistry. In combination with the energy decomposition analysis (EDA) method, we have nowadays a very useful tool to quantitatively understand the physical factors that govern the activation barriers of reactions within organic and organometallic chemistry. In this Perspective article, we present selected illustrative examples of the application of this method to pericyclic reactions (Diels-Alder and double group transfer reactions) to show that this methodology nicely complements other more traditional, widely used theoretical methods.

Is it possible to achieve a complete desaturation of cycloalkanes promoted by o-benzyne?

The competence between the possible reaction pathways in the reaction of different cycloalkanes and o-benzyne has been computationally explored. Systems able to undergo complete desaturation by consecutive double hydrogen atom transfer reactions are suggested.

Introducing “carborazine” as a novel heterocyclic aromatic species

Aromatic carborazine is a heterocyclic ring consisting of C, B and N atoms at opposite faces.

Aromaticity of metallabenzenes and related compounds

In this review, we focus on the aromaticity of a particular family of organometallic compounds known as metallabenzenes, which are characterized by the formal replacement of a CH group in benzene by an isolobal transition metal fragment.

Energy decomposition analysis of gauche preference in 2-haloethanol, 2-haloethylamine (halogen = F, Cl), their protonated forms and anti preference in 1-chloro-2-fluoroethane

Small, electronegative elements contribute more electrostatic and orbital stabilization to the anti → gauche isomerization, and greater steric repulsion. The first and the latter actually oppose our traditional view of conformational equilibria.

Energy decomposition analysis approaches and their evaluation on prototypical protein–drug interaction patterns

The partitioning of the interaction energy into chemical components such as electrostatics, polarization, and charge transfer is possible with energy decomposition analysis approaches. We review and evaluate these for biomolecular applications.

Aromaticity of azines through dyotropic double hydrogen transfer reaction

Density functional theory calculations have been performed at B3LYP/6-31+G (d) level to quantify the aromaticities of mono- to triazines through dyotropic double hydrogen transfer (DDHT) reaction. The reaction was chosen such that the azines are products of double hydrogen dyotropic rearrangement, and activation barriers and energies of the reactions were functions of the aromaticities of azines. Small activation barriers and high energies of reactions were characteristic of the reactions delivering highly aromatic azines. Synchronicity, reaction energies and energies of activation have been analyzed, and the aromaticity values obtained thereof were compared with the aromaticity values from other geometric and magnetic criteria. Energies of activation were found superior to the energies of reaction for the determination of the aromaticities. Aromaticities of most of the azines were comparable to the aromaticity of benzene. Activation barriers and reaction energies for the dyotropic reactions delivering contiguous or polynitrogeneous azines had thermodynamic contributions arising from the contiguous nature of azines, in addition to the aromaticity related thermodynamic contributions. Moreover, the aromaticity values of azines are also affected by the fusion of azine to the reaction center. When corrected for these factors, the aromaticities of azines using energies of activation for DDHT correlated nicely with the aromaticities of azines reported in the literature through NICS(0)πzz and some other energetic methods.

Origin of Reactivity Trends of Noble Gas Endohedral Fullerenes Ng2@C60 (Ng = He to Xe)

We have computationally studied the factors governing the enhanced Diels-Alder reactivity of noble gas endohedral fullerenes Ng2@C60 when Ng = Ar, Kr, and Xe as compared to Ng = none, He, and Ne. To this end, we have employed the activation strain model of reactivity in combination with the energy decomposition analysis (EDA) method in DFT calculations on the Diels-Alder cycloaddition reaction between Ng2@C60 and 1,3-butadiene. Our results indicate that when heavier noble gas dimers are introduced inside the C60 cage, dramatic effects on both the geometrical and electronic structure of the fullerenic cage occur leading to a remarkable enhanced interaction between the deformed reactants along the entire reaction coordinate.

Ene-ene-yne reactions: activation strain analysis and the role of aromaticity

The trend in reactivity of the thermal cycloisomerization reactions of 1,3-hexadien-5-ynes, A=B-C=D-E≡F, were explored and analyzed by using density functional theory at the M06-2X/def2-TZVPP level. These reactions proceed through formally aromatic transition states to form a bent-allene intermediate with relatively high activation barriers. Activation-strain analyses show that the major factor controlling this Hopf cyclization is the geometrical strain energy associated with the rotation of the terminal [A] group. This rotation is necessary for achieving a favorable HOMO-LUMO overlap with the yne-moiety [F] associated with the formation of the new A-F single bond. In addition, the relationship between the aromaticity of the corresponding cyclic transition states (all six-membered rings) and the computed activation barriers were analyzed. The calculations also indicate that the aromatization of the bent-allene structures takes place through two consecutive 1,2-hydrogen shifts, the second one exhibiting negligible energy barriers.