Introduction: aromaticity

Flat building blocks for flat silicene

Silicene is the silicon equivalent of graphene, which is composed of a honeycomb carbon structure with one atom thickness and has attractive characteristics of a perfect two-dimensional π-conjugated sheet. However, unlike flat and highly stable graphene, silicene is relatively sticky and thus unstable due to its puckered or crinkled structure. Flatness is important for stability, and to obtain perfect π-conjugation, electron-donating atoms and molecules should not interact with the π electrons. The structural differences between silicene and graphene result from the differences in their building blocks, flat benzene and chair-form hexasilabenzene. It is crucial to design flat building blocks for silicene with no interactions between the electron donor and π-orbitals. Here, we report the successful design of such building blocks with the aid of density functional theory calculations. Our fundamental concept is to attach substituents that have sp-hybrid orbitals and act as electron donors in a manner that it does not interact with the π orbitals. The honeycomb silicon molecule with BeH at the edge designed according to our concept, clearly shows the same structural, charge distribution and molecular orbital characteristics as the corresponding carbon-based molecule.

Highly-conducting molecular circuits based on antiaromaticity

Aromaticity is a fundamental concept in chemistry. It is described by Hückel’s rule that states that a cyclic planar π-system is aromatic when it shares 4n+2 π-electrons and antiaromatic when it possesses 4n π-electrons. Antiaromatic compounds are predicted to exhibit remarkable charge transport properties and high redox activities. However, it has so far only been possible to measure compounds with reduced aromaticity but not antiaromatic species due to their energetic instability. Here, we address these issues by investigating the single-molecule charge transport properties of a genuinely antiaromatic compound, showing that antiaromaticity results in an order of magnitude increase in conductance compared with the aromatic counterpart. Single-molecule current–voltage measurements and ab initio transport calculations reveal that this results from a reduced energy gap and a frontier molecular resonance closer to the Fermi level in the antiaromatic species. The conductance of the antiaromatic complex is further modulated electrochemically, demonstrating its potential as a high-conductance transistor.

Antiaromatic molecules are predicted to have unusual charge transport properties, but are notoriously unstable and reactive. Here, the authors successfully fabricate an antiaromatic molecular circuit, based on a macrocyclic complex, displaying much higher conductance than its aromatic counterpart.

Electron conjugation versus π-π repulsion in substituted benzenes: why the carbon-nitrogen bond in nitrobenzene is longer than in aniline

Gas-phase electron diffraction experiments show that the C-N bond in aniline (1.407 Å) is significantly shorter than in nitrobenzene (1.486 Å). It is known that the amino group is electron-donating and the nitro group is electron-withdrawing, and both substitution groups can effectively conjugate with benzene. Thus, it is puzzling why the C-N bond in nitrobenzene is even longer than the single C-N bond in methylamine (1.472 Å). In this work, we performed computations by strictly localizing the π electrons with the block-localized wavefunction (BLW) method, which is a variant of ab initio valence bond theory. Geometry optimizations of electron-localized states, where the conjugation over the C-N bond is quenched, show that the conjugation in nitrobenzene is only half of the conjugation in aniline. But even in optimal electron-localized states, the C-N bond in nitrobenzene is still 0.074 Å longer than in aniline. As a consequence, it is indeed not the π conjugation which is responsible for the disparity of the C-N bond distances in these systems. Instead, we demonstrated that the π-π repulsion, which is contributed by both Pauli exchange and electrostatic interaction, plays the key role in this "abnormal" behavior. Notably, the π resonance within the nitro group generates a considerable dipole, which repels the π electrons in the benzene ring. The deactivation of the resonance within the nitro group significantly shortens the C-N bond by 0.06 Å. The unfavorable π-π electrostatic repulsion is further exemplified by N2O4. In fact, the destabilizing π-π repulsion is ubiquitous but largely neglected in conjugated systems where only the stabilizing conjugation is the focus. Experimental phenomena such as the C-N bond distances in aniline and nitrobenzene result from the balance of both stabilizing and destabilizing forces.

Cycloreversion of the CO2 trimer: a paradigmatic pseudopericyclic [2 + 2 + 2] cycloaddition reaction

Very recently, the CO₂ trimer has been experimentally synthesized, isolated and characterized. This process opens new ways for the withdrawal and storage of this greenhouse gas. The trimer is reported to be stable up to -40 °C, with a lifetime of about 40 min at this temperature. At these or under harsher thermal conditions it reverts to the three monomers. The mechanism of this reaction has been theoretically studied and the electronic character of the associated transition state has been analyzed from a variety of perspectives (energetic, magnetic, electron localization and delocalization functions) which indicate that it has paradigmatic pseudopericyclic character. To allow for a comparative study, the isoelectronic fragmentations of cyclohexane into three units of ethylene and of benzene into three units of acetylene have been included in this work. The study of a similar series of formally forbidden-four-centered [2 + 2] cycloreversions confirmed the pseudopericyclic nature of these reactions when the CO₂ dimer or trimer is involved.

CC/B–N substitution in five membered heterocycles. A computational analysis

Novel five-membered azaboroles are aromatic, stable under neutral conditions, isomer stabilization energy is explained using σ-bond and aromatic stabilization energies.

Applications of Systematic Molecular Scaffold Enumeration to Enrich Structure-Activity Relationship Information

Establishing structure–activity relationships (SARs) in hit identification during early stage drug discovery is important in accelerating hit confirmation and expansion. We describe the development of EnCore, a systematic molecular scaffold enumeration protocol using single atom mutations, to enhance the application of objective scaffold definitions and to enrich SAR information from analysis of high-throughput screening output. A list of 43 literature medicinal chemistry compound series, each containing a minimum of 100 compounds, published in the Journal of Medicinal Chemistry was collated to validate the protocol. Analysis using the top representative Level 1 scaffolds this list of literature compound series demonstrated that EnCore could mimic the scaffold exploration conducted when establishing SAR. When EnCore was applied to analyze an HTS library containing over 200 000 compounds, we observed that over 70% of the molecular scaffolds matched extant scaffolds within the library after enumeration. In particular, over 60% of the singleton scaffolds with only one representative compound were found to have structurally related compounds after enumeration. These results illustrate the potential of EnCore to enrich SAR information. A case study using literature cyclooxygenase-2 inhibitors further demonstrates the advantage of EnCore application in establishing SAR from structurally related scaffolds. EnCore complements literature enumeration methods in enabling changes to the physicochemical properties of molecular scaffolds and structural modifications to aliphatic rings and linkers. The enumerated scaffold clusters generated would constitute a comprehensive collection of scaffolds for scaffold morphing and hopping.

Induction of unidirectional π-electron rotations in low-symmetry aromatic ring molecules using two linearly polarized stationary lasers

A new laser-control scenario of unidirectional π-electron rotations in a low-symmetry aromatic ring molecule having no degenerate excited states is proposed. This scenario is based on dynamic Stark shifts of two relevant excited states using two linearly polarized stationary lasers. Each laser is set to selectively interact with one of the two electronic states, the lower and higher excited states are shifted up and down with the same rate, respectively, and the two excited states become degenerate at their midpoint. One of the four control parameters of the two lasers, i.e. two frequencies and two intensities, determines the values of all the other parameters. The direction of π-electron rotations, clockwise or counter-clockwise rotation, depends on the sign of the relative phase of the two lasers at the initial time. An analytical expression for the time-dependent expectation value of the rotational angular momentum operator is derived using the rotating wave approximation (RWA). The control scenario depends on the initial condition of the electronic states. The control scenario with the ground state as the initial condition was applied to toluene molecules. The derived time-dependent angular momentum consists of a train of unidirectional angular momentum pulses. The validity of the RWA was checked by numerically solving the time-dependent Schrödinger equation. The simulation results suggest an experimental realization of the induction of unidirectional π-electron rotations in low-symmetry aromatic ring molecules without using any intricate quantum-optimal control procedure. This may open up an effective generation method of ring currents and current-induced magnetic fields in biomolecules such as amino acids having aromatic ring molecules for searching their interactions.

Pnictogen-Silicon Analogues of Benzene

Since the discovery of the first "inorganic benzene" (borazine, B3N3H6), the synthesis of other noncarbon derivatives is an ongoing challenge in Inorganic Chemistry. Here we report on the synthesis of the first pnictogen-silicon congeners of benzene, the triarsa- and the triphospha-trisilabenzene [(PhC(NtBu)2)3Si3E3] (E = P (1a), As (1b)) by a simple metathesis reaction. These compounds are formed by the reaction of [Cp″2Zr(η(1:1)-E4)] (E = P, As; Cp″ = C5H3tBu2) with [PhC(NtBu)2SiCl] in toluene at room temperature along with the silicon pnictogen congeners of the cyclobutadiene, [(PhC(NtBu)2)2Si2E2] (E = P (2a), As (2b)), which is unprecedented for the arsenic system 2b. All compounds were comprehensively characterized, and density functional theory calculations were performed to verify the stability and the aromatic character of the triarsa- and the triphospha-trisilabenzene.

Innocent BN bond substitution in anthracene derivatives

Extended azaborine heterocycles are promising biomedical and electronic materials. Herein we report the synthesis of a novel family of azaborine anthracene derivatives and their structural, electrochemical and spectroscopic characterization. We observe that the properties of these materials are remarkably similar to the parent hydrocarbons, suggesting the innocence of the CC to BN bond substitution. Our results support the prospective stability to long-term usage of extended azaborines and the feasibility of using such materials in device applications.

Aromaticity of the doubly charged [8]circulenes

Magnetically induced current densities and current pathways have been calculated for a series of fully annelated dicationic and dianionic tetraphenylenes, which are also named [8]circulenes.

Search for aromatic anions in the P2E3− (E = N, P, As, Sb, Bi) series

Nine aromatic five-membered rings were computationally characterized in the P₂E₃⁻ (E = N, P, As, Sb, Bi) series of anions.

Topological definition of ring currents

A definition of ring currents in a velocity vector field is proposed according to topological criteria: ring currents are axial vortices confined in, or rotating beyond, a separatrix, i.e., the boundary which marks the limits of the vortex.

Understanding conductivity in molecular switches: a real space approach in octaphyrins

In recent years, expanded porphyrins have emerged as a promising class of π-conjugated switches whose conductance is studied from the electron density.

Understanding the molecular switching properties of octaphyrins

Triggering Hückel–Möbius topological and aromaticity switches in octaphyrins by protonation and redox reactions.

Celebrating the 150th anniversary of the Kekulé benzene structure

This themed issue presents a collection of articles on ‘Electron delocalization and aromaticity: Celebrating the 150th Anniversary of the Kekulé Benzene Structure’.

Low bandgap semiconducting polymers for polymeric photovoltaics

This review highlights the design rules for low bandgap semiconducting polymers, with the overview of their applications in polymer solar cells and polymer photodetectors.

A Study of π-π Stacking Interactions and Aromaticity in Polycyclic Aromatic Hydrocarbon/Nucleobase Complexes

We analysed the interactions and aromaticity electron-density delocalisation observed in π-π complexes between the phenalenyl radical and acenaphthylene, and the DNA and RNA nucleobases (adenine, guanine, cytosine, thymine and uracil). Interaction energies are obtained at the M06-2X/6-311++G(2df,p) computational level for gas phase and PCM-water conditions. For both the phenalenyl radical and acenaphthylene, the complexes formed with guanine are the most stable ones. Atoms in molecules and natural bond orbital results reveal weak π-π interactions between both interacting moieties, characterized by bond critical points between C⋅⋅⋅C and C⋅⋅⋅N atoms. Nucleus independent chemical shifts (NICS) indicate the retention of the aromatic character of the monomers in the outer region of the complex. The fluctuation indexes reveal a loss of electron delocalisation upon complexation for all cases except guanine complexes. Nevertheless, the interface region shows large negative NICS values, which is not associated with an increase of the aromaticity or electron-density delocalisation, but with magnetic couplings of both molecules, leading to an unrealistic description of the aromatic behaviour in that region.

An Aromaticity Scale Based on the Topological Analysis of the Electron Localization Function Including σ and π Contributions

In this work, the average bifurcation value of the electron localization function (ELF) of both σ (ELFσ) and π (ELFπ) contributions was used to construct an aromaticity scale for chemical compounds. We have validated the scale with a series of well-known molecules and then used it to evaluate global aromaticity on aluminum based clusters, which present σ aromaticity and π antiaromaticity. The proposed scaled predicts an overall antiaromatic character for the Al4(4)(-) moiety.

First-Principles Prediction of Enthalpies of Formation for Polycyclic Aromatic Hydrocarbons and Derivatives

In this article, the first-principles prediction of enthalpies of formation is demonstrated for 669 polycyclic aromatic hydrocarbon (PAH) compounds and a number of related functionalized molecules. It is shown that by extrapolating density functional theory calculations to a large basis set limit and then applying a group based correction scheme that good results may be obtained. Specifically, a mean unsigned deviation and root mean squared deviation from the experimental enthalpies of formation data of 5.0 and 6.4 kJ/mol, respectively, are obtained using this scheme. This computational scheme is economical to compute and straightforward to apply, while yielding results of reasonable reliability. The results are also compared for a smaller set of molecules to the predictions given by the G3B3 and G3MP2B3 variants of the Gaussian-3 model chemistry with a mean unsigned deviation and root mean squared deviation from the experimental enthalpies of formation of 4.5 and 4.8 kJ/mol, respectively.

Triplet State Aromaticity: NICS Criterion, Hyperconjugation, and Charge Effects

Aromaticity, one of the most important concepts in organic chemistry, has attracted considerable interest from both experimentalists and theoreticians. It remains unclear which NICS index is best to evaluate the triplet-state aromaticity. Here, we carry out thorough density functional theory (DFT) calculations to examine this issue. Our results indicate that among the various computationally available NICS indices, NICS(1)zz is the best for the triplet state. The correlations can be improved from 0.840 to 0.938 when only neutral species are considered, demonstrating the significant effect of the charge on the triplet-state aromaticity. In addition, calculations suggest that five-membered cyclic species with "hyperconjugative" aromaticity (and antiaromaticity) in the S0 state will become antiaromatic (and aromatic) in the T1 state, indicating an important role of hyperconjugation. Finally, a moderate correlation (r(2) =0.708) is identified between the NICS(1)zz values and spin distributions.

Linearly Polymerized Benzene Arrays As Intermediates, Tracing Pathways to Carbon Nanothreads

How might fully saturated benzene polymers of composition [(CH)6]n form under high pressure? In the first approach to answering this question, we examine the stepwise increase in saturation of a one-dimensional stack of benzene molecules by enumerating the partially saturated polymer intermediates, subject to constraints of unit cell size and energy. Defining the number of four-coordinate carbon atoms per benzene formula unit as the degree of saturation, a set of isomers for degree-two and degree-four polymers can be generated by either thinking of the propagation of partially saturated building blocks or by considering a sequence of cycloadditions. There is also one 4 + 2 reaction sequence that jumps directly from a benzene stack to a degree-four polymer. The set of degree-two polymers provides several useful signposts toward achieving full saturation: chiral versus achiral building blocks, certain forms of conformational freedom, and also dead ends to further saturation. These insights allow us to generate a larger set of degree-four polymers and enumerate the many pathways that lead from benzene stacks to completely saturated carbon nanothreads.