Aromaticity changes along the reaction coordinate connecting the cyclobutadiene dimer to cubane and the benzene dimer to hexaprismane

Aromaticity and reactivity are two deeply connected concepts. Most of the thermally allowed cycloadditions take place through aromatic transition states, while transition states of thermally forbidden reactions are usually less aromatic, if at all. In this work, we perform a numerical experiment to discuss the change of aromaticity that occurs along the reaction paths that connect two antiaromatic units of cyclobutadiene to form cubane and two aromatic rings of benzene to yield hexaprismane. It is found that the aromaticity profile along the reaction coordinate of the [4+4] cycloaddition of two antiaromatic cyclobutadiene molecules goes through an aromatic highest energy point and finishes to an antiaromatic cubane species. Up to our knowledge, this represents the first example of a theoretically and thermally forbidden reaction path that goes through an intermediate aromatic region. In contrast, the aromaticity profile in the [6+6] cycloaddition of two aromatic benzene rings show a slow steady decrease of aromaticity from reactants to the highest energy point and from this to the final hexaprismane molecule a plunge of aromaticity is observed. In both systems, the main change of aromaticity occurs abruptly near the highest energy point, when the distance between the centers of the two rings is about 2.2 Å.     

Interplay among aromaticity, magnetism, and nonlinear optical response in all-metal aromatic systems

All-metal aromatic molecules are the latest inclusion in the family of aromatic systems. Two different classes of all-metal aromatic clusters are primarily identified: one is aromatic only in the low spin state, and the other shows aromaticity even in high-spin situations. This observation prompts us to investigate the effect of spin multiplicity on aromaticity, taking Al(4)(2-), Te(2)As(2)(2-), and their copper complexes as reference systems. Among these clusters, it has been found that the molecules that are aromatic only in their singlet state manifest antiaromaticity in their triplet state. The aromaticity in the singlet state is characterized by the diatropic ring current circulated through the bonds, which are cleaved to generate excess spin density on the atoms in the antiaromatic triplet state. Hence, in such systems, an antagonistic relationship between aromaticity and high-spin situations emerges. On the other hand, in the case of triplet aromatic molecules, the magnetic orbitals and the orbitals maintaining aromaticity are different; hence, aromaticity is not depleted in the high-spin state. The nonlinear optical (NLO) behavior of the same set of clusters in different spin states has also been addressed. We correlate the second hyperpolarizability and spin density in order to judge the effect of spin multiplicity on third-order NLO response. This correlation reveals a high degree of NLO behavior in systems with excess spin density. The variance of aromaticity and NLO response with spin multiplicity is found to stem from a single aspect, the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), and eventually the interplay among aromaticity, magnetism, and NLO response in such materials is established. Hence, the HOMO-LUMO energy gap becomes the cornerstone for tuning the interplay. This correlation among the said properties is not system-specific and thus can be envisaged even beyond the periphery of all-metal aromatic clusters. Such interplay is of crucial importance in tailoring novel paradigm of multifunctional materials.     

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.     

The aromatic fluctuation index (FLU): a new aromaticity index based on electron delocalization

In this work, the aromatic fluctuation index (FLU) that describes the fluctuation of electronic charge between adjacent atoms in a given ring is introduced as a new aromaticity measure. This new electronic criterion of aromaticity is based on the fact that aromaticity is related to the cyclic delocalized circulation of pi electrons. It is defined not only considering the amount of electron sharing between contiguous atoms, which should be substantial in aromatic molecules, but also taking into account the similarity of electron sharing between adjacent atoms. For a series of rings in 15 planar polycyclic aromatic hydrocarbons, we have found that, in general, FLU is strongly correlated with other widely used indicators of local aromaticity, such as the harmonic-oscillator model of aromaticity, the nucleus independent chemical shift, and the para-delocalization index (PDI). In contrast to PDI, the FLU index can be applied to study the aromaticity of rings with any number of members and it can be used to analyze both the local and global aromatic character of rings and molecules.     

Solvent-dependent aromatic versus antiaromatic conformational switching in meso-(heptakis)pentafluorophenyl [32]heptaphyrin

We have investigated conformational switching dynamics of meso-heptakis(pentafluorophenyl) [32]heptaphyrin(1.1.1.1.1.1.1) in various solvents using steady-state, time-resolved, and temperature-dependent spectroscopy. Absorption and fluorescence spectra of [32]heptaphyrin are quite sensitive to solvent environments. In nonpolar toluene, the antiaromatic figure-of-eight conformation is dominant, as seen in the X-ray crystallography, based on broad and weak absorption bands without any fluorescence and moderate paratropic ring current. On the other hand, a well-resolved sharp absorption spectrum, strong fluorescence, and diatropic ring current in the 1H NMR spectrum in slightly polar THF indicate that most of [32]heptaphyrin molecules take significantly distorted M?bius conformation with aromatic character. By using transient absorption (TA) spectroscopy, the lowest singlet excited-state lifetimes have been revealed to decay biexponentially with the time constants of 5 and 65 ps for figure-of-eight and M?bius conformations, respectively. Based on these results along with vertical excitation energy calculations, we are able to assign two conformers as Hückel antiaromatic and M?bius aromatic species, respectively; it shoulf be noted that the aromaticity/antiaromaticity does not change with temperature variation. Interestingly, in moderately polar solvent, ethyl ether, we find out that these two conformational isomers coexist with a dynamic equilibrium, as revealed by excitation-wavelength-dependent TA, temperature-dependent absorption and 1H NMR spectra. Through our findings, we have demonstrated that the conformational switching dynamics between Hckel antiaromatic and M?bius aromatic conformers in [32]heptaphyrin(1.1.1.1.1.1.1) are strongly affected by solvent medium as well as temperature.     

Peculiar antiaromatic inorganic molecules of tetrapnictogen in Na+Pn4- (Pn = P,As, Sb) and important consequences for hydrocarbons

Although aromaticity has been observed in inorganic and all-metal species, the concept of antiaromaticity has not been extended beyond organic molecules. Here, we present theoretical and experimental evidence that the 6 -electron tetrapnictogen dianions in Na+Pn42- (Pn = P, As, Sb) undergo a transition from being aromatic to antiaromatic upon electron detachment, yielding the first inorganic antiaromatic Na+Pn4- molecules. Two types of antiaromatic structures were characterized, the conventional rectangular species and a new peculiar quasiplanar rhombus species. Aromaticity and antiaromaticity in the tetrapnictogen molecules were derived from molecular orbital analyses and verified by experimental photodetachment spectra of Na+Pn42-. On the basis of our findings for the tetrapnictogen clusters, we predicted computationally that the organic C4H4- anion also possesses two antiaromatic structures: rectangular and rhombus. Moreover, only the rhombus antiaromatic minimum was found for the radical NC3H4, thus extending the peculiar rhombus antiaromatic structure first uncovered in inorganic clusters into organic chemistry.     

Ground- and excited-state aromaticity and antiaromaticity in benzene and cyclobutadiene

The aromaticity and antiaromaticity of the ground state (S 0), lowest triplet state (T 1), and first singlet excited state (S 1) of benzene, and the ground states (S 0), lowest triplet states (T 1), and the first and second singlet excited states (S 1 and S 2) of square and rectangular cyclobutadiene are assessed using various magnetic criteria including nucleus-independent chemical shifts (NICS), proton shieldings, and magnetic susceptibilities calculated using complete-active-space self-consistent field (CASSCF) wave functions constructed from gauge-including atomic orbitals (GIAOs). These magnetic criteria strongly suggest that, in contrast to the well-known aromaticity of the S 0 state of benzene, the T 1 and S 1 states of this molecule are antiaromatic. In square cyclobutadiene, which is shown to be considerably more antiaromatic than rectangular cyclobutadiene, the magnetic properties of the T 1 and S 1 states allow these to be classified as aromatic. According to the computed magnetic criteria, the T 1 state of rectangular cyclobutadiene is still aromatic, but the S 1 state is antiaromatic, just as the S 2 state of square cyclobutadiene; the S 2 state of rectangular cyclobutadiene is nonaromatic. The results demonstrate that the well-known "triplet aromaticity" of cyclic conjugated hydrocarbons represents a particular case of a broader concept of excited-state aromaticity and antiaromaticity. It is shown that while electronic excitation may lead to increased nuclear shieldings in certain low-lying electronic states, in general its main effect can be expected to be nuclear deshielding, which can be substantial for heavier nuclei.     

Can Substituted Cyclopentadiene Become Aromatic or Antiaromatic?

Cyclopentadiene derivatives with electronegative (F, Cl) or electropositive (H(3)Si, Me(3)Si) bis-5,5-substituents were studied at the B3LYP/6-311G* level of theory. It was found that there is no special stabilization or destabilization for any of the derivatives; the energetic effects that were previously attributed to aromatic stabilization or antiaromatic destabilization are the result of interactions in the reference systems. A nucleus-independent chemical shift (NICS) scan study at the HF-GIAO/6-311+G* theoretical level of these and similar derivatives suggest that they all show different magnitudes of diamagnetic ring current. None of the derivatives shows a paramagnetic ring current. Thus, cyclopentadienes are neither aromatic nor antiaromatic. It is also concluded that a diamagnetic ring current is perhaps necessary but certainly not a sufficient condition for aromaticity. The NICS scan procedure describes the type of ring current in the system, whereas a single isotropic NICS value (i.e., NICS(1)) may wrongly assign the type of ring current. It is shown that neither NICS(1) nor the NICS scan procedure can be used as a single aromaticity criterion.     

Unique Role of Heterole-Fused Structures in Aromaticity and Physicochemical Properties of 7,8-Dehydropurpurins

Porphyrins with a fused five-membered ring, such as 7,8-dehydropurpurins, have appeared as an emerging class of unique porphyrinoids. Their altered absorption spectra, relatively short lifetimes of excited states, and small HOMO–LUMO gaps arise from the harmony of the antiaromatic 20π-circuit and the aromatic 18π-circuit. In this regard, the electronic properties of 7,8-dehydropurpurins are expected to be controlled by modulating the contribution of the antiaromatic π-circuit to the whole aromaticity. Here the comparison of pyrrole- and phosphole-fused 7,8-dehydropurpurins is reported in terms of their aromaticity and physicochemical properties. The spectroscopic investigation revealed the larger contribution of the antiaromatic 24π-circuit in pyrrole-fused 7,8-dehyrdopurpurins than in phosphole-fused 7,8-dehydropurpurins. The DFT calculations also supported the feasibility of tuning the aromaticity of 7,8-dehydropurpurins by heterole-fused structures. Thus, the introduction of heterole-fused structures into porphyrinoids is a universal strategy to get new insight into aromaticity and their intrinsic properties in cyclic π-conjugated molecules.     

Dianion and dication of tetrabenzo[5.7]fulvalene. Greater antiaromaticity than aromaticity in comparable systems

The dianion, 5(2-), and dication, 5(2+), of tetrabenzo[5.7]fulvalene represent an aromaticity/antiaromaticity continuum in which the fluorenyl system changes from aromatic in 5(2-) to antiaromatic in 5(2+). Conversely, the antiaromatic dibenzotropylium system of 5(2-) becomes an aromatic system in 5(2+), allowing an examination of aromaticity/antiaromaticity within the same carbon framework. Dianion 5(2-) was prepared and characterized by (1)H NMR spectroscopy. The fluorenyl system of 5(2-) showed the downfield shifts expected for an aromatic system, while the dibenzotropylium system showed the paratropic shifts expected for an antiaromatic system. The conclusions from (1)H NMR spectroscopy were supported by NICS(1) zz calculations for each system. Comparison of the (1)H NMR spectrum and NICS(1) zz of 5(2-) with those of 5(2+) supported the assignments of aromaticity/antiaromaticity for each system. Aromaticity/antiaromaticity were further examined through comparison of the degree of bond length alternation, which showed that the bond length alternation was slightly greater for the antiaromatic ring systems than for the aromatic systems. However, when structures of 5(2-) and 5(2+) with no bond length alternation were examined, there was a dramatic increase in the degree of antiaromaticity for the antiaromatic ring systems as evaluated through NICS. This result suggests that a decrease in bond length alternation results in an increase in antiaromaticity as well as an increase in aromaticity. The magnitude of the antiaromaticity of the fluorenyl system in 5(2+) was greater than the magnitude of the aromaticity in the fluorenyl system of 5(2-), with similar effects shown by the analogous tropylium systems. This is consistent with the behavior of the antiaromatic dication of tetrabenzo[5.5]fulvalene, compared to that of its aromatic dianion, and also with the behavior of the cyclopentadienyl cation/anion and tropylium cation/anion.     

Gauge-origin-independent magnetizabilities of solvated molecules using the polarizable continuum model

We present an implementation of the polarizable continuum model in its integral equation formulation for the calculation of the magnetizabilities of solvated molecules. The gauge-origin independence of the calculated magnetizabilities and the fast basis set convergence are ensured through the use of London atomic orbitals. Our implementation can use Hartree-Fock and multiconfigurational self-consistent-field (MCSCF) wave functions as well as density-functional theory including hybrid functionals such as B3LYP. We present the results of dielectric continuum effects on water and pyridine using MCSCF wave functions, as well as dielectric medium effects on the magnetizability of the aromatic amino acids as a model for how a surrounding protein environment affects the magnetizability of these molecules. It is demonstrated that the dielectric medium effects on the magnetizability anisotropies of the aromatic amino acids may be substantial, being as large as 25% in the case of tyrosine.     

Origin of the Electron Transport Properties of Aromatic and Antiaromatic Single Molecule Circuits

Antiaromatic molecules have been predicted to exhibit increased electron transport properties when placed between two nanoelectrodes compared to their aromatic analogues. While some studies have demonstrated this relationship, others have found no substantial increase. We use atomistic simulations to establish a general relationship between the electronic spectra of aromatic, antiaromatic, and quinoidal molecules and illustrate its implications for electron transport. We compare the electronic properties of a series of aromatic-antiaromatic counterparts and show that antiaromaticity effectively p-dopes the aromatic electronic spectra. As a consequence, the conducting properties of aromatic-antiaromatic analogues are closely related. For similar attachment points to the electrodes, an interference feature is expected in the HOMO-LUMO gap of one whenever it is absent in the other one. We demonstrate how the relative conductance of aromatic-antiaromatic pairs can be tuned and even reversed through the choice of chemical linker groups. Our work provides a general picture relating connectivity, (anti)aromaticity, and quantum interference and establishes new design rules for single molecule circuits.     

To what extent can aromaticity be defined uniquely?

Statistical analyses of quantitative definitions of aromaticity, ASE (aromatic stabilization energies), RE (resonance energies), Lambda (magnetic susceptibility exaltation), NICS, HOMA, I5, and A(J), evaluated for a set of 75 five-membered pi-electron systems: aza and phospha derivatives of furan, thiophene, pyrrole, and phosphole (aromatic systems), and a set of 30 ring-monosubstituted compounds (aromatic, nonaromatic, and antiaromatic systems) revealed statistically significant correlations among the various aromaticity criteria, provided the whole set of compounds is involved. Hence, broadly considered, the various manifestations of aromaticity are related and aromaticity can be regarded statistically as a one-dimensional phenomenon. In contrast, when comparisons are restricted to some regions or groups of compounds, e.g., aromatic compounds with ASE > 5 kcal/mol or polyhetero-five-membered rings, the quality of the correlations can deteriorate or even vanish. In practical applications, energetic, geometric, and magnetic desriptors of aromaticity do not speak with the same voice. Thus, in this sense, the phenomenon of aromaticity is regarded as being statistically multidimensional.     

Expanded porphyrins and aromaticity

meso-Aryl-substituted expanded porphyrins that are porphyrin homologues consisting of more than five pyrrolic units are a nice platform to realize diverse aromatic and antiaromatic species as well as stable radical species. They are also an ideal series to realize topologically twisted molecules with distinct M?bius aromaticity and antiaromaticity.     

QTAIM N-center delocalization indices as descriptors of aromaticity in mono and poly heterocycles

The implementation of the n-center electron delocalization indices, n-DIs, and n-order electron localization indices, n-LIs, within the framework of the quantum theory of atoms in molecules, QTAIM, is performed. n-DIs are shown to be very useful to study the local aromaticity in monocyclic and polycyclic compounds. Total and pi n-DIs from n=4 to 7 were computed for a series of typical 4, 5, 6, and 7-center aromatic and antiaromatic rings. For n>or=5 the pi n-DI accounts for the 95% of the total n-DI and can be employed alone to measure the aromaticity. A scaling factor on the n-DIs is required in order to compare the aromaticity of [5c-6e] and [6c-6e] rings, the same correction allows to estimate the relative aromatic stabilization of polycyclic compounds using the sum of its values for individual rings. This is called Effective Scaled Electron Delocalization, ESED. The comparison with other aromaticity indices reflects a good correlation between ESED and both resonance energies, and HOMA indices. The most important differences between scaled pi n-DIs and NICS(0) indices are found for compounds that contain rings with different number of centers or pi electrons.     

Probing the Most Aromatic and Antiaromatic Pyrrolium Rings by Maximizing Hyperconjugation and Push–Pull Effect

Hyperconjugation, a weak interaction in organic chemistry, can have a strong effect on aromaticity, leading to the concept of hyperconjugative aromaticity, which was first proposed by Mulliken in 1939. However, most studies are limited to main group chemistry. Here we report the most aromatic and antiaromatic pyrrolium ring by maximizing the hyperconjugation caused by transition metal fragments and the push–pull effect. Our calculations reveal that the origin of the outperformance of transition metal substituents over main group ones on hyperconjugative aromaticity could be attributed to their higher highest occupied molecular orbitals (HOMOs). Among the group 11 transition metals, a silver substituent results in the best performance. All these findings highlight the magic of the transition metal (silver) and could be particularly helpful for the design of other aromatic and antiaromatic counterparts based on a nonaromatic parent species.     

Aromaticity Reversal in the Lowest Excited Triplet State of Archetypical Möbius Heteroannulenic Systems

The aromaticity reversal in the lowest triplet state (T₁) of a comparable set of Hückel/Möbius aromatic metalated expanded porphyrins was explored by optical spectroscopy and quantum calculations. In the absorption spectra, the T₁ states of the Möbius aromatic species showed broad, weak, and ill‐defined spectral features with small extinction coefficients, which is in line with typical antiaromatic expanded porphyrins. In combination with quantum calculations, these results indicate that the Möbius aromatic nature of the S₀ state is reversed to Möbius antiaromaticity in the T₁ state. This is the first experimental observation of aromaticity reversal in the T₁ state of Möbius aromatic molecules.     

Aromaticity: molecular-orbital picture of an intuitive concept

Geometry is one of the primary and most direct indicators of aromaticity and antiaromaticity: a regular structure with delocalized double bonds (e.g., benzene) is symptomatic of aromaticity, whereas a distorted geometry with localized double bonds (e.g., 1,3-cyclobutadiene) is characteristic of antiaromaticity. Here, we present a molecular-orbital (MO) model of aromaticity that explains, in terms of simple orbital-overlap arguments, why this is so. Our MO model is based on accurate Kohn-Sham DFT analyses of the bonding in benzene, 1,3-cyclobutadiene, cyclohexane, and cyclobutane, and how the bonding mechanism is affected if these molecules undergo geometrical deformations between regular, delocalized ring structures, and distorted ones with localized double bonds. We show that the propensity of the pi electrons is always, that is, in both the aromatic and antiaromatic molecules, to localize the double bonds, against the delocalizing force of the sigma electrons. More importantly, we show that the pi electrons nevertheless decide about the localization or delocalization of the double bonds. A key component of our model for uncovering and resolving this seemingly contradictory situation is to analyze the bonding in the various model systems in terms of two interpenetrating fragments that preserve, in good approximation, their geometry along the localization/delocalization modes.     

Visualizing electron delocalization in contorted polycyclic aromatic hydrocarbons

Electron delocalization in contorted polycyclic aromatic hydrocarbon (PAH) molecules was examined through 3D isotropic magnetic shielding (IMS) contour maps built around the molecules using pseudo-van der Waals surfaces. The resulting maps of electron delocalization provided an intuitive, yet detailed and quantitative evaluation of the aromatic, non aromatic, and antiaromatic character of the local and global conjugated cyclic circuits distributed over the molecules. An attractive pictural feature of the 3D IMS contour maps is that they are reminiscent of the Clar π-sextet model of aromaticity. The difference in delocalization patterns between the two faces of the electron circuits in contorted PAHs was clearly visualized. For π-extended contorted PAHs, some splits of the π system resulted in recognizable patterns typical of smaller PAHs. The differences between the delocalization patterns of diastereomeric chiral PAHs could also be visualized. Mapping IMS on pseudo-van der Waals surfaces around contorted PAHs allowed visualization of their superimposed preferred circuits for electron delocalization and hence their local and global aromaticity patterns.

Electron delocalization in contorted polycyclic aromatic hydrocarbon (PAH) molecules was examined through 3D isotropic magnetic shielding (IMS) contour maps built around the molecules using pseudo-van der Waals surfaces.