Aromaticity of Distorted Benzene Rings: Exploring the Validity of Different Indicators of Aromaticity

Nuclear quantum and H/D isotope effects on aromaticity: path integral molecular dynamics study

Aromaticity is an important concept in organic chemistry, and thus, many theoretical and experimental studies have been conducted so far. However, the majority of theoretical studies have concentrated on the aromaticity of the stationary point structures. Herein, the influence of nuclear quantum fluctuation (nuclear quantum effects: NQEs) and thermal fluctuation on the aromaticity of benzene have been analyzed by path integral molecular dynamics (PIMD) simulation. The PIMD simulations revealed that the NQEs affected not only the C-H bonds but also the C-C bonds. The HOMA and NICS calculations demonstrated that the aromaticity decreased due to the NQEs of carbon atoms, attributed to an increase in the contribution from specific vibrational modes strongly correlated with benzene's aromaticity.

Effect of Benzo-Annelation on Triplet State Energies in Polycyclic Conjugated Hydrocarbons

In a series of earlier studies, the effect of benzo-annelation was found to be a useful tool for tuning the aromaticity in polycyclic conjugated compounds to desired level. In this work we studied the (anti)aromaticity of benzo-annelated derivatives of three conjugated hydrocarbons (anthracene, fluoranthene and biphenylene) in their lowest lying singlet (S0) and triplet (T1) states by means of the energy effect (ef), harmonic oscillator model of aromaticity (HOMA), multicentre delocalization indices (MCI), magnetically induced current densities (MICDs) and nucleus independent chemical shifts (NICS). We showed that benzo-annelation is a topology-based effect which can be used to modify the T1 state excitation energies (E(T1)). A quantitative model was established being able to accurately predict the E(T1) based only on the numbers of angularly, linearly and geminally annelated benzene rings. In addition, it was demonstrated that the E(T1) can be directly related to the (anti)aromatic character of the central ring in the studied molecules in their S0 state.

Tubular Magnetic Shielding Scan (TMSS): A New Technique for Molecular Space Exploration. (i) The Case of Aromaticity of Benzene and [n]Paracyclophanes

Both traditional and novel techniques were employed in this work for magnetic shielding evaluation to shed new light on the magnetic and aromaticity properties of benzene and 12 [n]paracyclophanes with n = 3-14. Density functional theory (DFT) with the B3LYP functional and all-electron Jorge-ATZP and x2c-TZVPPall-s basis sets was utilized for geometry optimization and magnetic shielding calculations, respectively. Additionally, the 6-311+G(d,p) basis set was incorporated for the purpose of comparing the magnetic shielding results. In addition to traditional evaluations such as NICS/NICSzz-Scan, and 2D-3D σiso(r)/σzz(r) maps, two new techniques were implemented: bendable grids (BGs) and cylindrical grids (CGs) of ghost atoms (Bqs). BGs allow for the recording of magnetic shielding from the bent ring levels of [n]pCPs, while CGs provide tubular magnetic shielding scan (TMSS) maps detailing the magnetic shielding from a cylindrical region above and below the ring frame. Our findings suggest that smaller [n]pCPs with n < 6 exhibit deviations in the magnetic shielding above and below the ring, indicating a broken electron delocalization under the ring. In contrast, larger [n]pCPs tend to behave similarly to benzene in terms of magnetic shielding. Moreover, we found that shorter polymethylene chains of [n]pCPs exhibit significantly higher magnetic shielding interactions with the ring. Both of the above techniques offer new and promising tools for characterizing nonplanar aromatic compounds, thereby contributing to a deeper understanding of their magnetic and electronic properties.

Carbon-[n]Triangulenes and Sila-[n]Triangulenes: Which Are Planar?

Using our recently suggested concept of a quasi-molecule ("tile") and, in the case of the planarity here at stake, its generalization to larger than tetratomics, we explain why carbon [n]triangulenes tend to be planar, while hybrids, where just a few or even all a- or b-type carbon atoms are silicon-substituted (sila-[n]triangulenes), tend to be planar/nonplanar when compared with the unsubstituted carbon-[n]triangulenes. Because other spin states of the parent carbon- and sila-[n]triangulenes tend to correlate with the same tiles, it is conjectured that no structural changes are expected to depend on their spin state. Other polycyclic and sila-compounds are also discussed.

Synthesis of precisely functionalizable curved nanographenes via graphitization-induced regioselective chlorination in a mechanochemical Scholl Reaction

While the synthesis of nanographenes has advanced greatly in the past few years, development of their atomically precise functionalization strategies remains rare. The ability to modify the carbon scaffold translates to controlling, adjusting, and adapting molecular properties. Towards this end, here, we show that mechanochemistry is capable of transforming graphitization precursors directly into chlorinated curved nanographenes through a Scholl reaction. The halogenation occurs in a regioselective, high-yielding, and general manner. Density Functional Theory (DFT) calculations suggest that graphitization activates specific edge-positions for chlorination. The chlorine atoms allow for precise chemical modification of the nanographenes through a Suzuki or a nucleophilic aromatic substitution reaction. The edge modification enables modulation of material properties. Among the molecules prepared, corannulene-coronene hybrids and laterally fully π-extended helicenes, heptabenzo[5]superhelicenes, are particularly noteworthy.

Aromaticity controls the excited-state properties of host-guest complexes of nanohoops

π-Conjugated organic molecules have exciting applications as materials for batteries, solar cells, light-emitting diodes, etc. Among these systems, antiaromatic compounds are of particular interest because of their smaller HOMO-LUMO energy gap compared to aromatic compounds. A small HOMO-LUMO gap is expected to facilitate charge transfer in the systems. Here we report the ground and excited-state properties of two model nanohoops that are nitrogen-doped analogs of recently synthesized [4]cyclodibenzopentalenes - tetramers of benzene-fused aromatic 1,4-dihydropyrrolo[3,2-b]pyrrole ([4]DHPP) and antiaromatic pyrrolo[3,2-b]pyrrole ([4]PP). Their complexes with C60 fullerene show different behavior upon photoexcitation, depending on the degree of aromaticity. [4]DHPP acts as an electron donor, whereas [4]PP is a stronger electron acceptor than C60. The ultrafast charge separation combined with the slow charge recombination that we found for [4]PP⊃C60 indicates a long lifetime of the charge transfer state.

A benchmark study of aromaticity indexes for benzene, pyridine and the diazines - I. Ground state aromaticity

Five different aromaticity indexes are benchmarked for benzene, pyridine and the diazines in their ground states. A basis set study was performed using the Pople style, Karlsruhe and Dunning's correlation consistent basis sets. Ten different DFT functionals, including LSDA, PBE, PBE0, B3LYP, CAM-B3LYP, wB97XD, M06-2X, SOGGA11X, M11 and MN15 were benchmarked by comparison with CCSD, CASSCF and MP2. Large out-of-plane imaginary frequencies were observed for some of the optimized structures at the correlated wavefunction level of theory. It was found that the DFT functionals in general predict the para-delocalization index, multicenter index and aromatic fluctuation index to be approximately 70%, 50% and 45% larger, respectively, compared to the CCSD method. Comparisons of the DFT functionals showed that the wB97XD, CAM-B3LYP and M06-2X functionals performed the best. Furthermore, the basis set dependence of the DFT functionals was found to be large for the electron sharing indexes. Based on these findings, it is recommended to perform ground state calculations of aromaticity indexes at the wB97XD, CAM-B3LYP or M06-2X level of theory utilizing a simple basis set of triple-ζ quality.

Five different aromaticity indexes are benchmarked for benzene, pyridine and the diazines in their ground states.

Valence non-Lewis density as an approach to describe and measure aromaticity of organic and inorganic molecules

Based on the linear combination of atomic orbital-molecular orbital by the natural bond orbitals (NBO) theory, the attractive donor-acceptor superposition interaction between filled (Lewis-type) and vacant (non-Lewis-type) orbitals provide a general mechanism for quantal energy lowering. This interaction has a direct impact on the quantity of the second-order stabilization energy. Therefore, the valence non-Lewis density (VNLD) index, the electron density of unoccupied valence nonbonding and antibonding orbitals, is introduced as an approach to describe and measure aromaticity. This index is based on the frontier orbital concept. To investigate the validity of the proposed aromaticity index, we selected several test sets of organic and inorganic molecules such as different ring sizes in cyclic and heterocyclic hydrocarbons, and all-metal and semimetal clusters, and compared our findings with previous aromaticity analysis. According to the results, VNLD values are well correlated and anticipated the order of aromaticity with the formerly introduced criteria. Furthermore, VNLD reveals that the rings with π-sextet electrons localized in a ring are more aromatic than the other rings, thus, it is in line with Clar's rule. Our proposed aromaticity index has advantages such as, easy to obtain from NBO analysis, and does not require reference molecules which made it more applicable for realizing the aromaticity order in many organic and inorganic compounds.

A Highly Strained All-Phenylene Conjoined Bismacrocycle

Herein, we report the precise synthesis of a 3D highly strained all-phenylene bismacrocycle, termed conjoined (1,4)[10]cycloparaphenylenophane (SCPP[10]). This structure consists of a twisted benzene ring which is bridged twice by phenylene units anchored in two para-positions. The conjoined structure of SCPP[10] was confirmed in real space at the atomic scale by scanning tunneling microscopy. Theoretical calculations indicate that this bismacrocycle has a very high strain energy of 110.59 kcal mol^-1 and the largest interphenylene torsion angle of 46.07° caused by multiple repulsive interactions. Furthermore, a 1:2 host-guest complex of SCPP[10] and [6,6]-phenyl-C₆₁ -butyric acid methyl ester was investigated, which represents the first peanut-shaped 1:2 host-guest complex based on bismacrocycles.

Quantum chemical determination of molecular geometries and spectral investigation of 4-ethoxy-2, 3-difluoro benzamide

The present work reports the application of density functional theory (DFT) at B3LYP with various basis sets which provide the relationship between the structural and spectral properties of 4-ethoxy-2, 3-difluoro benzamide (4EDFB). A Complete vibrational analysis has been performed at the density functional theory (DFT) method with various basis sets in the ground state. The results of vibrational wave numbers are in good agreement with the experimental spectra (Infrared and Raman). Energy gap of the molecule is evaluated using frontier molecular orbital energies (HOMO-LUMO). The frontier energy gap value reveals the chemical reactivity and intermolecular charge transfer occur within the molecule. Global chemical descriptors provide the local and global softness and local reactivity parameters used to identify the nucleophilic and electrophilic behavior of a specific site within the compound. The dimer structure is performed to evaluate the intermolecular hydrogen bond (O-H-O). The title molecule is capable of receiving second harmonic generation (SHG) is due to high value of hyperpolarizability indicates the NLO activity of the molecule. Apart from NLO entities, aromaticity and the molecular electrostatic potential surface (MEP) explain the hydrogen bonding and provide the reactive behavior of the molecule. The Mulliken population analysis leads to redistribution of electron density in the ring.

Aromaticity and Electron Density of Hypericin

The influence of the substituents on the geometry of the central ring system of hypericin has been analyzed. Substitution that causes flattening of the hypericin central rings is connected with introducing the aromatic character of the empty rings. All the hypericin rings have an aromatic character illustrated by the Harmonic Oscillator Measure of Aromaticity (HOMA), Nucleus Independent Chemical Shift (NICS), Fluctuation Index (FLU), and Ellipticity Index (EL) indices. Quantum Theory of Atoms in Molecules (QTAIM) and Natural Bond Orbital (NBO) analyses performed on 7,14-dihydrophenanthro[1,10,9,8-opqra]perylene, its substituted analogues, and hypericin show an influence of this substitution on electron density of the central rings.

Dual Activation of Aromatic Diels-Alder Reactions

The unusually fast Diels–Alder reactions of [5]cyclophanes were analyzed by DFT at the BLYP‐D3(BJ)/TZ2P level of theory. The computations were guided by an integrated activation‐strain and Kohn–Sham molecular orbital analysis. It is revealed why both [5]metacyclophane and [5]paracyclophane exhibit a significant rate enhancement compared to their planar benzene analogue. The activation strain analyses revealed that the enhanced reactivity originates from 1) predistortion of the aromatic core resulting in a reduced activation strain of the aromatic diene, and/or 2) enhanced interaction with the dienophile through a distortion‐controlled lowering of the HOMO–LUMO gap within the diene. Both of these physical mechanisms and thus the rate of Diels–Alder cycloaddition can be tuned through different modes of geometrical distortion (meta versus para bridging) and by heteroatom substitution in the aromatic ring. Judicious choice of the bridge and heteroatom in the aromatic core enables effective tuning of the aromatic Diels–Alder reactivity to achieve activation barriers as low as 2 kcal mol⁻¹, which is an impressive 35 kcal mol⁻¹ lower than that of benzene.

Molecular Photoswitching Aided by Excited-State Aromaticity

Central to the development of optoelectronic devices is the availability of efficient synthetic molecular photoswitches, the design of which is an arena where the evolving concept of excited-state aromaticity (ESA) is yet to make a big impact. The aim of this minireview is to illustrate the potential of this concept to become a key tool for the future design of photoswitches. The paper starts with a discussion of challenges facing the use of photoswitches for applications and continues with an account of how the ESA concept has progressed since its inception. Then, following some brief remarks on computational modeling of photoswitches and ESA, the paper describes two different approaches to improve the quantum yields and response times of switches driven by E/Z photoisomerization or photoinduced H-atom/proton transfer reactions through simple ESA considerations. It is our hope that these approaches, verified by quantum chemical calculations and molecular dynamics simulations, will help stimulate the application of the ESA concept as a general tool for designing more efficient photoswitches and other functional molecules used in optoelectronic devices.

Aromaticity of benzene derivatives: an exploration of the Cambridge Structural Database

The harmonic oscillator model of aromaticity (HOMA) index, one of the most popular aromaticity indices for solid-state benzene rings in the Cambridge Structural Database (CSD), has been analyzed. The histograms of HOMA for benzene, for benzene derivatives with one formyl, nitro, amino or hydroxy group as well as the histograms for the derivatives with two formyl, nitro, amino or hydroxy groups inortho,metaandparapositions were investigated. The majority of the substituted benzene derivatives in the CSD are characterized by a high value of HOMA, indicating fully aromatic character; however, the distribution of the HOMA value from 1 to about 0 indicates decreasing aromaticity down to non-aromatic character. Among the benzene derivatives investigated, a significant decrease in aromaticity can be related to compounds with diamino and dinitro groups in themetaposition.

The electron density of delocalized bonds (EDDB) applied for quantifying aromaticity

In this study the recently developed electron density of delocalized bonds (EDDB) is used to define a new measure of aromaticity in molecular rings. The relationships between bond-length alternation, electron delocalization and diatropicity of the induced ring current are investigated for a test set of representative molecular rings by means of correlation and principal component analyses involving the most popular aromaticity descriptors based on structural, electronic, and magnetic criteria. Additionally, a qualitative comparison is made between EDDB and the magnetically induced ring-current density maps from the ipsocentric approach for a series of linear acenes. Special emphasis is given to the comparative study of the description of cyclic delocalization of electrons in a wide range of organic aromatics in terms of the kekulean multicenter index KMCI and the newly proposed EDDB^k index.

New electron delocalization tools to describe the aromaticity in porphyrinoids

There are several possible pathways in the macrocycle of large porphyrinoids and, among aromaticity indices, only AV_minis capable of recognizing the most aromatic one.

Magnetic resonance energy and topological resonance energy

Ring-current diamagnetism of a polycyclic π-system is closely associated with thermodynamic stability due to the individual circuits. Magnetic resonance energy (MRE), derived from the ring-current diamagnetic susceptibility, was explored in conjunction with graph-theoretically defined topological resonance energy (TRE). For many aromatic molecules, MRE is highly correlative with TRE with a correlation coefficient of 0.996. For all π-systems studied, MRE has the same sign as TRE. The only trouble with MRE may be that some antiaromatic and non-alternant species exhibit unusually large MRE-to-TRE ratios. This kind of difficulty can in principle be overcome by prior geometry-optimisation or by changing spin multiplicity. Apart from the semi-empirical resonance-theory resonance energy, MRE is considered as the first aromatic stabilisation energy (ASE) defined without referring to any hypothetical polyene reference.

The role of the long-range exchange corrections in the description of electron delocalization in aromatic species

In this article, we address the role of the long-range exchange corrections in description of the cyclic delocalization of electrons in aromatic systems at the density functional theory level. A test set of diversified monocyclic and polycyclic aromatics is used in benchmark calculations involving various exchange-correlation functionals. A special emphasis is given to the problem of local aromaticity in acenes, which has been a subject of long-standing debate in the literature. The presented results indicate that the noncorrected exchange-correlation functionals significantly overestimate cyclic delocalization of electrons in heteroaromatics and aromatic systems with fused rings, which in the case of acenes leads to conflicting local aromaticity predictions from different criteria. © 2017 Wiley Periodicals, Inc.

The nature of the Au–N bond in gold(iii) complexes with aromatic nitrogen-containing heterocycles: the influence of Au(iii) ions on the ligand aromaticity

Aromatic character of nitrogen-containing rings is decreased upon gold(iii) coordination.

Aromaticity and Through-Space Interaction between Aromatic Rings in [2.2]Paracyclophanes

The HOMA index calculated for [2.2]paracyclophanes in the solid state reveals a slight decrease of aromaticity. Interactions between aromatic rings of [2.2]paracyclophane have been investigated using AIM and NCI analysis in both crystal and optimized [2.2]paracyclophane structures. AIM analysis showed that the C···C bond path between the two aromatic rings is present only in few [2.2]paracyclophanes. The NCI approach visualized the dispersion and repulsive interactions between the aromatic rings of every [2.2]paracyclophane. Combination of AIM and the NCI approach is necessary for determining and identifying nonbonded interactions in [2.2]paracyclophanes.

Perimeter ring currents in benzenoids from Pauling bond orders

It is shown that the ring currents in perimeter hexagonal rings of Kekulean benzenoids, as estimated within the Randić conjugated-circuit model, can be calculated directly without tedious pairwise comparison of Kekulé structures or Kekulé counting for cycle-deleted subgraphs. Required are only the Pauling bond orders of perimeter bonds and the number of Kekulé structures of the benzenoid, both readily available from the adjacency matrix of the carbon skeleton. This approach provides easy calculation of complete current maps for benzenoids in which every face has at least one bond on the perimeter (as in the example of cata-condensed benzenoids), and allows qualitative evaluation of the main ring-current contributions to (1)H chemical shifts in general benzenoids. A combined Randić-Pauling model for correlation of ring current and bond length through bond order is derived and shown to be consistent with resilience of current under bond alternation.

The effect of Li doping on the nonlinear optical properties of [2.2]paracyclophane

The similar molecules [2.2]paracyclophane (22PCP) and 1,1,2,2,9,9,10,10-octafluoro[2.2]paracyclophane (8F22PCP) have both generated considerable synthetic interest since they were first prepared. In this work, the nonlinear optical properties of 22PCP, 8F22PCP, and the related Li-doped systems 22PCP-Li and 8F22PCP-Li (which have a Li atom above 22PCP and 8F22PCP, respectively) were investigated. An analysis of natural bond orbital charges showed that there is greater charge transfer from the Li atom to the benzene rings in 8F22PCP-Li than in 22PCP-Li. The variation in the calculated nucleus independent chemical shift (NICS) value as a function of the distance from the lower benzene ring towards the upper benzene ring was found to be W-shaped for both 22PCP and 22PCP-Li. Moreover, whereas all of the NICS values of 22PCP and 22PCP-Li were markedly negative, all of the NICS values of 8F22PCP and 8F22PCP-Li were either positive or only moderately negative. Calculations of the electro-optical properties of these systems showed that the first hyperpolarizability of 22PCP-Li was noticeably larger than that of 8F22PCP-Li. According to the two-level model, the larger first hyperpolarizability of 22PCP-Li is due to its smaller transition energy.

Scalar and Spin-Orbit Relativistic Corrections to the NICS and the Induced Magnetic Field: The case of the E12(2-) Spherenes (E = Ge, Sn, Pb)

Can relativistic effects modify the NICS and the B(ind) values? In this manuscript we evaluate the relativistic corrections incorporated via the zeroth-order regular approximation to the calculations of nucleus-independent chemical shifts and the induced magnetic field (B(ind)) in the E12(2-) spherenes (E = Ge, Sn, Pb). We found that both electron delocalization descriptors are strongly affected by the relativistic corrections. For instance, for plumbaspherene, the difference in values from the nonrelativistic to the relativity-included calculation is almost 40 ppm! Our results show that the changes observed in the NICS and B(ind) values in the title cages are a consequence of the treatment of the relativistic effects. If these effects are included as scalar or spin-orbit calculations, then we can establish the next trend: Ge12(2-) is a nonaromatic species, Sn12(2-) is a low aromatic species, and Pb12(2-) is strongly aromatic, according to calculated NICS and B(ind) values. Thus, any prediction of electron delocalization in molecules containing heavy elements without considering an adequate treatment for relativistic effects may lead to an erroneous chemical interpretation.

Magnetic criteria of aromaticity

This review describes the current state of magnetic criteria of aromaticity. The introduction contains the fundamentals of ring currents in aromatic and antiaromatic systems, followed by a brief description of experimental and computational tools: NMR, diamagnetic susceptibility exaltation, current density analyses (CDA) and nucleus independent chemical shifts (NICS). This is followed by more comprehensive chapters: NMR - focusing on the work of R. Mitchell - NICS and CDA - describing the progress and development of the methods to their current state and presenting some examples of representative work.

The non-covalent nature of the molecular structure of the benzene molecule

The benzene molecule is one of the most emblematic systems in chemistry, with its structural features being present in numerous different compounds. We have carried out an analysis of the influence of quantum mechanical interference on the geometric features of the benzene molecule, showing that many of the characteristics of its equilibrium geometry are a consequence of non-covalent contributions to the energy. This result implies that quasi-classical reasoning should be sufficient to predict the defining aspects of the benzene structure such as its planarity and equivalence of its bond lengths.

Quantifying aromaticity with electron delocalisation measures

Aromaticity cannot be measured directly by any physical or chemical experiment because it is not a well-defined magnitude. Its quantification is done indirectly from the measure of different properties that are usually found in aromatic compounds such as bond length equalisation, energetic stabilisation, and particular magnetic behaviour associated with induced ring currents. These properties have been used to set up the myriad of structural-, energetic-, and magnetic-based indices of aromaticity known to date. The cyclic delocalisation of mobile electrons in two or three dimensions is probably one of the key aspects that characterise aromatic compounds. However, it has not been until the last decade that electron delocalisation measures have been widely employed to quantify aromaticity. Some of these new indicators of aromaticity such as the PDI, FLU, ING, and INB were defined in our group. In this paper, we review the different existing descriptors of aromaticity that are based on electron delocalisation properties, we compare their performance with indices based on other properties, and we summarise a number of applications of electronic-based indices for the analysis of aromaticity in interesting chemical problems.

Revisiting Aromaticity and Chemical Bonding of Fluorinated Benzene Derivatives

The electron delocalization of benzene (C6H6) and hexafluorobenzene (C6F6) was analyzed in terms of the induced magnetic field, nucleus-independent chemical shift (NICS), and ring current strength (RCS). The computed out-of-plane component of the induced magnetic field at a distance (r) greater than or equal to 1.0 Å above the ring center correlates well (R (2)>0.99) with the RCS value. According to these criteria, fluorination has two effects on the C6 skeleton; concomitantly, the resonant effects diminish the π electron delocalization and the inductive effects decrease the charge density at the ring center and therefore reduce the magnitude of the paratropic current generated in this region. The equilibrium between both effects decreases aromaticity in the fluorinated benzene derivatives. These results can be extrapolated to determine the aromaticity of any derivative within the series of fluorinated benzene derivatives (C6H(6-n)F n , where n=1-5).