Aromaticity of benzenoid hydrocarbons with inserted –B=B– and –BH–BH– groups: a comparison

Structures of selected polycyclic conjugated hydrocarbons with –B=B– and –BH–BH– moieties inserted in different places were calculated at the B3LYP/6-311++G** level and their aromatic properties evaluated. HOMA, NICS(0), NICS(1)zz, Λ and PDI indices were used for studying their aromatic properties. Both optimized planar (as in parent hydrocarbons) and non-planar structures were taken into account. It is shown that insertion of both types of boron groups disturbs and decreases the aromaticity of the corresponding hydrocarbons. The decreasing effect of the –BH–BH– group is much stronger. What is quite intriguing is that it appears that non-planar structures of the studied compounds have a little higher aromaticity than the strictly planar ones. Mutual correlations between results obtained by different aromaticity indices are calculated and thoroughly discussed.     

HOMA parameters for the boron–boron bond: How the introduction of a BB bond influences the aromaticity of selected hydrocarbons

An extension of the harmonic oscillator model of aromaticity (HOMA) model for systems with boron–boron bonds is presented. For the first time, the parameters of the HOMA model are estimated using only theoretically calculated bond lengths. The HOMA parameters obtained make geometric aromaticity studies possible for a large number of compounds containing the boron–boron bond. The derived HOMA parameters have been used to investigate how the introduction of the boron–boron moiety in the structure of selected hydrocarbons modifies their aromaticity. The conclusion is that the insertion of a boron–boron bond usually strongly decreases the aromaticity of the boron-containing compounds in comparison to their parent hydrocarbons.     

A DFT study on the magnetostructural property of ferromagnetic heteroverdazyl diradicals with phenylene coupler

We have theoretically designed five different m-phenylene coupled high-spin bis-heteroverdazyl diradicals and their analogous p-phenylene coupled low-spin positional isomers. The geometry-based aromaticity index, harmonic oscillator model of aromaticity (HOMA) values for both the couplers (local HOMA), and the whole diradicals (global HOMA) have been calculated for all the diradicals. We also qualitatively relate these HOMA values with the intramolecular magnetic exchange coupling constants (J), calculated using a broken symmetry approach within unrestricted density functional theory framework. Structural aromaticity index HOMA of linkage specific benzene rings in our designed diradical systems shows that the aromatic character depends on the planarity of the molecule and it controls the sign and magnitude of J. The predicted J values are explained on the basis of spin polarization maps, average dihedral angles, and magnetic orbitals. The effect of the spin leakage phenomenon on magnetic exchange coupling constant and that on HOMA values of certain phosphaverdazyl systems has been explicitly discussed. In addition, a similar comparison is made between the calculated exchange coupling constants and corresponding HOMA values. The main novelty of this work stands on the consideration of the aromatic behavior by means of the geometrical index HOMA. We also estimate another aromaticity index, nucleus independent chemical shift (NICS) values for the phenylene coupler in each diradical to measure aromaticity and compare its value with that of HOMA. The ground state stabilities of these diradicals have also been compared.     

Tuning the Structure and Properties of N-doped Positively Charged Polycyclic Aromatic Hydrocarbons

A detailed study of the geometry, aromatic character, electronic and magnetic properties for a series of positively charged N-doped polycyclic aromatic hydrocarbons (PAHs) was performed. Magnetic properties of the examined molecules were analyzed by means of the magnetically induced current density calculated using the diamagnetic-zero version of the continuous transformation of origin of current density (CTOCD-DZ) method. The comparative study of the local aromaticity of the studied molecules was performed using several different indices: energy effect (ef), harmonic oscillator model of aromaticity (HOMA) index, six centre delocalization index (SCI) and nucleus independent chemical shifts (NICS). The presence of N-atoms in the inner rings was found to cause a planarity distortion in the studied N-doped systems. The geometric changes and charged nature of the studied N-doped systems do not significantly influence the current density and the local aromaticity distribution in comparison with the corresponding parent benzenoid hydrocarbons. The present study demonstrates how quantum chemical calculations can be used for rational design of novel PAHs and for fine tuning of their properties.     

Aromaticity and magnetotropicity of dicyclopenta-fused polyacenes

Most dicyclopenta-fused polyacenes are predicted to be moderately aromatic and diatropic, although they have no (4n + 2)-site conjugated circuits. We applied our graph theory of aromaticity and magnetotropicity to these molecules and found that these anomalous properties arise from a set of non-conjugated circuits, which contribute collectively to aromaticity and diatropicity. This result indicates that the conjugated circuit model is not always applicable to such non-alternant hydrocarbons. Dianions of dicyclopenta-fused polyacenes are more aromatic than their respective neutral species, because they are iso-pi-electronic with aromatic polyacenes.     

Five-membered heterocycles. Part III. Aromaticity of 1,3-imidazole in 5+n hetero-bicyclic molecules

The aromaticity (in the form of HOMA index) of 1,3-imidazole ring and its bicyclic derivatives was studied on the basis of statistical data from Cambridge Structural Database and X-ray investigations performed by authors. As a starting point, aromaticity of the 1,3-imidazoles with exocyclic X substituent at C2 (XN, O or S) was calculated. Subsequently, the HOMA index was calculated for various 5+n bicyclic skeletons with N, O, and S as endocyclic heteroatoms in the second ring. For the isolated 1,3-imidazole ring, aromaticity depends on exocyclic substituent at C2 and decreases in sequence N, S, O. For bicyclic derivatives it was found that both rings are aromatic and coplanar in very few molecules (17% of investigated ones), and positive charge—located on endocyclic heteroatom—increases the aromaticity. For remaining compounds, presence of sp³ carbons excludes possibility of aromatic ring existence.     

An insight into the local aromaticities of polycyclic aromatic hydrocarbons and fullerenes

In this work we quantify the local aromaticity of six-membered rings in a series of planar and bowl-shaped polycyclic aromatic hydrocarbons (PAHs) and fullerenes. The evaluation of local aromaticity has been carried out through the use of structurally (HOMA) and magnetically (NICS) based measures, as well as by the use of a new electronically based indicator of aromaticity, the para delocalization index (PDI), which is defined as the average of all the Bader delocalization indices between para-related carbon atoms in six-membered rings. The series of PAHs selected includes C(10)H(8), C(12)H(8), C(14)H(8), C(20)H(10), C(26)H(12), and C(30)H(12), with benzene and C(60) taken as references. The change in the local aromaticity of the six-membered rings on going from benzene to C(60) is analyzed. Finally, we also compare the aromaticity of C(60) with that of C(70), open [5,6]- and closed [6,6]-C(60)NH systems, and C(60)F(18).     

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

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

Aromaticity of the most stable adenine and purine tautomers in terms of Hückel's 4N+2 principle

Electron structures of the fused rings in 7H and 9H tautomers of purine and adenine follow the 4N + 2 rule; the values of pEDA, HOMA, NICS and FLU indices document their aromatic character. In the 1H and 3H tautomers, these rings contain five or seven π electrons, hence they do not follow this rule and consequently exhibit lower aromaticity. This also applies to the aromaticity of whole molecules.     

Density functional theoretical investigation of the aromatic nature of BN substituted benzene and four ring polyaromatic hydrocarbons

We have studied the topological and local aromaticity of BN-substituted benzene, pyrene, chrysene, triphenylene and tetracene molecules. The nucleus-independent chemical shielding (NICS), harmonic oscillator model of aromaticity (HOMA), para-delocalization index (PDI) and aromatic fluctuation index (FLU) have been calculated to quantify aromaticity in terms of magnetic and structural criteria. We find that charge separations due to the introduction of heteroatoms largely affect both the local and topological aromaticity of these molecules. Our studies show that the presence of any kind of heteroatom in the ring not only reduces the local delocalization in the six membered ring, but also affects strongly the topological aromaticity. In fact, the relative orders of the topological and local aromaticity depend strongly on the position of the heteroatoms in the structure. In general, more ring shared BN containing molecules are less aromatic than the less ring shared BN molecules. In addition our results provide evidence that the structural stability of the molecule is dominated by the σ bond rather than the π bond.     

On the stability of methyl cations with aromatic substituents. Theoretical approach

The extreme stability of methyl cations surrounded by certain aromatic fragments was studied using density functional theory methods. The pK_R+, aromaticity, ionization potential and dipole moment of tris[6-(dimethylamino)-1-azulenyl]methyl [1] and 2,6,10-tris(dimethylamino)trioxatriangulenium [2], as well as the corresponding for crystal violet [3] and tert-butanol [4], these last as the reference precursors of stable carbenium, were calculated and the results were found to be in agreement with experimental results previously reported. Furthermore, the calculations suggest that the large stability is a direct consequence of inductive effects as well as the intrinsic aromaticity of the branches as measured by the HOMA (harmonic oscillator model of aromaticity) method.     

Local aromaticity in benzo- and benzocyclobutadieno-annelated anthracenes

Abstract  

Considerations based on the energetics of cyclic conjugation in individual rings indicate that benzocyclobutadieno-annelation has the opposite effect on local aromaticity in benzenoid hydrocarbons to benzo-annelation. This finding is now tested and corroborated by density functional theory (DFT) calculations of the geometry of all benzo- and benzocyclobutadieno-annelated congeners of anthracene. The harmonic oscillator model of aromaticity (HOMA) and some similar (geometry-based) indices of local aromaticity are found to have the same dependence on the modes of annelation as the molecular-graph-based energy effects.     

Aromaticity and Diatropicity of Polyacenequinododimethides

Typical polyacenequinododimethides exist only in a single classical structure. These hydrocarbons are moderately aromatic and diatropic, although they have no aromatic conjugated circuits. This apparent dichotomy was resolved with our graph theory of aromaticity and magnetotropicity. Many nonconjugated circuits were found to contribute collectively to aromaticity and diatropicity. For individual molecules, local aromaticity increases with distance from the exo‐methylene groups. This fact indicates that the conjugated‐circuit model is not always applicable to semibenzenoid hydrocarbons such as polyacenequinododimethides.     

Local aromaticity of the six-membered rings in pyracylene. A difficult case for the NICS indicator of aromaticity

In this work, we have analyzed the local aromaticity of the six-membered rings (6-MRs) of planar and pyramidalized pyracylene species through the structurally based harmonic oscillator model of aromaticity (HOMA), the electronically based para-delocalization index (PDI), and the magnetic-based nucleus independent chemical shift (NICS) measurements, as well as with maps of ring current density. According to ring currents and PDI and HOMA indicators of aromaticity, there is a small reduction of local aromaticity in the 6-MRs of pyracylene with a bending of the molecule. In the case of NICS, the results depend on whether the NICS value is calculated at the center of the ring (NICS(0)) or at 1 A above (NICS(1)(out)) or below (NICS(1)(in)) the ring plane. While NICS(1)(out) values also indicate a slight decrease of aromaticity with bending, NICS(0) and NICS(1)(in) wrongly point out a large increase of aromaticity upon distortion. We have demonstrated that the NICS(0) reduction in the 6-MRs of pyracylene upon bending is due to (a) a strong reduction of the paratropic currents in 5-MRs and (b) the fact that, due to the distortion, the paratropic currents point their effects in other directions.     

The influence of aza-substitution on azole aromaticity

A homogeneous set of values for the aromaticity indices ASE (Aromatic Stabilisation Energy), HOMA (Harmonic Oscillator Model of Aromaticity) and NICS(1) (Nucleus-Independent Chemical Shift) for azoles has been investigated using multiple linear regression analysis. Statistically-significant relationships were found between the aromaticity indices and the number of nitrogen atoms at positions 2/5 and 3/4 of the ring. Aza-derivatives of pyrrole, furan and thiophene all gave similar relationships. For all three indices aza-substitution at positions 2 and/or 5 increases aromaticity. However, aza-substitution at positions 3 and/or 4 decreases classical aromaticity (ASE and HOMA) but increases magnetic aromaticity (NICS(1)). These indices appear to be measuring different properties of the azoles. The influence of aza-substitution on these different aspects of aromaticity is tentatively rationalized in terms of either bond length equalization or uniformity of π electron distribution.     

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.     

Analysis of electron delocalization in aromatic systems: individual molecular orbital contributions to para-delocalization indexes (PDI)

Our research group has recently defined two new aromaticity indexes based on the analysis of electron delocalization in aromatic species using the quantum theory of atoms-in-molecules. One of these indexes is the para-delocalization index (PDI) that measures the electronic delocalization between para-related carbon atoms in six-membered rings. In this paper, we show that this index can be partitioned into individual molecular orbital contributions. We have applied this PDI decomposition to several polycyclic aromatic hydrocarbons showing that this partitioning provides new insight into the origin of aromaticity.