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.     

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⁶⁻ anion agreed well with the reported experimental NMR values for He@C60 and He@C60⁶⁻.     

Structure, reactivity and aromaticity of acenes and their BN analogues: a density functional and electrostatic investigation

Density functional calculations have been carried out on a series of linearly annelated acenes and their BN analogues. Even though borazine shows aromatic and reactivity behavior parallel with that of benzene, its condensed derivatives show patterns different from those of their hydrocarbon analogues. Nucleus independent chemical shift (NICS) values in acenes suggest that the aromaticity of the inner rings is more than that of benzene, whereas in BN-acenes there is no substantial change in the aromaticity of the individual rings. Molecular electrostatic potential (MESP) is employed to obtain further insights into the bonding and reactivity trends for these systems. The MESP topography patterns of acenes and BN-acenes are substantially different, with BN-acenes showing more localized pi electron features compared to those of acenes. The MESP values at the critical points (CPs) indicate overall lowering of aromaticity in these annelated systems. However, this change is gradual among the BN-acenes.     

Is Free Cyclooctatetraene Dianion an Aromatic System? A Quantum Chemistry Study

To investigate whether free cyclooctatetraene dianion (COT^2?) is aromatic, quantum chemistry methods were used to optimize its structure. Based on the optimized structures, the natural population analysis (NPA) charge, bond order, delocalization energy, nucleus‐independent chemical shift (NICS), and harmonic oscillator model of aromaticity (HOMA) values were computed by DFT‐B3LYP method with basis set 6‐311++G**, which shows that COT^2? is not aromatic as it is not planar and has different bond lengths and bond orders, smallest delocalization energy and positive NICS values. To further confirm the finding, the changes of NICS and energy against ring distortion angle were scanned. The COT^2? has positive NICS values all along the angle from 180° to 120° while other aromatic systems always have negative values. The energy scanning suggests that COT^2? should have the weakest capability to maintain its planar structure. All the calculations strongly indicate that COT^2? is not aromatic. This study also suggests that NICS scan might be a good approach to judge aromaticity.     

Multidimensionality of delocalization indices and nucleus independent chemical shifts in polycyclic aromatic hydrocarbons

The aromaticity and local-aromaticity of a large set of polycyclic aromatic hydrocarbons (PAHs) is studied using multicenter delocalization indices from generalized population analysis and the popular nucleus independent chemical shift (NICS) index. A method for the fast computation of the NICS values is introduced, using the so-called pseudo-pi-method. A detailed examination is made of the multidimensional nature of aromaticity. The lack of a good correlation between the NICS and the multicenter delocalization indices is reported and the grounds discussed. It is shown through a thorough statistical analysis that the NICS values arise not only from local aromaticity of the benzenoid rings, but also from other circuits. It is shown that the NICS indices do not reveal the individual aromatic nature of a specific ring, contrary to the delocalization indices.     

The theoretical study of aromaticity in N-heteroatom compounds

In this article, we used nucleus-independent chemical shifts (NICS), aromatic stabilization energies (ASE), and magnetic susceptibility isotropic, calculated with (density functional theory) B3LYP levels at the 6-311+G (d, p) basis set, to evaluate the aromaticity of a set of 12 six-member planar π-electron aromatic systems: mono- and multi-substituted benzenes by nitrogen atoms. NICS were calculated at the center of the rings, 0.5 Å and 1 Å above the molecular plane, and in order to reflect the π-electron effects, we compare the values of NICS(0.5). The result statistically revealed significant correlations among the above three criteria. The comparison result was obtained as benzene > mono-nitrogen benzene > di-nitrogen benzene > tri-nitrogen benzene > tetra-nitrogen benzene > penta-nitrogen benzene.     

Magnetic properties and aromaticity of o-, m-, and p-benzyne

The relative aromaticities of the three singlet benzyne isomers, 1,2-, 1,3-, and 1,4-didehydrobenzenes have been evaluated with a series of aromaticity indicators, including magnetic susceptibility anisotropies and exaltations, nucleus-independent chemical shifts (NICS), and aromatic stabilization energies (all evaluated at the DFT level), as well as valence-bond Pauling resonance energies. Most of the criteria point to the o-benzyne相似文献   

Application of AIM parameters at ring critical points for estimation of pi-electron delocalization in six-membered aromatic and quasi-aromatic rings

The AIM parameters at the ring critical point (the electron density and its Laplacian, the total electron energy density and both its components, potential and kinetic electron energy densities), have been intercorrelated with aromaticity indices: the geometry-based HOMA and the magnetism-based NICS, NICS(1), and NICS(1)(zz). A set of 33 phenylic rings having possibly a diversified aromatic character, and a set of 20 quasi-rings formed by intramolecular hydrogen and lithium bonds, have been taken into consideration. It has been found that the density of total electron energy, H, may serve as a new quantitative characteristic of pi-electron delocalization. The dependences between H values and aromaticity indices are correlated (cc(H/HOMA)=0.99, cc(H/NICS(1)zz)=0.95).     

Which NICS aromaticity index for planar pi rings is best?

[STRUCTURE: SEE TEXT] Five increasingly sophisticated aromaticity indexes, based on nucleus-independent chemical shifts (NICS), were evaluated against a uniform set of aromatic stabilization energies (ASE) for 75 mono- and polyheterocyclic five-membered rings. While acceptable statistical correlations were given by all of the NICS methods, the most fundamentally grounded index, NICS(0)pizz (based on the pi contribution to the out-of-plane zz tensor component), performed best statistically (cc=0.980) and in practice. The easily computable NICS(1)zz index is a useful alternative (cc=0.968).     

An insight on the aromatic changes in closed shell icosagen,tetrel, and pnictogen phenalenyl derivatives

A computational study of the aromatic and antiaromatic characteristics of closed shell charged phenalenyl (PLY⁺¹ and PLY^?1) upon replacement of the central carbon atom by icosagen (B, Al and Ga), tetrel (Si and Ge) and pnictogen (N, P and As) atoms comprising systems in which the icosagen and pnictogen derivatives considered are neutral while the tetrel ones are anions or cations, has been carried out at the B3LYP/6–311++G(d,p) computational level. By substitution, two different kinds of structures have been obtained, one planar (N and B) and another one bowl-shaped depending on the size of the central atom. In terms of aromaticity, the substitution of the central C atom causes a loss of the aromatic character in all cases as indicated by nucleus-independent chemical shifts (NICS) profiles and NICS values on the 0.001 au isosurface. Regarding the charge, PLY⁺¹ presents larger electron delocalisation than PLY^?1, phenomenon associated with aromaticity. Furthermore, the current density maps for those planar systems corroborate NICS findings, showing anticlockwise currents in PLY⁺¹ (like in benzene) but clockwise in PLY-N⁰ and PLY-B⁰, indicating aromatic and antiaromatic behaviour, respectively.     

Local aromaticity in polyacenes manifested by individual proton and carbon shieldings: DFT mapping of aromaticity

Exponential dependencies between locally calculated geometric and magnetic indexes of aromaticity, harmonic oscillator model of aromaticity (HOMA) and nucleus independent chemical shifts (NICS)(0), NICS(1) and NICS(1)zz, and the number of conjugated benzene rings in linear acenes, from benzene to decacene were observed at B3LYP/6-311+G** level of theory. Correlations between HOMA and NICS indexes showed exponential dependencies and were fitted with simple three-parameter function. Similar correlations between both indexes of aromaticity and proton and carbon nuclear isotropic shieldings of individual acene rings were observed. Contrary to proton data, the predicted ¹³C nuclear isotropic shieldings of carbon atoms belonging to inner rings in polyacenes were less shielded, indicating lower aromaticity and therefore, higher reactivity.     

Local aromaticities in large polyacene molecules

There has been controversy on the relative aromaticities of individual rings in a large polyacene molecule. Nucleus-independent chemical shift (NICS) values suggest that the highly reactive inner rings might be more aromatic than the outer ones and even more aromatic than benzene. We evaluated the bond resonance energies (BREs) and hypothetical geometry-independent pi-electron currents for a series of linear polyacenes and noticed that for large polyacene molecules the inner rings are never more aromatic than the outer ones. Global HOMA (harmonic oscillator model of aromaticity) values are highly correlative with percentage topological resonance energies (% TREs) but not with average NICS values. Magnetic properties, such as NICS and ring-current intensity, are highly dependent on molecular geometry and so must be carefully related to aromaticity.     

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.     

运用密度泛函活性理论的信息论方法研究苯并富烯衍生物的芳香性

文献中已有越来越多的芳香性体系被发现,同时也有越来越多的芳香性指标被提出来,但是如何解释芳香化合物稳定性的起源以及理解芳香性的本质仍然是当今理论化学中一个悬而未决的难题。运用我们新近提出的密度泛函活性理论信息论方法,不久前我们曾对一系列富烯衍生物进行了系统研究并得到了一个全新的认识。本文进一步探讨苯并富烯衍生物的芳香性行为,目的在于考察一个或多个苯环与富烯连接之后其芳香性发生变化的情况。运用香农熵,费舍尔信息,Ghosh-Berkowitz-Parr熵,Onicescu信息能,信息增益,以及相对Rényi熵六个信息量,和四种芳香指标,ASE,HOMA,FLU和NICS,我们系统地研究了信息量和芳香性指标在单、双、三苯并富烯衍生物中的相关性。我们发现,不管是否有苯环与富烯相连,芳香指标和信息量的交叉相关性都是一样的。这表明,虽然苯环本身具有芳香性,但苯环与富烯相连并不能改变富烯的芳香性与反芳香性本质。苯并富烯衍生物与富烯衍生物的芳香性和反芳香性一致。苯并富烯衍生物的芳香性和反芳香性完全取决于富烯本身的芳香性和反芳香性。这些结果为认识和理解复杂体系芳香性和反芳香性起源和本质将提供有益的启示。     

On the correlation between the bifurcation values of the electron localization function and the nuclear independent chemical shift in aromatic and heteroaromatic compounds

In this work, for a representative set of 24 aromatic molecules, that includes hydrocarbons, hetero and charged rings, we explore the correlations between two aromaticity indexes. One of these indexes is based on the bifurcation values of the σ and π electron localization function (ELF), while the other one is the nuclear independent chemical shift (NICS). We observe that the quality of the possible correlations between these two kinds of indexes strongly depends on the kind of rings and on the particular indexes that are considered.     

The delocalization index as an electronic aromaticity criterion: application to a series of planar polycyclic aromatic hydrocarbons

This work introduces a new local aromaticity measure, defined as the mean of Bader's electron delocalization index (DI) of para-related carbon atoms in six-membered rings. This new electronic criterion of aromaticity is based on the fact that aromaticity is related to the cyclic delocalized distribution of pi-electrons. We have found that this DI and the harmonic oscillator model of aromaticity (HOMA) index are strongly correlated for a series of six-membered rings in eleven planar polycyclic aromatic hydrocarbons. The correlation between the DI and the nucleus-independent chemical shift (NICS) values is less remarkable, although in general six-membered rings with larger DI values also have more negative NICS indices. We have shown that this index can also be applied, with some modifications, to study of the aromaticity in five-membered rings.     

Ligand-stabilized aromatic three-membered gold rings and their sandwichlike complexes

Electronic structure calculations (DFT) suggest that ligand-stabilized three-membered gold(I) rings constituting the core structure in a series of cyclo-Au3L(n)H(3-n) (L = CH3, NH2, OH and Cl; n = 1, 2, 3) molecules exhibit aromaticity, which is primarily due to 6s and 5d cyclic electron delocalization over the triangular Au3 framework (s- and d-orbital aromaticity). The aromaticity of the novel triangular gold(I) isocycles was verified by a number of established criteria of aromaticity. In particular, the nucleus-independent chemical shift, NICS(0), the upfield changes in the chemical shifts for Li+, Ag+, and Tl+ cations over the Au3 ring plane, and their interaction with electrophiles (e.g., H+, Li+, Ag+, and Tl+) are indicative for the aromaticity of the three-membered gold(I) rings. Interestingly, unlike the respective substituted derivatives of cyclopropenium cation and the bora-cyclopropene carbacyclic analogues, the aromatic Au3 rings, although exhibit comparable diatropicity, react with electrophiles in a different way affording 1:1 and 2:1 sandwichlike complexes. The bonding in the three-membered gold(I) rings is characterized by a common ring-shaped electron density, more commonly seen in aromatic organic molecules and in "all-metal" aromatics, such as the cyclo-[Hg3]4- tetraanion. Moreover, the cation-pi interactions in the 1:1 and 1:2 sandwichlike complexes formed upon reacting the Au3 rings with electrophiles, depending on the nature of the cation, are predicted to be predominantly electrostatic (Li+, Tl+) or covalent (H+, Ag+). The 1:2 complexes constitute a new class of sandwichlike complexes, which are expected to have novel properties and applications.     

Sigma-aromaticity and sigma-antiaromaticity in saturated inorganic rings

The magnetic properties of a series of inorganic saturated rings, (SiH2)n, (GeH2)n, (NH)n, (PH)n, (AsH)n, On, Sn, and Sen (n = 3-6), exhibit zigzag behavior with ring size resembling that of aromatic and antiaromatic Hückel pi-systems and (CH2)n rings. Computed GIAO-SCF nucleus-independent chemical shifts (NICS) and localized (LMO) NICS analysis indicate that the sigma-ring electrons are chiefly responsible for this zigzag behavior. This evidence for sigma-aromaticity is further supported by theoretical strain energy (TSE). The Hückel 4n + 2/4n aromaticity/antiaromaticity rule for pi-electron systems applies well to the smaller saturated rings.