All-Metal Aromaticity: Revisiting the Ring Current Model among Transition Metal Clusters

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.

Conflicting Aromaticity in Trirhodium(I) Rosarin

Aromaticity is one of the most important and widely used concepts in chemistry. Among the various experimentally discovered and theoretically predicted compounds that possess different types of aromaticity, conflicting aromaticity, where aromatic and antiaromatic electron delocalization is present in one molecule simultaneously, remains one of the most controversial and elusive concepts, although theoretically predicted 15 years ago. In this work, we synthesized a novel conflicting aromatic trirhodium complex that contains a σ-aromatic metal fragment surrounded by the π-antiaromatic organic ligand and characterized it by nuclear magnetic resonance spectroscopy, high-resolution mass spectrometry, and X-ray single crystal structure analysis. Experimental characterization and quantum chemical calculations confirm the unique conflicting aromaticity of the synthesized trirhodium molecule. Thus, this novel conflicting aromatic molecule expands the family of aromatic compounds. This discovery will enable researchers to develop and understand the phenomena of conflicting aromaticity in chemistry.

Multiconfigurational actinide nitrides assisted by double Möbius aromaticity

Understanding the bonding nature between actinides and main-group elements remains a key challenge in actinide chemistry due to the involvement of f orbitals. Herein, we propose a unique "aromaticity-assisted multiconfiguration" (AAM) model to elucidate the bonding nature in actinide nitrides (An2N2, An = Ac, Th, Pa, U). Each planar four-membered An2N2 with equivalent An-N bonds possesses four delocalized π electrons and four delocalized σ electrons, forming a new family of double Möbius aromaticity that contributes to the molecular stability. The unprecedented aromaticity further supports actinide nitrides to exhibit multiconfigurational characters, where the unpaired electrons (2, 4 or 6 in naked Th2N2, Pa2N2 or U2N2, respectively) either are spin-free and localized on metal centres or form metal-ligand bonds. High-level multiconfigurational computations confirm an open-shell singlet ground state for actinide nitrides, with small energy gaps to high spin states. This is consistent with the antiferromagnetic nature observed experimentally in uranium nitrides. The novel AAM bonding model can be authenticated in both experimentally identified compounds containing a U2N2 motif and other theoretically modelled An2N2 clusters and is thus expected to be a general chemical bonding pattern between actinides and main-group elements.

The σ+π dual aromaticity of typical bi-tetrazole ring molecule TKX-50

Two complexes of dihydroxylammonium 5,5'-bistetrazole-1,1'-diolate (TKX-50) were employed to evaluate the aromaticity of their tetrazole rings via deep analysis such as the electronic structure, the ZZ component of the natural chemical shielding tensor (NICSZZ) and component orbitals, localized orbital locator purely contributed by σ-orbitals (LOL-σ) and localized orbital locator purely contributed by π-orbitals (LOL-π), the anisotropy of the induced current density (AICD) and the ZZ component of iso-chemical shielding surface (ICSSZZ) of these tetrazole rings thereof. The conclusion shows: that all tetrazole rings and bi-tetrazole rings in complexes have strong σ and a comparable strength π double aromaticity; all these magnetic shields almost symmetrically increase from the central axis to the tetrazole ring atoms; tetrazole rings in complex II show a little stronger dual aromaticity than that in complex I mainly due to the different orientation of the fragment 2 encompassing two hydroxylamine groups resulting in different effects on the contributions of σ orbitals and π orbitals to total aromaticity of tetrazole rings thereof; the difference in aromaticity is fundamentally caused by the atoms O with stronger electron-withdrawing than atom N in fragment 2 interact with bi-tetrazole ring through O in complex I but through N in complex II.

Advances for Triangular and Sandwich-Shaped All-Metal Aromatics

Much experimental work has been contributed to all-metal σ, π and δ-aromaticity among transition metals, semimetallics and other metals in the past two decades. Before our focused investigations on the properties of triangular and sandwich-shaped all-metal aromatics, A. I. Boldyrev presented general discussions on the concepts of all-metal σ-aromaticity and σ-antiaromaticity for metallo-clusters. Schleyer illustrated that Nucleus-Independent Chemical Shifts (NICS) were among the most authoritative criteria for aromaticity. Ugalde discussed the earlier developments of all-metal aromatic compounds with all possible shapes. Besides the theoretical predictions, many stable all-metal aromatic trinuclear clusters have been isolated as the metallic analogues of either the σ-aromatic molecule's [H3]+ ion or the π-aromatic molecule's [C3H3]+ ion. Different from Hoffman's opinion on all-metal aromaticity, triangular all-metal aromatics were found to hold great potential in applications in coordination chemistry, catalysis, and material science. Triangular all-metal aromatics, which were theoretically proved to conform to the Hückel (4n + 2) rule and possess the smallest aromatic ring, could also play roles as stable ligands during the formation of all-metal sandwiches. The triangular and sandwich-shaped all-metal aromatics have not yet been specifically summarized despite their diversity of existence, puissant developments and various interesting applications. These findings are different from the public opinion that all-metal aromatics would be limited to further applications due to their overstated difficulties in synthesis and uncertain stabilities. Our review will specifically focus on the summarization of theoretical predictions, feasible syntheses and isolations, and multiple applications of triangular and sandwich shaped all-metal aromatics. The appropriateness and necessities of this review will emphasize and disseminate their importance and applications forcefully and in a timely manner.

Group 13 five- and six-membered rings with rare 2π aromaticity

According to Hückel's rule, cyclic species are aromatic if they have 4n + 2 (n = 0, 1, 2, etc.) π electrons. However, large aromatic rings (atom number > 4) with minimum 2π electrons (i.e., n = 0) are rather rare because of the structural instability stemming from the deficient π electrons compared to the ring size. To date, the largest 2π-aromatic ring is a five-membered Ga5 core reported in a recent experiment. Herein, density functional theory calculations predicted seven inverse-sandwich Na2(MH)5 and half-sandwich Ca(MH)5 or Ca(MH)6 (M = Al, Ga or In) structures. They all have planar central five- or six-membered Al/Ga/In rings, rather negative binding energies and large HOMO-LUMO gaps. Their dianionic metal rings exhibit obvious aromatic characters and appreciable diatropic ring currents due to the good delocalization of 2π electrons donated by the Na/Ca metals. Interestingly, they also have novel collective bonding with the Na/Ca atoms interacting with both the metal ring and surrounding H atoms. These metalloaromatic rings not only greatly enrich the precious family of 2π aromatics, but also increase the maximum ring atom number from five to six, thus paving the way for Hückel-type 2π-aromatic large rings.

Magnetic Antiaromaticity─Paratropicity─Does Not Necessarily Imply Instability

Magnetically induced ring currents are a conventional tool for the characterization of aromaticity. Dia- and paratropic currents are thought to be associated with stabilization (aromaticity) and destabilization (antiaromaticity), respectively. In the present work, I have questioned the validity of the paratropic currents as a measure of antiaromaticity among monocyclic hydrocarbons. I have shown that while reduced/oxidized radical ions of hydrocarbons sustain strong paratropic currents, they often gain extra stabilization via cyclic conjugation compared to their acyclic counterparts.

Planar pentacoordinate s-block metals

The presence of a delocalized π-bond is often considered an essential criterion for achieving planar hypercoordination. Herein, we show that σ-delocalization could be sufficient to make the planar configuration the most stable isomer in a series of planar pentacoordinate s-block metals. High-level ab initio computations reveal that the global minimum of a series of interalkali and interalkali-alkaline earth clusters (LiNa5, Li5Mg+, Na5Mg+, K5Ca+, CaRb5+, Rb5Sr+, and SrCs5+) adopts a singlet D5h structure with a planar pentacoordinate lithium or alkaline earth metal (AE = Mg, Ca, Sr). These clusters are unusual combinations to stabilize a planar pentacoordinate atom, as all their constituents are electropositive. Despite the absence of π-electrons, Hückel's rule is fulfilled by the six σ-electrons. Furthermore, the systems exhibit a diatropic ring current in response to an external magnetic field and a strong magnetic shielding, so they might be classified as σ-aromatic. Therefore, multicenter σ-bonds and the resulting σ-delocalization stabilize these clusters, even though they lack π-aromaticity.

Engineering magic number Au19 and Au20 cage structures using electron withdrawing atoms

Gold cages are a subset of gold nanoparticles and these structures are of major interest due to their favourable physiochemical properties. In order for these structures to be useful in applications, they must be chemically stable. The objective of this research is to transform non-magic number cage structures into magic number cage structures by the addition of electron-withdrawing groups on the cages. The electronic properties for Au₁₉X and Au₂₀X₂ (X = F, Cl, Br, I) are calculated and observed. It is expected that the electron-withdrawing groups will remove the electron density from the gold cages and leave them positively charged. We first optimize the geometries of the initial gold cages and verify the structures are a local minima. From there, we attach our halogens to the gold cages and optimize the structures to determine the NICS values and HOMO-LUMO gaps. NICS values were found to be more negative when a more electronegative halogen was used. Calculations have found that Au₁₉F and Au₂₀F₂ have the most negative NICS values, indicating greater spherical aromaticity. Iodine, being the least electronegative atom, had the most positive NICS value and smallest HOMO-LUMO gap. All calculations are compared to the magic cluster Au₁₈ which satisfies Hirsh's 2(N + 1)² rule for n = 2.

Planar Four-Membered Diboron Actinide Compound with Double Möbius Aromaticity

The Möbius rule predicts that a planar four-membered metallacycle can be aromatic with four mobile electrons, but such a simple ring has escaped recognition because it usually favors Hückel anti-aromaticity. Here, we report that a quasi-square four-membered actinide compound (Pa₂B₂) is doubly Möbius aromatic. Chemical bonding analyses reveal that this diboron protactinium molecule has four delocalized π electrons in addition to four delocalized σ electrons, satisfying the 4n Möbius rule for both σ and π components. Energetically, the block-localized wavefunction method, which is the simplest variant of ab initio valence bond theory, shows that the delocalization energy for the π and σ electrons reaches up to 65.0 and 72.3 kcal/mol, respectively, while the extra cyclic resonance energy (ECRE) amounts to 45 kcal/mol. The large positive ECRE values strongly confirm the unprecedented double Möbius aromaticity in Pa₂B₂. We anticipate that this new type of aromatic molecule can enrich the concept of Möbius aromaticity and open a new avenue for actinide compounds.

Planar σ-Aromaticity in Ga-Doped Au Clusters

Recently, the first example of Au-Ga clusters is synthesized and characterized, which can be described by the jellium model as a superatom with 8 valence electrons that come from the joint contribution of Au and Ga atoms, opening a whole new field for further research. Here, the structure features and stability of one Ga-doped Au cluster with magic number electrons (6 and 8) are analyzed in detail. Moreover, the valence electron fillings and chemical bonding of them are also further explored. It is found that Au₃Ga and Au₅Ga clusters present planar configurations, and they have higher stability than that of neighbor clusters. The AIMD simulations show that these two clusters still have a good thermal stability at 500 K. The molecular orbital analyses show that the Au₃Ga and Au₅Ga have three and one typical delocalization orbital throughout the whole planar spaces, respectively, following the planar σ-aromaticity rule. The ELF and LOL analyses are further performed, and the results are consistent with the molecular orbital analyses. The NICSzz-scan curves confirm that the Au₃Ga is more aromatic than the Au₅Ga, and the reason is that the former has more delocalized electrons than the latter. Our work opens up aromaticity studies in the Au-Ga clusters.

φ-Aromaticity in prismatic {Bi6}-based clusters

The occurrence of aromaticity in organic molecules is widely accepted, but its occurrence in purely metallic systems is less widespread. Molecules comprising only metal atoms (M) are known to be able to exhibit aromatic behaviour, sustaining ring currents inside an external magnetic field along M-M connection axes (σ-aromaticity) or above and below the plane (π-aromaticity) for cyclic or cage-type compounds. However, all-metal compounds provide an extension of the electrons' mobility also in other directions. Here, we show that regular {Bi₆} prisms exhibit a non-localizable molecular orbital of f-type symmetry and generate a strong ring current that leads to a behaviour referred to as φ-aromaticity. The experimentally observed heterometallic cluster [{CpRu}₃Bi₆]^-, based on a regular prismatic {Bi₆} unit, displays aromatic behaviour; according to quantum chemical calculations, the corresponding hypothetical Bi₆^2- prism shows a similar behaviour. By contrast, [{(cod)Ir}₃Bi₆] features a distorted Bi₆ moiety that inhibits φ-aromaticity.

Core-electron contributions to the molecular magnetic response

Orbital contributions to the magnetic response depend on the method used to compute them. Here, we show that dissecting nuclear magnetic shielding tensors using natural localized molecular orbitals (NLMOs) leads to anomalous core contributions. The arbitrariness of the assignment might significantly affect the interpretation of the magnetic response of nonplanar molecules such as C₆₀ or [14]helicene and the assessment of their aromatic character. We solve this problem by computing the core- and σ-components of the induced magnetic field (and NICS) and the magnetically induced current density by removing the valence electrons (RVE). We estimate the core contributions to the magnetic response by performing calculations on the corresponding highly charged molecules, such as C₆H₆³⁰⁺ for benzene, using gauge-including atomic orbitals and canonical molecular orbitals (CMOs). The orbital contributions to nuclear magnetic shielding tensors are usually estimated by employing a natural chemical shielding (NCS) analysis in NLMO or CMO bases. The RVE approach shows that the core contribution to the magnetic response is small and localized at the nuclei, contrary to what NCS calculations suggest. This may lead to a completely incorrect interpretation of the magnetic σ-orbital response of nonplanar structures, which may play a major role in the overall magnetic shielding of the system. The RVE approach is thus a simple and inexpensive way to determine the magnetic response of the core- and σ-electrons.

The smallest 4f-metalla-aromatic molecule of cyclo-PrB2− with Pr–B multiple bonds

The concept of metalla-aromaticity proposed by Thorn–Hoffmann (Nouv. J. Chim. 1979, 3, 39) has been expanded to organometallic molecules of transition metals that have more than one independent electron-delocalized system. Lanthanides, with highly contracted 4f atomic orbitals, are rarely found in multiply aromatic systems. Here we report the discovery of a doubly aromatic triatomic lanthanide-boron molecule PrB₂⁻ based on a joint photoelectron spectroscopy and quantum chemical investigation. Global minimum structural searches reveal that PrB₂⁻ has a C_2v triangular structure with a paramagnetic triplet ³B₂ electronic ground state, which can be viewed as featuring a trivalent Pr(III,f²) and B₂⁴⁻. Chemical bonding analyses show that this cyclo-PrB₂⁻ species is the smallest 4f-metalla-aromatic system exhibiting σ and π double aromaticity and multiple Pr–B bonding characters. It also sheds light on the formation of the rare B₂⁴⁻ tetraanion by the high-lying 5d orbitals of the 4f-elements, completing the isoelectronic B₂⁴⁻, C₂²⁻, N₂, and O₂²⁺ series.

We report the smallest 4f-metalla-aromatic molecule of PrB₂⁻ exhibiting σ and π double aromaticity and multiple Pr–B bond characters.

Aromaticity of Heterocirculenes

This review summarizes the results on the aromaticity of a series of synthesized and hypothetical neutral heterocirculene molecules and their double charged ions. The aromaticity of heterocirculenes is a direct reflection of their electronic structure responsible for the specific optoelectronic and photophysical properties. We show how the presence of a heteroatom in the outer macrocycle affects the aromaticity of hetero[8]circulenes. In addition, we also describe the change in aromaticity and strain energy for a series of the “lower” (n < 8) and “higher” (n > 8) hetero[n]circulenes. It was demonstrated that the loss of planarity with increased strain leads to an increased antiaromaticity of the lower hetero[n]circulenes, whereas higher hetero[n]circulenes demonstrate a more pronounced aromatic nature because of the small departure from planarity of each heteroarene ring in hetero[n]circulene molecule. Finally, we discuss the aromatic nature of the first examples of π-extended hetero[8]circulenes.

Inter/Intramolecular Cascade of 1,6-Enynes Catalyzed by All-Metal Aromatic Tripalladium Complexes and Carboxylic Acids

Trinuclear all-metal aromatic clusters are an original class of molecules with a cyclic and planar metal core. Characterized by peculiar metal-metal delocalized bonds, they represent a new frontier in transition-metal catalysis. We report a study on C-C-forming reactions of polyunsaturated substrates catalyzed by trinuclear all-metal aromatic palladium clusters. The synthesis of two new families of tricyclic compounds was obtained with a broad functional group tolerance under mild reaction conditions. A peculiar regio- and diastereoselectivity characterized the method, demonstrating that trinuclear palladium complexes are complementary to their popular mononuclear peers. Furthermore, preliminary studies on the mechanism of these polycyclization reactions revealed unique features of the homogeneous catalytic system.

Exploring water adsorption and reactivity in a series of doped aluminum cluster anions

We present a systematic density functional study of central- and surface-doped aluminum cluster anions Al₁₂X^- (X = Mg, B, Ga, Si, P, Sc-Zn), their interactions and reactivity with water. Adsorption of water molecules on central-doped clusters is governed by the cluster electron affinity. Doping introduces a dramatic change in the cluster electronic structure by virtue of different ordering and occupation of super-atomic shells, which leads to the creation of complementary active sites controlling the reactivity with water. Surface doping creates unequal charge distribution on the cluster surface, resulting in the adsorption and reactivity of surface-doped clusters being dominated by electrostatic effects. These results demonstrate the strong influence of the doping position on the nature of the interaction and reactivity of the cluster, and contribute to a better understanding of doping effects.

The B2 Structural Motif as a Tool for Modulating Ring Currents in Monocyclic Li Clusters

Magnetically induced current densities, calculated at the M06-2X/def2-TZVP level using the diamagnetic-zero version of the continuous transformation of origin of current density (CTOCD-DZ) method, were employed to study the aromaticity in Li3B2− and Li4B2. It was found that the Li3/Li4 rings in Li3B2− and Li4B2 remarkably resemble the monocyclic Li3+ and Li42+ clusters. Unlike the parent Li3+ and Li42+ systems that sustain negligibly weak global current density circulation, the Li3B2− and Li4B2 clusters exhibit a strong diatropic current density. The present work demonstrates how structural modifications introduced by the B2 unit can be used for modulating the current density in cyclic Li-based clusters.

Unraveling the stability of cyclobutadiene complexes using aromaticity markers

Cyclobutadiene (CBD) is the paradigmatic antiaromatic molecule but is known to form highly stable aromatic complexes, e.g. CBD-Fe(CO)3. This intriguing reversal of aromaticity from antiaromatic to aromatic terrain during the complexation process cannot be appropriately handled with single-reference-based theoretical techniques. We explore this aromaticity reversal, for the first time, by a detailed aromaticity analysis using magnetically induced current densities (MICD) and nucleus independent chemical shifts (NICS) using genuine ab initio multi-reference wavefunction-based theory. We trace the dramatic change of aromaticity for a prototypical cyclobutadiene complex, CBD-CH+ (CH+Fe(CO)3), considering a 3D potential energy surface for two independent parameters, namely the approach of CH+ and the automerization cross-section of cyclobutadiene. The 3D potential energy surfaces indicate the presence of a conical intersection/avoided crossing between the ground and the first excited state. The plot of aromaticity indices and the corresponding numerical values show that the change of aromaticity indices is drastic around the conical intersection/avoided crossing and automerization of cyclobutadiene plays a crucial role in the formation of cyclobutadiene complexes. Computations on analogous CBD-Be and CBD-CO systems (Be/COFe(CO)3) emphasize the generality of the conclusions drawn from the CBD-CH+ system.

Cubic aromaticity in ligand-stabilized doped Au superatoms

The magnetic response of valence electrons in doped gold-based M@Au₈L₈ ^q superatoms (M = Pd, Pt, Ag, Au, Cd, Hg, Ir, and Rh; L = PPh₃; and q = 0, +1, +2) is studied by calculating the gauge including magnetically induced currents (GIMIC) in the framework of the auxiliary density functional theory. The studied systems include 24 different combinations of the dopant, total cluster charge, and cluster structure (cubic-like or oblate). The magnetically induced currents (both diatropic and paratropic) are shown to be sensitive to the atomic structure of clusters, the number of superatomic electrons, and the chemical nature of the dopant metal. Among the cubic-like structures, the strongest aromaticity is observed in Pd- and Pt-doped M@Au₈L₈ ⁰ clusters. Interestingly, Pd- and Pt-doping increases the aromaticity as compared to a similar all-gold eight-electron system Au₉L₈ ⁺¹. With the recent implementation of the GIMIC in the deMon2k code, we investigated the aromaticity in the cubic and butterfly-like M@Au₈ core structures, doped with a single M atom from periods 5 and 6 of groups IX-XII. Surprisingly, the doping with Pd and Pt in the cubic structure increases the aromaticity compared to the pure Au case not only near the central atom but encompassing the whole metallic core, following the aromatic trend Pd > Pt > Au. These doped (Pd, Pt)@Au₈ nanoclusters show a closed shell 1S²1P⁶ superatom electronic structure corresponding to the cubic aromaticity rule 6n + 2.

Relating nucleus independent chemical shifts with integrated current density strengths

Indices based on the nucleus independent chemical shift (NICS) are the most frequently used in analysis of magnetic aromaticity. The magnetically induced current density, on the other hand, is a key concept in defining magnetic aromaticity. The integrated current strength (current strength susceptibility) was found to be a very useful tool in aromaticity studies. There is widely accepted notion that the properly chosen NICS-based index can provide information on the current density strength and direction in a molecule of interest. In this work, a detailed numerical testing of the relationship between the integrated bond current strength and the most employed NICS indices was performed for a set of 43 monocyclic aromatic molecules. Based on the statistical data analysis, the relationship between the bond current strength and its π and σ electron components, on one side, and the isotropic NICS (NICSiso and NICSπ,iso) and zz-component of the NICS tensor (NICSzz and NICSπ,zz), on the other side, was examined. It was found that between the NICSπ,zz(1) and π-electron bond current strenghts there is very good linear correlation. Quite surprisingly, it was revealed that the NICSiso(1) and NICSzz(1) are not correlated with the π electron bond current strengths. On the other hand, a reasonably good linear correlation was found between the NICSzz(1) and total bond current strengths.

On the NICS limitations to predict local and global current pathways in polycyclic systems

Here, we analyze the possibility of predicting local and global current densities in a series of bicyclic hydrocarbons with 4n and 4n + 2 π-electrons from the nucleus-independent chemical shifts (NICS) computations.

Theoretical Insight into the Structural Nonplanarity in Aromatic Fused-Ring Metallabenzenes

Metalla-aromatics have attracted considerable attention due to their fascinating structural and reactive properties as well as their potential as prospective functional materials. Metallabenzenes and their fused-ring counterparts are significant members of metalla-aromatics, while their crystal structures often display seemly counterintuitive nonplanar geometry. The geometric bending of metallabenzenes has been attributed to the unfavorable antibonding interactions in the σ-space orbitals rather than the general opinion regarding the π-space orbitals of an aromatic compound. However, the origin of the geometric bending in fused-ring metallabenzenes remains elusive. In this work, we elucidated that such a "σ-control mechanism" still holds for fused-ring metallabenzenes by performing systematic calculations for a plethora of metallabenzenes and fused-ring metallabenzenes. Furthermore, we found that a more bent geometry can be achieved for fused-ring metallabenzenes than their corresponding metallabenzenes by fusing the aromatic rings at the ortho-position of a metal center to induce extra repulsion. The more significant bending in fused-ring metallabenzenes also favors the aromaticity enhancement. These findings not only provide mechanistic insight into the unexpected geometric distortion in both metallabenzenes and fused-ring metallabenzenes but also pave the way to design and develop bent metalla-aromatics with enhanced metalla-aromaticity, which hold great potential as aromatic functional materials.

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.

Magnetically-induced current density investigation in carbohelicenes and azahelicenes

A computational study of the magnetically-induced current (MIC) density has been carried out for a variety of ortho fused polycyclic aromatic molecules at the density functional theory level with the gauge including magnetically induced current (GIMIC) method. With this method, the aromatic character of each ring in a homologous series of carbohelicenes with an increasing number of fused benzene rings is assessed and compared with other aromaticity criteria such as the Nucleus Independent Chemical Shift [NICS(0), NICSzz(0)] and Bond Length Alternation (BLA) parameters. All criteria indicate that the two outer rings are the most aromatic ones [i.e. higher induced current, more negative NICS(0) and NICSzz(0) values, and smaller BLA values]. For the large helicenes (n > 10), the current drops along the following four rings and then rises again. Additionally, we have proven that this behavior is not due to a difference of the local magnetic field coming from a difference of orientation of the ring with respect to the external magnetic field (oriented along the helical axis). Upon fusing additional benzene rings to form hexa-peri-hexabenzo[7]helicene, some rings (B, D, and F) are a lot less aromatic (even non-aromatic) than the others. The NICS(0) and NICSzz(0) values exaggerate this behavior because they are all positive values, which is a signature of antiaromaticity. Then, when substituting one, three, or four benzene rings with pyrrole ones to form mono-aza-[7]-helicene, tri-aza-[7]-helicene, and tetra-aza-[7]-helicene, remarkable changes in the electronic structures of the helicenes are observed. Indeed, the induced currents are always smaller in the pyrrole rings than in the benzene ones. This has been further investigated using the streamline and the color map representations, which indicate that the diatropic current density passing through the plane cutting the C-N bonds in the pyrrole rings is stronger but is more localized than the current density passing through the plane cutting the C-C bonds in the benzene rings. This gives a positive but smaller total induced current for the pyrrole rings than for the benzene ones. For these systems, the NICS(0) and even more the NICSzz(0) values are not fully reliable to probe the local aromaticity, contrary to the induced current. Indeed, the NICSzz(0) values for the tri-aza-[7]-helicene molecule range from -14.23 to 1.14 ppm, which cannot lead to the same conclusion as the induced current values (10.03 to 12.87 nA T-1).

Multiple-Valence Aluminum and the Electronic and Geometric Structure of Al nO m Clusters

Electronic stability in aluminum clusters is typically associated with either closed electronic shells of delocalized electrons or a +3 oxidation state of aluminum. To investigate whether there are alternative routes toward electronic stability in aluminum oxide clusters, we used theoretical methods to examine the geometric and electronic structure of Al _nO _m (2 ≤ n ≤ 7; 1 ≤ m ≤ 10) clusters. Electronically stable clusters with large HOMO-LUMO (highest occupied molecular orbital and lowest unoccupied molecular orbital) gaps were identified and could be grouped into two categories. (1) Al_{2 n}O_{3 n} clusters with a +3 oxidation state on the aluminum and (2) planar clusters including Al₄O₄, Al₅O₃, Al₆O₅, and Al₆O₆. The structures of the planar clusters have external Al atoms bound to a single O atom. Their electronic stability is explained by the multiple-valence Al sites, with the internal Al atoms having an oxidation state of +3, whereas the external Al atoms have an oxidation state of +1.

Norcorrole as a Delocalized, Antiaromatic System

Nickel norcorrole provides an unusual example of a molecule that is strongly antiaromatic according to the magnetic criterion, but which exhibits, according to high-quality DFT calculations, a symmetric, delocalized structure with no difference in bond length between adjacent C_meso-C_α bonds. A fragment molecular orbital analysis suggests that these discordant observations are a manifestation of the high stability of the dipyrrin fragments, which retain their electronic and structural integrity even as part of the norcorrole ring system.

Interpreting Aromaticity and Antiaromaticity through Bifurcation Analysis of the Induced Magnetic Field

In all molecules, a current density is induced when the molecule is subjected to an external magnetic field. In turn, this current density creates a particular magnetic field. In this work, the bifurcation value of the induced magnetic field is analyzed in a representative set of aromatic, non-aromatic and antiaromatic monocycles, as well as a set of polycyclic hydrocarbons. The results show that the bifurcation value of the ring-shaped domain adequately classifies the studied molecules according to their aromatic character. For aromatic and nonaromatic molecules, it is possible to analyze two ring-shaped domains, one diatropic (inside the molecular ring) and one paratropic (outside the molecular ring). Meanwhile, for antiaromatic rings, only a diatropic ring-shaped domain (outside the molecular ring) is possible to analyze, since the paratropic domain (inside the molecular ring) is irreducible with the maximum value (attractor) at the center of the molecular ring. In some of the studied cases, i. e., in heteroatomic species, bifurcation values do not follow aromaticity trends and present some inconsistencies in comparison to ring currents strengths, showing that this approximation provides only a qualitative estimation about (anti)aromaticity.

Electron Delocalization in Planar Metallacycles: Hückel or Möbius Aromatic?

In this work the relationship between the formal number of π-electrons, d-orbital conjugation topology, π-electron delocalization and aromaticity in d-block metallacycles is investigated in the context of recent findings concerning the correlation of π-HOMO topology and the magnetic aromaticity indices in these species. It is demonstrated that for π-electron rich d-metallacycles the direct link between aromaticity, the number of π-electrons and the frontier π-orbital topology does not strictly hold and for such systems it is very difficult to unambiguously associate their aromaticity with the "4n+2" (Hückel) and "4n" (Möbius) rules. It is also shown that the recently proposed electron density of delocalized bonds (EDDB) method can successfully be used not only to quantify and visualize aromaticity in such difficult cases, but also - in contrast to magnetic aromaticity descriptors - to provide a great deal of information on the real role of d-orbitals in metallacycles without the ambiguity of bookkeeping of electrons in the π-subsystem of the molecular ring. Interestingly, some of the metallacycles studied cannot be classified exclusively as Hückel or Möbius because they have a hybrid Hückel-Möbius or even quasi-aromatic nature.

Magnetically induced current density in triple-layered beryllium–boron clusters

Magnetically induced current densities reveal the double aromatic character of the examined Be–B clusters.

Diaqua-β-octaferrocenyltetraphenylporphyrin: a multiredox-active and air-stable 16π non-aromatic species

Herein the synthesis and properties of the first β-octaferrocenyltetraphenylporphyrin, {TPPFc₈(H₂O)₂}, in its extraordinary stable and non-aromatic 16π form are reported, showing seven separate reversible redox events.

Impact of heteroatoms (S, Se, and Te) on the aromaticity of heterocirculenes

Computations reveal the structural and energetic aspects of aromaticity in heterocirculenes.

Aromaticity of acenes: the model of migrating π-circuits

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

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

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

Calculations of current densities and aromatic pathways in cyclic porphyrin and isoporphyrin arrays

Magnetically induced current density susceptibilities have been studied for a number of cyclic ethyne- and butadiyne-bridged porphyrin and isoporphyrin arrays. The current density susceptibilities have been calculated using the gauge-including magnetically induced current (GIMIC) method, which is interfaced to the TURBOMOLE quantum chemistry code. Aromatic properties and current pathways have been analyzed and discussed by numerical integration of the current density susceptibilities passing selected chemical bonds yielding current strength susceptibilities. Despite the interrupted π-framework, zinc(ii) isoporphyrin sustains a ring current of ca. 10 nA T^-1. Porphyrin and isoporphyrin dimers sustain a significant current strength at the linker, whereas the larger porphyrinoid arrays sustain mainly local ring currents. Isoporphyrin dimers with saturated meso carbons have strong net diatropic ring-current strengths of 20 nA T^-1 fulfilling Hückels aromaticity rule. Porphyrin trimers and tetramers exhibit almost no current strength at the linker. The porphyrin moieties maintain their strong net diatropic ring current.