Characterization of resonance structures in aromatic rings of benzene and its heavier-element analogues

We present the experimental visualization of the valence-electron-density distribution in benzene and its kinetically stabilized heavier-element analogues, i.e., 1,2-disilabenzene and 1,2-digermabenzene. The valence-electron-density-distribution (EDD) analysis on the 1,2-disila- and 1,2-digermabenzenes revealed that these contain incompletely delocalized π electrons on their cyclic conjugation systems, making them less aromatic compared to benzene. Based on the results of this EDD analysis in combination with anisotropy of the current-induced density (ACID) calculations, considerable contributions from the characteristic resonance structures of 1,2-disila- and 1,2-digermabenzenes with cleaved EE bonds can be expected.

Evaluation of Slight Changes in Aromaticity through Electronic and Density Functional Reactivity Theory-Based Descriptors

Aromaticity is a useful tool to rationalize the structure, stability, and reactivity in several compounds. Although aromaticity is not directly an observable, it is well accepted that electronic delocalization around the molecular ring is a key stabilizing feature of aromatic compounds. This contribution presents a systematic evaluation of the capability of delocalization and reactivity criteria to describe aromaticity in a set of fluorinated benzenes. The aromaticity indices are compared with quantities obtained from the magnetic criteria of aromaticity, i.e., the strength of the ring current induced by an external magnetic field and the popular NICS _zz (1) index. In this evaluation, the indices based on delocalization criteria used are aromatic fluctuation index (FLU), para-delocalization index (PDI), PDI_π, and the multicenter delocalization index (MCI). In addition, indices based on the bifurcation values of scalar functions are derived from electron density such as electron localization function (the π contribution, ELF_π) and the π contribution of the localized orbital locator (LOL_π). Furthermore, reactivity indices based on chemical reactivity and the information-theoretic (reactivity) approach are para-linear response (PLR), Shannon entropy, Fisher information, and Ghosh-Berkowitz-Parr (GBP) entropy. The results obtained show that the delocalization-based indicators present a high sensitivity to slight changes in aromaticity and that the reactivity criterion can be considered as a complementary tool for the study of this phenomenon, even when these changes are minimal. These results encourage the use of multiple indicators for a complete understanding of aromaticity in various chemical compounds.

Energetic and Geometric Characteristics of Substituents, Part 3: The Case of NO2 and NH2 Groups in Their Mono-Substituted Derivatives of Six-Membered Heterocycles

Substituted heterocyclic arenes play important roles in biochemistry, catalysis, and in the design of functional materials. Exemplary six-membered heteroaromatic molecules, that differ from benzene by inclusion of one heteroatom, are pyridine, phosphorine, arsabenzene, and borabenzene. This theoretical study concerns the influence of the heteroatom present in these molecules on the properties of substituents of two types: electron-donating (ED) NH2 group and electron-accepting (EA) NO2 group, attached at the 2-, 3-, or 4-position. The effect is evaluated by the energy of interaction (Erel) between the substituent and the substituted system and electronic properties of the substituents described by the charge of the substituent active region (cSAR) index. In addition, several geometric descriptors of the substituent and heteroaromatic ring, as well as changes in the aromaticity, are considered. The latter are assessed using the Electron Density of Delocalized Bonds (EDDBs) property of delocalized π electrons. The obtained results show that the electronegativity (EN) of the heteroatom has a profound effect on the EA/ED properties of the substituents. This effect is also reflected in the geometry of studied molecules. The Erel parameter indicates that the relative stability of the molecules is highly related to the electronic interactions between the substituent and the heteroarene. This especially applies to the enhancement or weakening of π-resonance due to the EN of the heteroatom. Additionally, in the 2-heteroarene derivatives, specific through-space ortho interactions contribute to the heteroatom effects.

Revival of Hückel Aromatic (Poly)benzenoid Subunits in Triplet State Polycyclic Aromatic Hydrocarbons by Silicon Substitution

By employing density functional theory (DFT) calculations we show that mono- and disilicon substitution in polycyclic aromatic hydrocarbons, having two to four benzene units, quenches their triplet state antiaromaticity by creating Hückel aromatic (poly)benzenoid subunit(s) and weakly antiaromatic, or almost nonaromatic silacycle. Therefore, such systems are predicted to be globally aromatic in both the ground state and the first excited triplet state. Putting the silicon atom(s) into various positions of a hydrocarbon provides an opportunity to tune the singlet-triplet energy gaps. They depend on the global aromaticity degree which, in turn, depends on the type of aromatic carbocyclic subunit(s) and the extent of their aromaticity. On the basis of the set of studied compounds, some preliminary rules on how to regulate the extent of global, semiglobal and local aromaticity are proposed. The results of this work extend the importance of Hückel aromaticity concept to excited triplet states which are usually characterized by the Baird type of (anti)aromaticity.

Triplet state (anti)aromaticity of some monoheterocyclic analogues of benzene, naphthalene and anthracene

By employing DFT calculations, we show the influence of heteroatom substitution on the triplet state (anti)aromaticity of benzene, naphthalene and anthracene.

Heteroatom effects on aromaticity of five-membered rings in acenaphthylene analogs

The pattern of cyclic conjugation was thoroughly studied in the series of N- and P-acenaphthylene derivatives using several different aromaticity indices: the energy effect (ef), multicenter delocalization index (MCI), harmonic oscillator model of aromaticity (HOMA) index, and nucleus independent chemical shifts (NICS). The Kekulé-structure-based reasoning predicts that there would be no cyclic conjugation in the "empty" five-membered heteroatom-containing rings in the studied molecules. It was found that according to the ef, MCI, and HOMA values, the extent of cyclic conjugation in the pentagonal rings is strongly influenced by the number and mutual arrangement of the hexagonal rings. In addition, it was revealed that in some of the examined molecules, the intensity of cyclic conjugation in the "empty" pentagons is even stronger than that of some hexagonal rings within the same molecule. The obtained results refute what one would expect based on "chemical intuition," which is usually strongly rooted to the Kekulé structures.Graphical abstract.

Aromaticity of unsaturated BEC4 heterocycles (E = N, P, As, Sb, O, S, Se, Te)

A compendium of pnictogen and chalcogen substituted boron heterocycles were assessed for their aromatic character by first principles density functional theory. Group-15 and Group-16 elements were placed at the ortho-, meta-, and para-positions of six-membered rings relative to boron to assess their impact on the aromaticity of the unsaturated heterocycles. Aromaticity was analyzed by a multidimensional approach using nuclear independent chemical shifts, gauge-including magnetically induced current, as well as natural bond orbital and natural resonance theory analyses. Based on these methods, we observe a general decline of aromaticity in heavier pnictaborines while the chalcogen analogues maintain relatively strong aromatic character. These general trends result from complementary π-π* natural bond order interactions that sustain resonance within the ring of each heterocycle establishing a pattern of cyclic delocalization. Consequently, natural resonance theory displays strong resonance, which is corroborated with the signed modulus of ring current, toroidal vortices of current maps, and elevated average induced current throughout the ring. The 1,3-configurations for pnictaborines and chalcogenaborines are generally more aromatic compared to the 1,2- and 1,4-isomers, which contain π-holes that limit diatropism within the heterocycles. However, an energetic trend favors the 1,2-heterocycles in both groups, with a few exceptions driven in large-part by π-donation of the lone pair on the heteroatom to the pz orbital on the adjacent boron resulting in stabilization. The importance of planarity for high aromaticity is demonstrated, especially in the pnictaborine isomers where pyramidalization at the pnictogen is favored, while bond regularity seems a less important criterion.

Interpreting Aromaticity and Antiaromaticity through Bifurcation Analysis of the Induced Magnetic Field

In all molecules, a current density is induced when the molecule is subjected to an external magnetic field. In turn, this current density creates a particular magnetic field. In this work, the bifurcation value of the induced magnetic field is analyzed in a representative set of aromatic, non-aromatic and antiaromatic monocycles, as well as a set of polycyclic hydrocarbons. The results show that the bifurcation value of the ring-shaped domain adequately classifies the studied molecules according to their aromatic character. For aromatic and nonaromatic molecules, it is possible to analyze two ring-shaped domains, one diatropic (inside the molecular ring) and one paratropic (outside the molecular ring). Meanwhile, for antiaromatic rings, only a diatropic ring-shaped domain (outside the molecular ring) is possible to analyze, since the paratropic domain (inside the molecular ring) is irreducible with the maximum value (attractor) at the center of the molecular ring. In some of the studied cases, i. e., in heteroatomic species, bifurcation values do not follow aromaticity trends and present some inconsistencies in comparison to ring currents strengths, showing that this approximation provides only a qualitative estimation about (anti)aromaticity.

The Paternò–Büchi reaction – a comprehensive review

This review represents the most complete description of the scientific results obtained on a photochemical reaction described 110 years ago by an Italian scientist.

Impact of heteroatoms (S, Se, and Te) on the aromaticity of heterocirculenes

Computations reveal the structural and energetic aspects of aromaticity in heterocirculenes.

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.

Aromatic stabilization of functionalized corannulene cations

The first comprehensive theoretical investigation of aromaticity in functionalized corannulene cations of general formula [CH3-C20H10](+) was accomplished. The experimentally known system [CH3-hub-C20H10](+) was augmented by two other possible isomers, namely, rim- and spoke-ones. Changes in aromaticity, when going from neutral corannulene to its functionalized cations, were monitored with the help of descriptors of different nature such as structure-based HOMA, topological PDI and FLU, and magnetic NICS. A highly efficient tool for analysis and visualization of delocalization and conjugation named ACID was also utilized. In the final step, a complete set of (1)H and (13)C chemical shifts was calculated and compared with the available experimental data. Conservation of aromaticity of 6-membered rings along with vanishing anti-aromatic character of central 5-membered rings was found to be the main reason for the exceptional stability of the hub-isomer. At the same time, functionalization of the corannulene moiety at the rim- or spoke-site resulted in dramatic elimination of aromaticity of 6-membered rings, whereas anti-aromatic character of the central ring remained. Altogether, it led to much lower stability of these isomers in comparison with that of the hub-one.

Lewis acid–base 1,2-addition reactions: synthesis of pyrylium borates from en-ynoate precursors

Treatment of methyl (Z)-2-alken-4-ynoates with the strong Lewis acid tris(pentafluorophenyl) borane, B(C₆F₅)₃, yield 2,5,6-substituted zwitterionic pyrylium borate species via an intramolecular 6-endo-dig cyclisation reaction.

New Approximation to the Third-Order Density. Application to the Calculation of Correlated Multicenter Indices

In this work we present the formulas for the calculation of exact three-center electron sharing indices (3c-ESI) and introduce two new approximate expressions for correlated wave functions. The 3c-ESI uses the third-order density, the diagonal of the third-order reduced density matrix, but the approximations suggested in this work only involve natural orbitals and occupancies. In addition, the first calculations of 3c-ESI using Valdemoro's, Nakatsuji's and Mazziotti's approximation for the third-order reduced density matrix are also presented for comparison. Our results on a test set of molecules, including 32 3c-ESI values, prove that the new approximation based on the cubic root of natural occupancies performs the best, yielding absolute errors below 0.07 and an average absolute error of 0.015. Furthemore, this approximation seems to be rather insensitive to the amount of electron correlation present in the system. This newly developed methodology provides a computational inexpensive method to calculate 3c-ESI from correlated wave functions and opens new avenues to approximate high-order reduced density matrices in other contexts, such as the contracted Schrödinger equation and the anti-Hermitian contracted Schrödinger equation.

Study of Electron Delocalization in 1,2-, 1,3-, and 1,4-Azaborines Based on the Canonical Molecular Orbital Contributions to the Induced Magnetic Field and Polyelectron Population Analysis

The electron delocalization in 1,2-azaborine, 1,3-azaborine, and 1,4-azaborine is studied using canonical molecular orbital contributions to the induced magnetic field (CMO-IMF) method and polyelectron population analysis (PEPA). Contour maps of the out-of-plane component of the induced magnetic field (Bz(ind)) of the π system show that the three azaborines, in contrast with borazine, sustain much of benzene's π-aromatic character. Among them, 1,3-azaborine exhibits the strongest π delocalization, while 1,4-azaborine is the weakest. Contour maps of Bz(ind) for individual π orbitals reveal that the differentiation of the magnetic response among the three isomers originates from the π-HOMO orbitals, whose magnetic response is governed by rotational allowed transitions to unoccupied orbitals. The low symmetry of azaborines enables a paratropic response from HOMO to unoccupied orbitals excitations, with their magnitude depending on the shape of interacting orbitals. 1,3-Azaborine presents negligible paratropic contributions to Bz(ind) from HOMO to unoccupied orbitals transitions, where 1,2- and 1,4-azaborine present substantial paratropic contributions, which lead to reduced diatropic response. Natural bond orbital (NBO) analysis employing PEPA shows that only the 1,3-azaborine contains π-electron fully delocalized resonance structures.

Enhanced in vitro cytotoxicity and cellular uptake of DNA bases functionalized gold nanoparticles in HeLa cell lines

DNA base functionalized AuNPs are internalized to a greater extent with high cellular uptake and enhanced cytotoxicity in HeLa cells.

Spectroscopic characterization of C-4 substituted 3,5-dichloro-4H-1,2,6-thiadiazines

The effect of C-4 substitution with electron withdrawing groups on the structure and electronic properties of three 3,5-dichloro-4H-1,2,6-thiadiazines is characterized using resonance Raman, absorption and photoluminescence spectroscopies.

A comparative study of the aromaticity of pyrrole, furan, thiophene, and their aza-derivatives

The relative aromaticity of pyrrole, furan, thiophene, and their aza-derivatives has been examined using TRE (topological resonance energy), MRE (magnetic resonance energy), ring current (RC), and ring current diamagnetic susceptibility (χG) methods. The results obtained were compared with results obtained by others who used the energetic method ASE (aromatic stabilization energy), the geometric method HOMA (harmonic oscillator model of aromaticity), and the magnetic method NICS(1) (nucleus-independent chemical shift). The impact of nitrogen atoms on the aromaticity of the aza-derivatives of pyrrole, furan, and thiophene is discussed. An excellent correlation was found between the energetic (TRE, MRE) and magnetic (RC and χG) criteria of aromaticity for all compounds. It was expected that inclusion of a heteroatom would decrease the aromaticity relative to the cyclopentadienyl anion. Our results show that the type of the first heteroatom, which donates two electrons to the system, as well as the number of nitrogen atoms and their positions in the molecule have a strong effect on aromaticity. In general, aromaticity is enhanced when the nitrogen atom is adjacent to the first heteroatom. The magnitude of aromaticity is related closely with the uniformity of distribution of π-electrons in the molecule.

Entropy versus aromaticity in the conformational dynamics of aromatic rings

Comparison of the results of Car-Parrinello molecular dynamics simulations of isolated benzene, pyrimidine and 1,2,4-triazine molecules reveals that the unusually low population of planar geometry of the benzene ring is caused by entropy effects despite its high aromaticity. The decrease in symmetry of the molecule results in smaller changes in entropy and Gibbs free energy due to out-of-plane deformations of the ring, leading to an increase in the population of planar geometry of the ring. This leads to differences in the topology of potential energy and Gibbs free energy surfaces.

Effect of the H-bonding on aromaticity of purine tautomers

Four tautomers of purine (1-H, 3-H, 7-H, and 9-H) and their equilibrium H-bonded complexes with F(-) and HF for acidic and basic centers, respectively, were optimized by means of the B3LYP/6-311++G(d,p) level of theory. Purine tautomer stability increases in the following series: 1-H < 3-H < 7-H < 9-H, consistent with increasing aromaticity. Furthermore, the presence of a hydrogen bond with HF does not change this order. For neutral H-bonded complexes, the strongest and the weakest intermolecular interactions occur (-14.12 and -10.49 kcal/mol) for less stable purine tautomers when the proton acceptor is located in the five- and six-membered rings, respectively. For 9-H and 7-H tautomers the order is reversed. The H-bond energy for the imidazole complex with HF amounts to -14.03 kcal/mol; hence, in the latter case, the fusion of imidazole to pyrimidine decreases its basicity. The ionic H-bonds of N(-)···HF type are stronger by ~10 kcal/mol than the neutral N···HF intermolecular interactions. The hydrogen bond N(-)···HF energies in pyrrole and imidazole are -32.28 and -30.03 kcal/mol, respectively, and are substantially stronger than those observed in purine complexes. The aromaticity of each individual ring and of the whole molecule for all tautomers in ionic complexes is very similar to that observed for the anion of purine. This is not the case for neutral complexes and purine as a reference. The N···HF bonds perturb much more the π-electron structure of five-membered rings than that of the six-membered ones. The H-bonding complexes for 7-H and 9-H tautomers are characterized by higher aromaticity and a much lower range of HOMA variability.