Partitioning of π-Electrons in Rings for Clar Structures of Benzenoid Hydrocarbons

Quantitative Resonance Theory Based on the Clar Sextet Model

Despite its great explanatory power in understanding the chemistry of polycyclic aromatic hydrocarbons (PAHs) and related systems, the Clar sextet rule still remains an intuitive and qualitative model with notable exceptions in some cases. Here we develop a quantitative theory of chemical resonance based on semilocalized Clar-type resonance structures (named the Clar resonators) consisting of variable numbers of Clar sextets and C═C bonds. The constructed wave functions of the Clar resonators are used to expand the actual wave function of the π-conjugated system obtained from a DFT or Hartree-Fock calculation. The resultant weights and one-electron energies of the Clar resonators can serve as a quantitative measure of the importance of these resonators. Implementing the theory in our open-source python code EzReson and applying it to over a thousand PAH molecules of different sizes and shapes, we show that the weight of the Clar resonators increases exponentially with increasing number of sextets and that their energy decreases linearly with the latter, thus confirming the general validity of the Clar rule. On the basis of such a large-scale resonance analysis, we propose three extended Clar rules, along with a unified quantitative model, that are able to evaluate the importance of all Clar resonators and the ring aromaticity for PAHs. Using the present theories, we uncover the essential role that the minor Clar resonators may play in correctly understanding the resonance stabilization and local aromaticity of rings, which was totally overlooked in the original Clar model.

A critical approach toward resonance-assistance in the intramolecular hydrogen bond interaction of 3,5-diiodosalicylic acid: a spectroscopic and computational investigation

The photophysics of a prospective drug molecule, 3,5-diiodosalicylic acid (3,5-DISA), having a wide spectrum of biological and medicinal applications, have been investigated using spectroscopic techniques and computational analyses. The remarkably large Stokes' shifts in various solvents from 3,5-DISA has been intertwined with the occurrence of an excited-state intramolecular proton transfer (ESIPT) reaction. Concurrently, the emergence of an intriguing dual emission feature in less interacting solvents is also reported and the spectral response of 3,5-DISA toward the variation of medium acidity/basicity has been exploited to decipher the nature of various species present in different solvents. Our experimental results, unveiling the occurrence of an ESIPT reaction in 3,5-DISA, have been aptly substantiated from computational studies in which the operation of ESIPT has been explored from structural as well as energetics (analysis of potential energy surface (PES)) perspectives. A major focus of the present study is on the evaluation of the intramolecular hydrogen bond (IMHB) interaction in 3,5-DISA, including the application of various methodologies to estimate the IMHB energy and subsequently, an in-depth analysis of the IMHB interaction reveals its partially covalent nature through the application of advanced quantum chemical tools, e.g., the natural bond orbital (NBO) method. In this context, the interplay between the aromaticity of the benzene nucleus and the IMHB energy has been rigorously explored, showing indications for the occurrence of resonance-assisted hydrogen bonding (RAHB) in 3,5-DISA. To this end, the geometric as well as magnetic criteria of aromaticity have been analyzed.

Quantum chemical exploration of the intramolecular hydrogen bond interaction in 2-thiazol-2-yl-phenol and 2-benzothiazol-2-yl-phenol in the context of excited-state intramolecular proton transfer: a focus on the covalency in hydrogen bond

The present work demonstrates a computational exploration of the intramolecular H-bond (IMHB) interaction in two model heterocyclic compounds - 2-thiazol-2-yl-phenol (2T2YP) and 2-benzothiazol-2-yl-phenol (2B2YP) by meticulous application of various quantum chemical tools. Major emphasis is rendered on the analysis of IMHB interaction by calculation of electron density ρ(r) and Laplacian ∇(2)ρ(r) at the bond critical point using the Atoms-In-Molecule methodology. Topological features based on ρ(r) suggest that at equilibrium geometry the IMHB interaction develops certain characteristics typical of a covalent interaction. The interplay between aromaticity and Resonance-Assisted H-Bond (RAHB) has also been discussed using both geometrical and magnetic criteria. The occurrence of IMHB interaction in 2T2YP and 2B2YP has also been criticized under the provision of the Natural Bond Orbital (NBO) analysis. The ESIPT phenomenon in the molecular systems is also critically addressed on the lexicon of potential energy surface (PES) analysis.

Conjugated circuits currents in hexabenzocoronene and its derivatives formed by joining proximal carbons

We report on calculated CC bond currents for a dozen derivatives of hexabenzocoroenene in which one or more proximal carbon atoms at the molecular periphery have been bridged. The approach that we use is graph-theoretical in nature, following our outline of this method in 2003, which is based on finding all conjugated circuits in all Kekulé valence structures of these molecules. To the π-electrons having 4n + 2 π-electrons are assigned anticlockwise π-electron currents and to conjugated circuits having 4n π-electrons are assigned π-electron currents. One may summarize the results reported in this work by stating that CC bond currents in the compounds considered decrease on going from peripheral rings to the central ring of the molecule, and also that CC bond currents decrease by insertion of bridges to proximal peripheral benzenoid rings.

Using Clar sextets for two- and three-dimensional aromatic systems

After a brief history of the aromaticity concept, the use of Clar sextet circles is reviewed for explaining various aspects of planar aromatic systems (benzenoids, heterocycles) and of tridimensional carbon aggregates. When folding a graphene sheet for obtaining nanotubes, nanotori, or nanocones, the congruence of Clar sextet circles allows the classification of all such aggregates into two classes (congruent or incongruent) with marked differences in properties; this is in agreement with the well-known condition h - k ≡ 0 (mod 3), equivalent to congruence, in terms of the chiral vectors h and k for graphene sheets.

A new approach for aromaticity criterion based on electrostatic field gradient

In this research, electric field gradient (EFG), which is the first derivative of electric field, is applied for evaluation of aromaticity of 89 cyclic organic compounds. In our calculations, DFT procedure (b3lyp) with basis set 6-311++G has been employed, and the obtained electronic structures and frequency test has been done for optimized geometries. The aromaticity evaluated for these compounds by EFG procedure is successfully compared with other well-known indices in literature, especially with nuclear independent chemical shift (NICS). These comparisons show that EFG method of assessment of aromaticity can be used as a rather valid procedure for this purpose. Flexibility and simplicity of EFG make this method a rather easy procedure for assessment of aromaticity. Since EFG method of aromaticity evaluation does not need specific programming and it can be done by known software such as Gaussian, therefore, the availability for everyone to calculate desired aromaticity by this method is one of attractive feature of it similar to NICS.

Partitioning of π-electrons in rings of aza-derivatives of naphthalene

A recently proposed method for calculating the ?-electron contents (EC) of rings of heteroatom-containing polycyclic conjugated molecules was applied to the aza-derivatives of naphthalene. The main finding was that a nitrogen atom in position ? (resp.?) diminishes (resp. increases) the EC-value of the respective ring. Such a regularity in the displacement of ?-electrons can be (qualitatively) rationalized by means of resonance-theoretical reasoning. .

Chemical graph theory andn-center electron delocalization indices: A study on polycyclic aromatic hydrocarbons

Relations between aromaticity indices derived from chemical graph theory and those based on 6-center electron delocalization are investigated for a series of polybenzenoid hydrocarbons. Aromatic stabilization obtained by means of the effective scaled electron delocalization is highly correlated to the resonance energy, RE, obtained both from SCF MO calculations and conjugated ring circuits model. Local aromaticity of benzene rings is discussed using two different criteria, in one of them aromaticity is just given by the cyclic pi-electron conjugation of the ring, whereas terms involving more than one ring are also considered in the other one. Indices derived from chemical graph theory and those obtained from the 6-center electron delocalization give rise to the same local aromaticity. Moreover, 6-center electron delocalization provides more quantitative information.

Topology-Driven Physicochemical Properties of π-Εlectron Systems. 1. Does the Clar Rule Work in Cyclic π-Electron Systems with the Intramolecular Hydrogen or Lithium Bond?

The Clar model predicting stability and electron distribution in benzenoid hydrocarbons seems to be also a good predictor for related properties in lithium o-acylphenolates and to a lesser extent in the phenols themselves. This conclusion is based on analysis of geometry changes in the analogues of naphthalene, phenanthrene, anthracene, and triphenylene, where benzene rings are systematically replaced with a quasi-ring formed as a beta-ketoenol or beta-ketoenolate complex with a lithium cation, i.e., where CH-CH-CH fragments are replaced with O...Li...O or O...H...O fragments. These systems were optimized at the MP2/6-31G(2d,p) level of theory. The energy of bond separation reactions in line with aromaticity indices HOMA and NICSs supported the above statement.