Quantitative Assessment of Aromaticity and Antiaromaticity Utilizing Vibrational Spectroscopy

Porous Nanographenes, Graphene Nanoribbons, and Nanoporous Graphene Selectively Synthesized from the Same Molecular Precursor

We demonstrate a family of molecular precursors based on 7,10-dibromo-triphenylenes that can selectively produce different varieties of atomically precise porous graphene nanomaterials through the use of different synthetic environments. Upon Yamamoto polymerization of these molecules in solution, the free rotations of the triphenylene units around the C-C bonds result in the formation of cyclotrimers in high yields. In contrast, in on-surface polymerization of the same molecules on Au(111) these rotations are impeded, and the coupling proceeds toward the formation of long polymer chains. These chains can then be converted to porous graphene nanoribbons (pGNRs) by annealing. Correspondingly, the solution-synthesized cyclotrimers can also be deposited onto Au(111) and converted into porous nanographenes (pNGs) via thermal treatment. Thus, both processes start with the same molecular precursor and end with a porous graphene nanomaterial on Au(111), but the type of product, pNG or pGNR, depends on the specific coupling approach. We also produced extended nanoporous graphenes (NPGs) through the lateral fusion of highly aligned pGNRs on Au(111) that were grown at high coverage. The pNGs can also be synthesized directly in solution by Scholl oxidative cyclodehydrogenation of cyclotrimers. We demonstrate the generality of this approach by synthesizing two varieties of 7,10-dibromo-triphenylenes that selectively produced six nanoporous products with different dimensionalities. The basic 7,10-dibromo-triphenylene monomer is amenable to structural modifications, potentially providing access to many new porous graphene nanomaterials. We show that by constructing different porous structures from the same building blocks, it is possible to tune the energy band gap in a wide range.

Unraveling the effect of aromaticity for the dynamics of excited states of single benzene fluorophores

The photophysical properties of a series of recently synthesized single benzene fluorophores were investigated using ensemble density functional theory calculations. The energetic stability of the ground and excited state species were counterposed against the aromaticity index derived from local vibrational modes. It was found that the large Stokes shift of the fluorophores (up to ca. 5800 cm  - 1 ) originates from the effect of electron donating and electron withdrawing substituents rather than π -delocalization and related (anti-)aromaticity. On the basis of nonadiabatic molecular dynamics simulations, the absence of fluorescence from one of the regioisomers was explained by the occurrence of easily accessible S  1 /S  0 conical intersections below the vertical excitation energy level. It is demonstrated in the manuscript that the analysis of local mode force constants and the related aromaticity index represent a useful tool for the characterization of π -delocalization effects in π -conjugated compounds.

Structural stability and the low-lying singlet and triplet states of BN-n-acenes, n = 1-7

The chemical stability and the low-lying singlet and triplet excited states of BN-n-acenes (n = 1-7) were studied using single reference and multireference methodologies. From the calculations, descriptors such as the singlet-triplet splitting, the natural orbital (NO) occupations and aromaticity indexes are used to provide structural and energetic analysis. The boron and nitrogen atoms form an isoelectronic pair of two carbon atoms, which was used for the complete substitution of these units in the acene series. The structural analysis confirms the effects originated from the insertion of a uniform pattern of electronegativity difference within the molecular systems. The covalent bonds tend to be strongly polarized which does not happen in the case of a carbon-only framework. This effect leads to a charge transfer between neighbor atoms resulting in a more strengthened structure, keeping the aromaticity roughly constant along the chain. The singlet-triplet splitting also agrees with this stability trend, maintaining a consistent gap value for all molecules. The BN-n-acenes molecules possess a ground state with monoconfigurational character indicating their electronic stability. The low-lying singlet excited states have charge transfer character, which proceeds from nitrogen to boron.

On-Surface Synthesis of C144 Hexagonal Coronoid with Zigzag Edges

Coronoids as polycyclic aromatic macrocycles enclosing a cavity have attracted a lot of attention due to their distinctive molecular and electronic structures. They can be also regarded as nanoporous graphene molecules whose electronic properties are critically dependent on the size and topology of their outer and inner peripheries. However, because of their synthetic challenges, the extended hexagonal coronoids with zigzag outer edges have not been reported yet. Here, we report the on-surface synthesis of C144 hexagonal coronoid with outer zigzag edges on a designed precursor undergoing hierarchical Ullmann coupling and cyclodehydrogenation on the Au(111) surface. The molecular structure is unambiguously characterized by bond-resolved noncontact atomic force microscopy imaging. The electronic properties are further investigated by scanning tunneling spectroscopy measurements, in combination with the density functional theory calculations. Moreover, the values of the harmonic oscillator model of aromaticity are derived from calculations that suggest that the molecular structure is ideally represented by Clar's model. Our results provide approaches toward realizing a hexagonal coronoid with zigzag edges, potentially inspiring fabrication of hexagonal zigzag coronoids with multiple radical characters in the future.

Mechanistic Details of the Sharpless Epoxidation of Allylic Alcohols—A Combined URVA and Local Mode Study

In this work, we investigated the catalytic effects of a Sharpless dimeric titanium (IV)–tartrate–diester catalyst on the epoxidation of allylalcohol with methyl–hydroperoxide considering four different orientations of the reacting species coordinated at the titanium atom (reactions R1–R4) as well as a model for the non-catalyzed reaction (reaction R0). As major analysis tools, we applied the URVA (Unified Reaction Valley Approach) and LMA (Local Mode Analysis), both being based on vibrational spectroscopy and complemented by a QTAIM analysis of the electron density calculated at the DFT level of theory. The energetics of each reaction were recalculated at the DLPNO-CCSD(T) level of theory. The URVA curvature profiles identified the important chemical events of all five reactions as peroxide OO bond cleavage taking place before the TS (i.e., accounting for the energy barrier) and epoxide CO bond formation together with rehybridization of the carbon atoms of the targeted CC double bond after the TS. The energy decomposition into reaction phase contribution phases showed that the major effect of the catalyst is the weakening of the OO bond to be broken and replacement of OH bond breakage in the non-catalyzed reaction by an energetically more favorable TiO bond breakage. LMA performed at all stationary points rounded up the investigation (i) quantifying OO bond weakening of the oxidizing peroxide upon coordination at the metal atom, (ii) showing that a more synchronous formation of the new CO epoxide bonds correlates with smaller bond strength differences between these bonds, and (iii) elucidating the different roles of the three TiO bonds formed between catalyst and reactants and their interplay as orchestrated by the Sharpless catalyst. We hope that this article will inspire the computational community to use URVA complemented with LMA in the future as an efficient mechanistic tool for the optimization and fine-tuning of current Sharpless catalysts and for the design new of catalysts for epoxidation reactions.

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.

Structural study of 1- and 2-naphthol: new insights into the non-covalent H-H interaction in cis-1-naphthol

Previous microwave studies of naphthol monomers were supplemented by measuring spectra of all ¹³C mono-substituted isotopologues of the cis- and trans-conformers of 1-naphthol and 2-naphthol in their natural abundances. The resulting data were utilized to determine substitution structures and so-called semi-experimental effective structures. Results from electronic structure calculations show that the OH group of cis-1-naphthol points ≈6° out of plane, which is consistent with the inertial defect data of cis- and trans-1-naphthol. The non-planarity of cis-1-naphthol is a result of a close-contact H-atom-H-atom interaction. This type of H-H interaction has been the subject of much controversy in the past and we provide here an in-depth theoretical analysis of it. The naphthol system is particularly well-suited for such analysis as it provides internal standards with its four different isomers. The methods used include quantum theory of atoms in molecules, non-covalent interactions, independent gradient model, local vibrational mode, charge model 5, and natural bond orbital analyses. We demonstrate that the close-contact H-H interaction is neither a purely attractive nor repulsive interaction, but rather a mixture of the two.

Self-Assembled Lanthanide Antenna Glutathione Sensor for the Study of Immune Cells

The small molecule 8-methoxy-2-oxo-1,2,4,5-tetrahydrocyclopenta[de]quinoline-3-carboxylic acid (2b) behaves as a reactive non-fluorescent Michael acceptor, which after reaction with thiols becomes fluorescent, and an efficient Eu³⁺ antenna, after self-assembling with this cation in water. This behavior makes 2b a highly selective GSH biosensor, which has demonstrated high potential for studies in murine and human cells of the immune system (CD4⁺ T, CD8⁺ T, and B cells) using flow cytometry. GSH can be monitored by the fluorescence of the product of addition to 2b (445 nm) or by the luminescence of Eu³⁺ (592 nm). 2b was able to capture baseline differences in GSH intracellular levels among murine and human CD4⁺ T, CD8⁺ T, and B cells. We also successfully used 2b to monitor intracellular changes in GSH associated with the metabolic variations governing the induction of CD4⁺ naïve T cells into regulatory T cells (T_REG).

Chemisorption-Induced Formation of Biphenylene Dimer on Ag(111)

We report an example that demonstrates the clear interdependence between surface-supported reactions and molecular-adsorption configurations. Two biphenyl-based molecules with two and four bromine substituents, i.e., 2,2'-dibromobiphenyl (DBBP) and 2,2',6,6'-tetrabromo-1,1'-biphenyl (TBBP), show completely different reaction pathways on a Ag(111) surface, leading to the selective formation of dibenzo[e,l]pyrene and biphenylene dimer, respectively. By combining low-temperature scanning tunneling microscopy, synchrotron radiation photoemission spectroscopy, and density functional theory calculations, we unravel the underlying reaction mechanism. After debromination, a biradical biphenyl can be stabilized by surface Ag adatoms, while a four-radical biphenyl undergoes spontaneous intramolecular annulation due to its extreme instability on Ag(111). Such different chemisorption-induced precursor states between DBBP and TBBP consequently lead to different reaction pathways after further annealing. In addition, using bond-resolving scanning tunneling microscopy and scanning tunneling spectroscopy, we determine with atomic precision the bond-length alternation of the biphenylene dimer product, which contains 4-, 6-, and 8-membered rings. The 4-membered ring units turn out to be radialene structures.

Porphyrinoids, a unique platform for exploring excited-state aromaticity

Recently, Baird (anti)aromaticity has been referred to as a description of excited-state (anti)aromaticity. With the term of Baird's rule, recent studies have intensively verified that the Hückel aromatic [4n + 2]π (or antiaromatic [4n]π) molecules in the ground state are reversed to give Baird aromatic [4n]π (or Baird antiaromatic [4n + 2]π) molecules in the excited states. Since the Hückel (anti)aromaticity has great influence on the molecular properties and reaction mechanisms, the Baird (anti)aromaticity has been expected to act as a dominant factor in governing excited-state properties and processes, which has attracted intensive scientific investigations for the verification of the concept of reversed aromaticity in the excited states. In this scientific endeavor, porphyrinoids have recently played leading roles in the demonstration of the aromaticity reversal in the excited states and its conceptual development. The distinct structural and electronic nature of porphyhrinoids depending on their (anti)aromaticity allow the direct observation of excited-state aromaticity reversal, Baird's rule. The explicit experimental demonstration with porphyrinoids has contributed greatly to its conceptual development and application in novel functional organic materials. Based on the significant role of porphyrinoids in the field of excited-state aromaticity, this review provides an overview of the experimental verification of the reversal concept of excited-state aromaticity by porphyrinoids and the recent progress on its conceptual application in novel functional molecules.

Multiple dehydrogenation of fluorene cation and neutral fluorene using the statistical molecular fragmentation model

The statistical molecular fragmentation (SMF) model was used to analyze the 306 fragmentation channels (containing 611 different species) that result from the fluorene (C13H10+) cation losing up to three hydrogen atoms (neutral radicals and/or a proton). Breakdown curves from such analysis permit one to extract experimentally inaccessible information about the fragmentation of the fluorene cation, such as the locations of the lost hydrogen atoms (or proton), yields of the neutral fragments, electronic states of the residues, and quantification of very low probability channels that would be difficult to detect. Charge localization during the fragmentation pathways was studied to provide a qualitative understanding of the fragmentation process. Breakdown curves for both the fluorene cation and neutral fluorene were compared. The SMF results match the rise and fall of the one hydrogen loss yield experimentally measured by imaging photoelectron-photoion coincidence spectroscopy using a VUV synchrotron.

Kekulene: On-Surface Synthesis, Orbital Structure, and Aromatic Stabilization

We revisit the question of kekulene's aromaticity by focusing on the electronic structure of its frontier orbitals as determined by angle-resolved photoemission spectroscopy. To this end, we have developed a specially designed precursor, 1,4,7(2,7)-triphenanthrenacyclononaphane-2,5,8-triene, which allows us to prepare sufficient quantities of kekulene of high purity directly on a Cu(111) surface, as confirmed by scanning tunneling microscopy. Supported by density functional calculations, we determine the orbital structure of kekulene's highest occupied molecular orbital by photoemission tomography. In agreement with a recent aromaticity assessment of kekulene based solely on C-C bond lengths, we conclude that the π-conjugation of kekulene is better described by the Clar model rather than a superaromatic model. Thus, by exploiting the capabilities of photoemission tomography, we shed light on the question which consequences aromaticity holds for the frontier electronic structure of a π-conjugated molecule.

Kekulene: On-Surface Synthesis, Orbital Structure, and Aromatic Stabilization

We revisit the question of kekulene's aromaticity by focusing on the electronic structure of its frontier orbitals as determined by angle-resolved photoemission spectroscopy. To this end, we have developed a specially designed precursor, 1,4,7(2,7)-triphenanthrenacyclononaphane-2,5,8-triene, which allows us to prepare sufficient quantities of kekulene of high purity directly on a Cu(111) surface, as confirmed by scanning tunneling microscopy. Supported by density functional calculations, we determine the orbital structure of kekulene's highest occupied molecular orbital by photoemission tomography. In agreement with a recent aromaticity assessment of kekulene based solely on C-C bond lengths, we conclude that the π-conjugation of kekulene is better described by the Clar model rather than a superaromatic model. Thus, by exploiting the capabilities of photoemission tomography, we shed light on the question which consequences aromaticity holds for the frontier electronic structure of a π-conjugated molecule.

Interplay of Ring Puckering and Hydrogen Bonding in Deoxyribonucleosides

The Cremer-Pople ring puckering analysis and the Konkoli-Cremer local mode analysis supported by the topological analysis of the electron density were applied for the first comprehensive analysis of the interplay between deoxyribose ring puckering and intramolecular H-bonding in 2'-deoxycytidine, 2'-deoxyadenosine, 2'-deoxythymidine, and 2'-deoxyguanosine. We mapped for each deoxyribonucleoside the complete conformational energy surface and the corresponding pseudorotation path. We found only incomplete pseudorotation cycles, caused by ring inversion, which we coined as pseudolibration paths. On each pseudolibration path a global and a local minimum separated by a transition state were identified. The investigation of H-bond free deoxyribonucleoside analogs revealed that removal of the H-bond does not restore the full conformational flexibility of the sugar ring. Our work showed that ring puckering predominantly determines the conformational energy; the larger the puckering amplitude, the lower the conformational energy. In contrast no direct correlation between conformational energy and H-bond strength was found. The longest and weakest H-bonds are located in the local minimum region, whereas the shortest and strongest H-bonds are located outside the global and local minimum regions at the turning points of the pseudolibration paths, i.e., H-bonding determines the shape and length of the pseudolibration paths. In addition to the H-bond strength, we evaluated the covalent/electrostatic character of the H-bonds applying the Cremer-Kraka criterion of covalent bonding. H-bonding in the puric bases has a more covalent character whereas in the pyrimidic bases the H-bond character is more electrostatic. We investigated how the mutual orientation of the CH₂OH group and the base influences H-bond formation via two geometrical parameters describing the rotation of the substituents perpendicular to the sugar ring and their tilting relative to the ring center. According to our results, rotation is more important for H-bond formation. In addition we assessed the influence of the H-bond acceptor, the lone pair (N, respectively O), via the delocalization energy. We found larger delocalization energies corresponding to stronger H-bonds for the puric bases. The global minimum conformation of 2'-deoxyguanosine has the strongest H-bond of all conformers investigated in this work with a bond strength of 0.436 which is even stronger than the H-bond in the water dimer (0.360). The application of our new analysis to DNA deoxyribonucleotides and to unnatural base pairs, which have recently drawn a lot of attention, is in progress.

Local aromaticity mapping in the vicinity of planar and nonplanar molecules

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

A Continuum from Halogen Bonds to Covalent Bonds: Where Do λ3 Iodanes Fit?

The intrinsic bonding nature of λ 3 -iodanes was investigated to determine where its hypervalent bonds fit along the spectrum between halogen bonding and covalent bonding. Density functional theory with an augmented Dunning valence triple zeta basis set ( ω B97X-D/aug-cc-pVTZ) coupled with vibrational spectroscopy was utilized to study a diverse set of 34 hypervalent iodine compounds. This level of theory was rationalized by comparing computational and experimental data for a small set of closely-related and well-studied iodine molecules and by a comparison with CCSD(T)/aug-cc-pVTZ results for a subset of the investigated iodine compounds. Axial bonds in λ 3 -iodanes fit between the three-center four-electron bond, as observed for the trihalide species IF 2 − and the covalent FI molecule. The equatorial bonds in λ 3 -iodanes are of a covalent nature. We explored how the equatorial ligand and axial substituents affect the chemical properties of λ 3 -iodanes by analyzing natural bond orbital charges, local vibrational modes, the covalent/electrostatic character, and the three-center four-electron bonding character. In summary, our results show for the first time that there is a smooth transition between halogen bonding → 3c–4e bonding in trihalides → 3c–4e bonding in hypervalent iodine compounds → covalent bonding, opening a manifold of new avenues for the design of hypervalent iodine compounds with specific properties.

Demonstration of Baird's rule complementarity in the singlet state with implications for excited-state intramolecular proton transfer

Greater aromaticity in the ground state leads to greater antiaromaticity in the excited state (and vice versa) which helps rationalize previously unexplained behavior of ESIPT fluorophores.

Spectroscopic Diagnosis of Excited-State Aromaticity: Capturing Electronic Structures and Conformations upon Aromaticity Reversal

Aromaticity, the special energetic stability derived from cyclic [4 n + 2]π-conjugated electronic structures, has been the topic of intense interest in chemistry because it plays a critical role in rationalizing molecular stability, reactivity, and physical/chemical properties. Recently, the pioneering work by Colin Baird on aromaticity reversal, postulating that aromatic (antiaromatic) character in the ground state reverses to antiaromatic (aromatic) character in the lowest excited triplet state, has attracted much scientific attention. The completely reversed aromaticity in the excited state provides direct insight into understanding the photophysical/chemical properties of photoactive materials. In turn, the application of aromatic molecules to photoactive materials has led to numerous studies revealing this aromaticity reversal. However, most studies of excited-state aromaticity have been based on the theoretical point of view. The experimental evaluation of aromaticity in the excited state is still challenging and strenuous because the assessment of (anti)aromaticity with conventional magnetic, energetic, and geometric indices is difficult in the excited state, which practically restricts the extension and application of the concept of excited-state aromaticity. Time-resolved optical spectroscopies can provide a new and alternative avenue to evaluate excited-state aromaticity experimentally while observing changes in the molecular features in the excited states. Time-resolved optical spectroscopies take advantage of ultrafast laser pulses to achieve high time resolution, making them suitable for monitoring ultrafast changes in the excited states of molecular systems. This can provide valuable information for understanding the aromaticity reversal. This Account presents recent breakthroughs in the experimental assessment of excited-state aromaticity and the verification of aromaticity reversal with time-resolved optical spectroscopic measurements. To scrutinize this intriguing and challenging scientific issue, expanded porphyrins have been utilized as the ideal testing platform for investigating aromaticity because they show distinct aromatic and antiaromatic characters with aromaticity-specific spectroscopic features. Expanded porphyrins exhibit perfect aromatic and antiaromatic congener pairs having the same molecular framework but different numbers of π electrons, which facilitates the study of the pure effect of aromaticity by comparative analyses. On the basis of the characteristics of expanded porphyrins, time-resolved electronic and vibrational absorption spectroscopies capture the changes in electronic structure and molecular conformations driven by the change in aromaticity and provide clear evidence for aromaticity reversal in the excited states. The approaches described in this Account pave the way for the development of new and alternative experimental indices for the evaluation of excited-state aromaticity, which will enable overarching and fundamental comprehension of the role of (anti)aromaticity in the stability, dynamics, and reactivity in the excited states with possible implications for practical applications.

A study on the aromaticity and ring currents of dithienopyridines and dithienobenzene

The global aromaticity of dithienopyridine and dithienobenzene isomers was investigated using the topological resonance energy (TRE) and percentage topological resonance energy (%TRE) methods. The effect of variations in the positions of sulfur and nitrogen atoms on [Formula: see text]-electron delocalization is analyzed. The local aromaticity of these isomers is described based on the bond resonance energy (BRE) and circuit resonance energy (CRE) methods. Our BRE and CRE results show that structure of the central six-membered rings has a strong effect on global aromaticity. The aromaticity of these dithienopyridine isomers is enhanced when a complete pyridine unit exists in their middle ring structure, while the aromaticity of the dithienobenzene isomers is enhanced when a complete benzene unit exists in their middle ring structure. For dithienopyridines, our results obtained using the TRE method correlate well with the Bird aromaticity index as reported in the literature. Our ring-current results show that all these compounds are diatropic systems.

The predictive power of aromaticity: quantitative correlation between aromaticity and ionization potentials and HOMO–LUMO gaps in oligomers of benzene, pyrrole, furan, and thiophene

The aromaticity of oligomers predicts their HOMO–LUMO gaps and ionization potentials.

The Many Facets of Chalcogen Bonding: Described by Vibrational Spectroscopy

A diverse set of 100 chalcogen-bonded complexes comprising neutral, cationic, anionic, divalent, and double bonded chalcogens has been investigated using ωB97X-D/aug-cc-pVTZ to determine geometries, binding energies, electron and energy density distributions, difference density distributions, vibrational frequencies, local stretching force constants, and associated bond strength orders. The accuracy of ωB97X-D was accessed by CCSD(T)/aug-cc-pVTZ calculations of a subset of 12 complexes and by the CCSD(T)/aug-cc-pVTZ //ωB97X-D binding energies of 95 complexes. Most of the weak chalcogen bonds can be rationalized on the basis of electrostatic contributions, but as the bond becomes stronger, covalent contributions can assume a primary role in the strength and geometry of the complexes. Covalency in chalcogen bonds involves the charge transfer from a lone pair orbital of a Lewis base into the σ* orbital of a divalent chalcogen or a π* orbital of a double bonded chalcogen. We describe for the first time a symmetric chalcogen-bonded homodimer stabilized by a charge transfer from a lone pair orbital into a π* orbital. New polymeric materials based on chalcogen bonds should take advantage of the extra stabilization granted by multiple chalcogen bonds, as is shown for 1,2,5-telluradiazole dimers.

D2BIA-flexible, not (explicitly) arbitrary and reference/structurally invariant-a very effective and improved version of the D3BIA aromaticity index

Although there are a multitude of aromaticity indexes, only a few have a widespread usage. All famous aromaticity indexes are limited: HOMA and FLU are reference-dependent; ELF is π-bond-dependent; PDI is structurally dependent and NICS is ring size dependent. These limitations stimulate the continuous search for better (i.e., having no dependency), more flexible (i.e., applied to any aromatic system) and more effective (i.e., with excellent correlations with other indexes) aromaticity indexes. The D3BIA was our first topological aromaticity index. It is flexible, reference-independent and effective for planar and caged aromatic molecules. However, one of its terms, the degree of degeneracy (δ), is arbitrary and difficult to carry out for new users. Thus, in this work, we show that D2BIA-an improved version of D3BIA-is a good candidate to be used widely, since it retains the strong points of D3BIA while avoiding its weak point. In particular cases where all studied systems have δ = 1 (e.g., for acenes), then D2BIA equals D3BIA. For our recent study with acenes, D3BIA (and, as a consequence, D2BIA) has (have) an excellent correlation with FLU according to the MP3 method. In this work, by using DFT calculations for a series involving several six-membered and five-membered heteroaromatic rings, only D2BIA and NICS have very good correlation. All other well known aromaticity indexes used in this work (FLU, HOMA and ELF) gave poor correlations. As to homoaromatic systems, only D2BIA vs NICS and D2BIA vs FLU plots have excellent correlations. HOMA has the worst results in this series. Thus, D2BIA proved to be flexible and effective for the analysis of heteroaromatic rings of different sizes and for caged homoaromatic systems. Moreover, D2BIA has better correlations than D3BIA for planar aromatic systems, and same correlations for caged-homoaromatic systems. Graphical abstract D2BIA-an effective and improved version of the D3BIA aromaticity index.

The Peculiar Role of the Au3 Unit in Aum Clusters: σ-Aromaticity of the Au5Zn+ Ion

The stability of small Au_m (m = 4-7) clusters is investigated by analyzing their energetic, geometric, vibrational, magnetic, and electron density properties. Gold clusters can be constructed from stable cyclic 3-center-2-electron (3c-2e) Au₃⁺ units (3-rings) with σ-aromaticity. The stabilization requires a flow of negative charge from internal 3-rings with electron-deficient bonding to peripheral 3-ring units with stronger Au-Au bonds. The valence-isoelectronic clusters Au₆ and Au₅Zn⁺ have similar electronic properties: Au₅Zn⁺ is a strongly σ-aromatic molecule. An understanding of the structure of Au_m clusters is obtained by deriving a Clar's Rule equivalent for polycyclic gold clusters: The structure with the larger number of rings with dominant 3c-2e character and a smaller degree of 3c-3e character occupies the global minimum of the Au_m potential energy surface.

Optical spectra, electronic structure and aromaticity of benzannulated N-heterocyclic carbene and its analogues of the type C6H4(NR)2E: (E = Si, Ge, Sn, Pb)

Aromaticity of carbene analogues of the type C₆H₄(NR)₂E: (E = C–Pb) has been established by Raman, UV-vis, ISE, NICS and ACID methods.

Substituent-sensitive fluorescence of sequentially N-alkylated tetrabenzotetraaza[8]circulenes

We explore the use of substituent-sensitive balance between fluorescence and non-radiative decay as a tool for optical tuning of N-butylated tetrabenzotetraaza[8]circulenes for organic light emitting diode applications.