Helical polycyclic hydrocarbons with open-shell singlet ground states and ambipolar redox behaviors

Organic π-conjugated polycyclic hydrocarbons (PHs) with an open-shell diradical character are attracting increasing interest due to their promising applications in organic electronics and spintronics. However, most of the open-shell PHs synthesized thus far are based on planar π-conjugated molecules. Herein, we report the synthesis and characterization of two new quinodimethane-embedded expanded helicenes H1 and H2. The helical structures of both molecules were revealed using X-ray crystallographic analysis. It was elucidated in detailed experimental and theoretical studies that they possess an open-shell singlet biradical structure in the ground state and show a small energy gap and amphoteric redox behavior. Both compounds can also be easily oxidized or reduced into relatively stable charged species. The dianions of H1 and H2 exhibit similar electronic structures to the respective isoelectronic structures of their all-benzenoid helical analogues according to NMR measurements and theoretical calculations. Moreover, the structures of the dication and dianion of H2 were identified by X-ray crystallographic analysis, revealing the effect of electron transfer on their backbones and aromaticity. This study thus opens up new avenues for both helical polycyclic π-systems and diradicaloids.

New Perspectives on Delocalization Pathways in Aromatic Molecular Chameleons

This study comprehensively analyzes the magnetically induced current density of polycyclic compounds labeled as "aromatic chameleons" since they can arrange their π-electrons to exhibit aromaticity in both the ground and the lowest triplet state. These compounds comprise benzenoid moieties fused to a central skeleton with 4n π-electrons and traditional magnetic descriptors are biased due to the superposition of local magnetic responses. In the S0 state, our analysis reveals that the molecular constituent fragments preserve their (anti)aromatic features in agreement with two types of resonant structures: one associated with aromatic benzenoids and the other with a central antiaromatic ring. Regarding the T1 state, a global and diatropic ring current is revealed. Our aromaticity study is complemented with advanced electronic and geometric descriptors to consider different aspects of aromaticity, particularly important in the evaluation of excited state aromaticity. Remarkably, these descriptors consistently align with the general features on the main delocalization pathways in polycyclic hydrocarbons consisting of fused 4n π-electron rings. Moreover, our study demonstrates an inverse correlation between the singlet-triplet energy difference and the antiaromatic character of the central ring in S0.

On the aromaticity and stability of benzynes in the ground and lowest-lying triplet excited states

In this work, we have revisited the aromaticity of benzyne isomers at the unrestricted density functional theory level (UDFT) using the energetic, magnetic, and delocalization criteria. In addition, this last criterion has also been analyzed employing complete active space (CASSCF) calculations. The results show conservation of aromaticity in these monocycles. Additionally it is observed that this trend is maintained in polycyclic aromatic hydrocarbon derivatives such as biradical didehydrophenanthrenes. Do these results imply a violation of Baird's rule? The answer is No, because this conservation in aromaticity is due to the loss of hydrogen atoms affects only the electronic σ skeleton and exerts a minor influence on the π cloud. Additionally, we have analyzed the relative stability of benzyne isomers and their relationship with experimental ΔES-T values. According to the literature, the stability of the benzynes in the singlet state is due to an effective interaction between the electrons of the biradical centers; however, this effect is completely reversed in the triplet state, which explains why the para isomer has the lowest ΔES-T gap.

Aromaticity and Magnetic Behavior in Benzenoids: Unraveling Ring Current Combinations

Nowadays, an active research topic is the connection between Clar's rule, aromaticity, and magnetic properties of polycyclic benzenoid hydrocarbons. In the present work, we employ a meticulous magnetically induced current density analysis to define the net current flowing through any cyclic circuit, connecting it to aromaticity based on the ring current concept. Our investigation reveals that the analyzed polycyclic systems display a prominent global ring current, contrasting with subdued semi-local and local ring currents. These patterns align with Clar's aromatic π-sextets only in cases where migrating π-sextet structures are invoked. The results of this study will enrich our comprehension of aromaticity and magnetic behavior in such systems, offering significant insights into coexisting ring current circuits in these systems.

Unlocking Heteroaromatic Ring Systems through Chalcogen Insertion into Boroles

In this work, we report the reactivity of various annulated borole derivatives toward chalcogen (O, S, and Se) insertion. Among a series of 9-borafluorenes with different boron substituents (Ph, Br, or o-carboranyl) and a mixed thiophene-benzene-fused derivative, only the 9-o-carboranyl-substituted 9-borafluorene yielded the complete set of chalcogen-containing heteroarenes, including the first 1,2-selenaborinine derivative. To evaluate the aromaticity of this heterocyclic analogue of phenanthrene, nucleus-independent chemical shift (NICS) values were computed and compared to those of its lighter group 16 congeners.

Exploring the influence of graphene on antiaromaticity of pentalene

Theoretical investigations on the influence of graphene fragments on the antiaromaticity of pentalene are conducted by employing multiple aromaticity descriptors based on magnetic, geometric and electronic criteria. NICS as a sole descriptor for analysing the antiaromaticity of pentalene on graphene fragments has to be carefully considered while looking through the other aromaticity indicators.

Gad S, M. El-Mehasseb I, E. El-Kelany K. Isosterism in pyrrole via azaboroles substitution, a theoretical investigation for electronic structural, stability and aromaticity

This work uses ab-initio CBS-QB3 and density functional theory (B3LYP) to analyze the structure, stability, and aromaticity of all isosteric nitrogen-boron pyrroles. The mono-NB unit substituted group of the isosteric NB pyrrole has four isosteres, whereas the multi-NB unit substituted group has two isosteres. These two groups make up all isosteric NB pyrrole. For structural, energetic, magnetic, and electron delocalization criteria, the results highlight the predominance of the PN3B2 isostere and its greater stability over other conformers. In addition, the global reactivity indices, ESP, HOMO-LUMO, and NBO charges have all been estimated to forecast the active side's electron donation and acceptance. These isosteres are categorized as weak electrophiles and marginal nucleophiles. NB-isosteres have poorer stability, HOMO-LUMO gap, and aromaticity than the parent (pyrrole). In general, NB compounds with more ring sharing are less aromatic than NB molecules with less ring sharing. The current study is anticipated to help in understanding of the chemistry of NB substituted molecules and their experimental identification and characterization.

Exploration of the potential energy surface in mixed Zintl clusters applying an automatic Johnson polyhedra generator: the case of arachno E6M24- (E = Si, Ge, Sn; M = Sb, Bi)

A new algorithm called Automatic Johnson Cluster Generator (AJCG) is presented, which, as its name indicates, allows the definition of the desired Johnson polyhedron to subsequently carry out all the possible permutations between the atoms that form this polyhedron. This new algorithm allows the exhaustive study of the structures' potential energy surface (PES). In addition, the AJCG algorithm is helpful for the study of three-dimensional compounds such as boranes or Zintl clusters and their structural derivatives with two or more different atoms. The automatic filling of vertices is particularly useful in mixed compounds because of the possibility of taking into account all possible configurations in the structure. As a test system, we investigated the arachno-type E6M24- (E = Si, Ge, Sn; M = Sb, Bi) structure which has eight vertices and complies with Wade-Mingos rules. Initially, we defined a bipyramidal structure (10 vertices), and filled the vertices with the atoms in all possible configurations. Since the selected system has eight atoms, the two remaining vertices were filled with pseudo atoms to complete the structure. After re-optimizing the initial population generated with AJCG, a large number of isomers with energy below 10 kcal mol-1 are identified. These results show that the most stable isomers possess homonuclear M-M bonds, except Sn6Bi24-. Although the overall putative minima differ at the PBE0-D3 and DLPNO-CCSD(T) levels, they are always competitive minima. In addition to using high-precision methodologies to correctly study relative energies, applying solvent effects in highly charged systems becomes mandatory. The aromatic character of these studied systems was demonstrated qualitatively with two- and three-dimensional mapping and quantitatively by calculating the value of the z-component of the induced magnetic field at the cage center, including scalar and spin-orbit correction for relativistic effects. The compounds studied have a high degree of aromaticity, which allows us to establish that despite structural modifications (i.e., from closo to arachno), the aromaticity is preserved.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Planar Elongated B12 Structure in M3B12 Clusters (M = Cu-Au)

Here, it is shown that the M₃B₁₂ (M = Cu-Au) clusters' global minima consist of an elongated planar B₁₂ fragment connected by an in-plane linear M₃ fragment. This result is striking since this B₁₂ planar structure is not favored in the bare cluster, nor when one or two metals are added. The minimum energy structures were revealed by screening the potential energy surface using genetic algorithms and density functional theory calculations. Chemical bonding analysis shows that the strong electrostatic interactions with the metal compensate for the high energy spent in the M₃ and B₁₂ fragment distortion. Furthermore, metals participate in the delocalized π-bonds, which infers an aromatic character to these species.

Structural Factors That Determine the Activity of the Xenobiotic Reductase B Enzyme from Pseudomonas putida on Nitroaromatic Compounds

Xenobiotic reductase B (XenB) catalyzes the reduction of the aromatic ring or nitro groups of nitroaromatic compounds with methyl, amino or hydroxyl radicals. This reaction is of biotechnological interest for bioremediation, the reuse of industrial waste or the activation of prodrugs. However, the structural factors that explain the binding of XenB to different substrates are unknown. Molecular dynamics simulations and quantum mechanical calculations were performed to identify the residues involved in the formation and stabilization of the enzyme/substrate complex and to explain the use of different substrates by this enzyme. Our results show that Tyr65 and Tyr335 residues stabilize the ligands through hydrophobic interactions mediated by the aromatic rings of these aminoacids. The higher XenB activity determined with the substrates 1,3,5-trinitrobenzene and 2,4,6-trinitrotoluene is consistent with the lower energy of the highest occupied molecular orbital (LUMO) orbitals and a lower energy of the homo orbital (LUMO), which favors electrophile and nucleophilic activity, respectively. The electrostatic potential maps of these compounds suggest that the bonding requires a large hydrophobic region in the aromatic ring, which is promoted by substituents in ortho and para positions. These results are consistent with experimental data and could be used to propose point mutations that allow this enzyme to process new molecules of biotechnological interest.

Spherical Aromaticity of Tetrahedral Pnictogens through Off-Nucleus Isotropic Magnetic Shielding

This work revealed the spherical aromaticity of some inorganic E₄ cages and their protonated E₄ H⁺ ions (E=N, P, As, Sb, and Bi). For this purpose, we employed several evaluations like (0D-1D) nucleus independent chemical shift (NICS), multidimensional (2D-3D) off-nucleus isotropic shielding σ_iso (r), and natural bond orbital (NBO) analysis. The magnetic calculations involved gauge-including atomic orbitals (GIAO) with two density functionals B3LYP and WB97XD, and basis sets of Jorge-ATZP, 6-311+G(d,p), and Lanl2DZp. The Jorge-ATZP basis set showed the best consistency. Our findings disclosed non-classical aromatic characters in the above molecules, which decreased from N to Bi cages. Also, the results showed more aromaticity in E₄ than E₄ H⁺ . The NBO analysis attributed the aromaticity in the above molecules to the residual density of the overlapping σ-bonding orbitals. So, the aromaticity in these molecules is unlike the classical aromaticity that is associated with electron delocalization. Scanning 1D σ_iso (r) variation along E-E bonds indicated a lowering in the shielding trend from N to Bi cages. The 3D results showed a similar decrease in the relative volumetric diffusion of the magnetic activity, whereas the volumetric ratio of V_1ppm /V_2ppm is almost constant for all the E₄ cages.

Aromaticity and Extrusion of Benzenoids Linked to [o-COSAN]- : Clar Has the Answer

Benzene and pyrene can be synthetically linked to [o-COSAN]^- keeping their aromaticity. In contrast, naphthalene and anthracene are extruded in the same reaction. We have proven that extrusion is only favorable if the number of Clar's π-sextets remains constant. Thus, Clar has the answer to whether an attached polycyclic aromatic hydrocarbon to [o-COSAN]^- is extruded or not.

Structures, Energetics, and Spectra of (NH) and (OH) Tautomers of 2-(2-Hydroxyphenyl)-1-azaazulene: A Density Functional Theory/Time-Dependent Density Functional Theory Study

Tautomerization of 2-(2-hydroxyphenyl)-1-azaazulene (2OHPhAZ) in the gas phase and ethanol has been studied using B3LYP, M06-2X, and ωB97XD density functional theory (DFT) with different basis sets. For more accurate data, energies were refined at CCSD(T)/6-311++G(2d,2p) in the gas phase. Nuclear magnetic resonance (NMR), aromaticity, Fukui functions, acidity, and basicity were also calculated and compared with experimental data. Time-dependent density functional theory (TDDFT)-solvation model based on density (TDDFT-SMD) calculations in acetonitrile have been utilized for the simulation of UV-vis electronic spectra. In addition, electronic structures of the investigated system have been discussed. The results reveal that the enol form (2OHPhAZ) is thermodynamically and kinetically stable relative to the keto tautomer (2OPhAZ) and different rotamers (2OHPhAZ-R1:R3) in the gas phase and ethanol. A comparison with the experiment illustrates a good agreement and supports the computational findings.

Persistent Planar Tetracoordinate Carbon in Global Minima Structures of Silicon-Carbon Clusters

Recently, we reported a series of global minima whose structures consist of carbon rings decorated with heavier group 14 elements. Interestingly, these structures feature planar tetracoordinate carbons (ptCs) and result from the replacement of five or six protons (H+) from the cyclopentadienyl anion (C5H5−) or the pentalene dianion (C8H62−) by three or four E2+ dications (E = Si–Pb), respectively. The silicon derivatives of these series are the Si3C5 and Si4C8 clusters. Here we show that ptC persists in some clusters with an equivalent number of C and Si atoms, i.e., Si5C5, Si8C8, and Si9C9. In all these species, the ptC is embedded in a pentagonal C5 ring and participates in a three-center, two-electron (3c-2e) Si-ptC-Si σ-bond. Furthermore, these clusters are π-aromatic species according to chemical bonding analysis and magnetic criteria.

Simple and efficient visualization of aromaticity: bond currents calculated from NICS values

Aromaticity is a fundamental concept in chemistry, underpinning the properties and reactivity of many organic compounds and materials. The ability to easily and accurately discern aromatic behavior is key to leveraging it as a design element, yet most aromaticity metrics struggle to combine accurate quantitative evaluation, intuitive interpretability, and user-friendliness. We introduce a new method, NICS2BC, which uses simple and inexpensive NICS calculations to generate information-rich and easily-interpreted bond-current graphs. We test the quantitative and qualitative characterizations afforded by NICS2BC for a selection of molecules of varying structural and electronic complexity, to demonstrate its accuracy and ease of analysis. Moreover, we show that NICS2BC successfully identifies ring-current patterns in molecules known to be difficult cases to interpret with NICS and enables deeper understanding of local aromaticity trends, demonstrating that our method adds additional insight.

Analysis of Local and Global Aromaticity in Si3C5 and Si4C8 Clusters. Aromatic Species Containing Planar Tetracoordinate Carbon

The minimum energy structures of the Si3C5 and Si4C8 clusters are planar and contain planar tetracoordinate carbons (ptCs). These species have been classified, qualitatively, as global (π) and local (σ) aromatics according to the adaptive natural density partitioning (AdNDP) method, which is an orbital localization method. This work evaluates these species’ aromaticity, focusing on confirming and quantifying their global and local aromatic character. For this purpose, we use an orbital localization method based on the partitioning of the molecular space according to the topology of the electronic localization function (LOC-ELF). In addition, the magnetically induced current density is analyzed. The LOC-ELF-based analysis coincides with the AdNDP study (double aromaticity, global, and local). Moreover, the current density analysis detects global and local ring currents. The strength of the global and local current circuit is significant, involving 4n + 2 π- and σ-electrons, respectively. The latter implicates the Si-ptC-Si fragment, which would be related to the 3c-2e σ-bond detected by the orbital localization methods in this fragment.