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.     

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.     

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.     

Electron delocalization and aromaticity in linear polyacenes: atoms in molecules multicenter delocalization index

Molecular aromaticity in the linear polyacenes is investigated using an atoms in molecules based six center index (SCI-AIM) which measures the electron delocalization. SCI-AIM values for the linear polyacenes indicate decreasing aromaticity going from outer to inner rings in the polyacene series. The SCI-AIM approach is compared to a Mulliken-like approach, and a critical comparison to the PDI index is made.     

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).     

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.     

Interatomic magnetizability: a QTAIM-based approach toward deciphering magnetic aromaticity

Interatomic magnetizability provides insight into the extent of electronic current density between two adjacent atomic basins. By studying a number of well-known aromatic, nonaromatic, and antiaromatic molecules, it is demonstrated that interatomic magnetizability (bond magnetizability) not only is able to verify the exact nature of aromaticity/antiaromaticity among different molecules, but also can distinguish the correct aromaticity order among sets of aromatic/antiaromatic molecules. The interatomic magnetizability is a direct measure of the current flux between two adjacent atomic basins and is the first QTAIM-derived index that evaluates aromaticity based on a response property, that is, magnetizability. Bond magnetizability is easy to compute, straightforward to interpret, and can be employed to evaluate the pure π- or σ-orbital contributions to magnetic aromaticity.     

Properties of aromaticity indices based on the one-electron density matrix

Proper normalization of two previously published indices yields aromaticity measures that, when computed within the Hückel molecular orbital (HMO) approximation, closely match the topological resonance energies per pi electron of aromatic annulenes and their ions. The normalized indices, which quantify aromaticity of individual rings in polycyclic systems, are equally applicable to homocyclic and heterocyclic compounds and can be readily computed from 1-matrices calculated at any level of electronic structure theory. However, only the index ING, derived from the Giambiagi formula, produces proper ordering of aromaticities of heterocyclic compounds, provided it is calculated from all-electron wavefunctions in conjunction with the atoms in molecule (AIM) partitioning. Its values are shown to be strongly affected by electron correlation effects. Because of its apparent inability to distinguish between anti- and nonaromatic systems, ING should only be employed for aromatic species.     

Electron delocalization and aromaticity variations in the stacked nucleic acid base pairs

A newly developed exchange-correlation functional (MPWB1K) in density functional theory has been applied to evaluate the electron delocalization of individual fragments in the stacking interaction between nucleic acid bases (NABs). Electronically and structural-based indices have been employed to investigate the aromaticity variation during stacking interaction. A quantitative study of NABs in their isolated and stacked forms reveals that stacking interaction causes a decrease in bond delocalization. It is shown that the decrease in the aromaticity is accompanied by local decrease in two-center delocalization indices within the pyrimidine rings. We found that the aromaticity exhibits a similar trend for NABs in both their isolated and stacked forms. Moreover, it is indicated that aromatic fluctuation index is more sensible index to delineate the aromaticity variation during stacking interaction.

     

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

文献中已有越来越多的芳香性体系被发现,同时也有越来越多的芳香性指标被提出来,但是如何解释芳香化合物稳定性的起源以及理解芳香性的本质仍然是当今理论化学中一个悬而未决的难题。运用我们新近提出的密度泛函活性理论信息论方法,不久前我们曾对一系列富烯衍生物进行了系统研究并得到了一个全新的认识。本文进一步探讨苯并富烯衍生物的芳香性行为,目的在于考察一个或多个苯环与富烯连接之后其芳香性发生变化的情况。运用香农熵,费舍尔信息,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.     

Molecular structure and bonding character of mono and divalent metal cations (Li+, Na+, K+, Be2+, Mg2+, and Ca2+) with substituted benzene derivatives: AIM,NBO, and NMR analyses

Cation–π complexes between several cations (Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺) and different π-systems such as para-substituted (F, Cl, OH, SH, CH₃, and NH₂) benzene derivatives have been investigated by UB3LYP method using 6-311++G** basis set in the gas phase and the water solution. The ions have shown cation–π interaction with the aromatic motifs. Vibrational frequencies and physical properties such as dipole moment, chemical potential, and chemical hardness of these compounds have been systematically explored. The natural bond orbital analysis and the Bader’s quantum theory of atoms in molecules are also used to elucidate the interaction characteristics of the investigated complexes. The aromaticity is measured using several well-established indices of aromaticity such as NICS, HOMA, PDI, FLU, and FLUπ. The MEP is given the visual representation of the chemically active sites and comparative reactivity of atoms. Furthermore, the effects of interactions on NMR data have been used to more investigation of the studied compounds.     

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.     

Interactions of the 5-fluorouracil anticancer drug with DNA pyrimidine bases: a detailed computational approach

The H-bonded complexes formed from interaction between 5-fluorouracil (FU) and DNA pyrimidine bases have been investigated by B3LYP method using 6-311++G** basis set in the gas phase and the water solution. Vibrational frequencies and physical properties such as dipole moment, chemical potential, and chemical hardness of these compounds have been systematically explored. The natural bond orbital analysis and the Bader’s quantum theory of atoms in molecules are also used to elucidate the interaction characteristics of the investigated complexes. The aromaticity is measured using several well-established indices of aromaticity such as NICS, HOMA, PDI, ATI, and FLU. The MEP is given the visual representation of the chemically active sites and comparative reactivity of atoms. Furthermore, the effects of interactions on NMR data have been used for further investigation of the studied compounds.