COMPAS-3: a dataset of peri-condensed polybenzenoid hydrocarbons.

We introduce the third installment of the COMPAS Project - a COMputational database of Polycyclic Aromatic Systems, focused on peri-condensed polybenzenoid hydrocarbons. In this installment, we develop two datasets containing the optimized ground-state structures and a selection of molecular properties of ∼39k and ∼9k peri-condensed polybenzenoid hydrocarbons (at the GFN2-xTB and CAM-B3LYP-D3BJ/cc-pvdz//CAM-B3LYP-D3BJ/def2-SVP levels, respectively). The manuscript details the enumeration and data generation processes and describes the information available within the datasets. An in-depth comparison between the two types of computation is performed, and it is found that the geometrical disagreement is maximal for slightly-distorted molecules. In addition, a data-driven analysis of the structure-property trends of peri-condensed PBHs is performed, highlighting the effect of the size of peri-condensed islands and linearly annulated rings on the HOMO-LUMO gap. The insights described herein are important for rational design of novel functional aromatic molecules for use in, e.g., organic electronics. The generated datasets provide a basis for additional data-driven machine- and deep-learning studies in chemistry.

Interpretable Deep-Learning Unveils Structure-Property Relationships in Polybenzenoid Hydrocarbons.

In this work, interpretable deep learning was used to identify structure-property relationships governing the HOMO-LUMO gap and the relative stability of polybenzenoid hydrocarbons (PBHs) using a ring-based graph representation. This representation was combined with a subunit-based perception of PBHs, allowing chemical insights to be presented in terms of intuitive and simple structural motifs. The resulting insights agree with conventional organic chemistry knowledge and electronic structure-based analyses and also reveal new behaviors and identify influential structural motifs. In particular, we evaluated and compared the effects of linear, angular, and branching motifs on these two molecular properties and explored the role of dispersion in mitigating the torsional strain inherent in nonplanar PBHs. Hence, the observed regularities and the proposed analysis contribute to a deeper understanding of the behavior of PBHs and form the foundation for design strategies for new functional PBHs.

A Crowning Achievement: The First Solution-Phase Synthesis of Circumcoronenes.

Circumcoronene is a highly symmetric polybenzenoid hydrocarbon that has fascinated organic chemists and materials scientists for decades. However, until recently, it had only been studied theoretically, because its synthesis in solution remained challenging. In a recent report, Wu and co-workers described the first successful solution-phase synthesis and isolation of crystalline circumcoronene derivatives, validating the long-predicted Clar structure. Their synthetic approach enables preparation of new polycyclic benzenoid hydrocarbons with multiple K-regions, which are difficult to obtain, and expands our understanding of the chemistry of these systems.

Localized Antiaromaticity Hotspot Drives Reductive Dehydrogenative Cyclizations in Bis- and Mono-Helicenes.

We describe reductive dehydrogenative cyclizations that form hepta-, nona-, and decacyclic anionic graphene subunits from mono- and bis-helicenes with an embedded five-membered ring. The reaction of bis-helicenes can either proceed to the full double annulation or be interrupted by addition of molecular oxygen at an intermediate stage. The regioselectivity of the interrupted cyclization cascade for bis-helicenes confirms that relief of antiaromaticity is a dominant force for these facile ring closures. Computational analysis reveals the unique role of the preexisting negatively charged cyclopentadienyl moiety in directing the second negative charge at a specific remote location and, thus, creating a localized antiaromatic region. This region is the hotspot that promotes the initial cyclization. Computational studies, including MO analysis, molecular electrostatic potential maps, and NICS(1.7)ZZ calculations, evaluate the interplay of the various effects including charge delocalization, helicene strain release, and antiaromaticity. The role of antiaromaticity relief is further supported by efficient reductive closure of the less strained monohelicenes where the relief of antiaromaticity promotes the cyclization even when the strain is substantially reduced. The latter finding significantly expands the scope of this reductive alternative to the Scholl ring closure.

The COMPAS Project: A Computational Database of Polycyclic Aromatic Systems. Phase 1: cata-Condensed Polybenzenoid Hydrocarbons.

Chemical databases are an essential tool for data-driven investigation of structure-property relationships and for the design of novel functional compounds. We introduce the first phase of the COMPAS Project─a COMputational database of Polycyclic Aromatic Systems. In this phase, we developed two data sets containing the optimized ground-state structures and a selection of molecular properties of ∼34k and ∼9k cata-condensed polybenzenoid hydrocarbons (at the GFN2-xTB and B3LYP-D3BJ/def2-SVP levels, respectively) and placed them in the public domain. Herein, we describe the process of the data set generation, detail the information available within the data sets, and show the fundamental features of the generated data. We analyze the correlation between the two types of computations as well as the structure-property relationships of the calculated species. The data and insights gained from them can inform rational design of novel functional aromatic molecules for use in, e.g., organic electronics, and can provide a basis for additional data-driven machine- and deep-learning studies in chemistry.

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.

Peptide-Metal Frameworks with Metal Strings Guided by Dispersion Interactions.

Despite impressive advances in the construction of metal-organic frameworks (MOFs), the formation of networks from peptidic ligands is difficult, though they are sought after for their modularity and biocompatibility. Herein we present a peptide-metal framework that consists of helical oligoproline ligands and Zn/K (or Zn/Rb). The crystalline network contains pleated nanosheets with the metal ions aligned in strings. This unprecedented architecture derives from under-appreciated London dispersion interactions between the oligoproline ligands that play in concert with the metal coordination to create the network. Hence, the secondary structure of the peptidic ligand represents an additional control element for the creation of new MOF architectures. We anticipate that our results will instruct the design of further peptidic MOFs and enable the generation of versatile biocompatible materials.

Nanoscale Phase Segregation in Supramolecular π-Templating for Hybrid Perovskite Photovoltaics from NMR Crystallography.

The use of layered perovskites is an important strategy to improve the stability of hybrid perovskite materials and their optoelectronic devices. However, tailoring their properties requires accurate structure determination at the atomic scale, which is a challenge for conventional diffraction-based techniques. We demonstrate the use of nuclear magnetic resonance (NMR) crystallography in determining the structure of layered hybrid perovskites for a mixed-spacer model composed of 2-phenylethylammonium (PEA+) and 2-(perfluorophenyl)ethylammonium (FEA+) moieties, revealing nanoscale phase segregation. Moreover, we illustrate the application of this structure in perovskite solar cells with power conversion efficiencies that exceed 21%, accompanied by enhanced operational stability.

Prediction of Spin Density, Baird-Antiaromaticity, and Singlet-Triplet Energy Gap in Triplet-State Polybenzenoid Systems from Simple Structural Motifs.

Triplet-state aromaticity has been recently proposed as a strategy for designing functional organic electronic compounds, many of which are polycyclic aromatic systems. However, in many cases, the aromatic nature of the triplet state cannot be easily predicted. Moreover, it is often unclear how specific structural manipulations affect the electronic properties of the excited-state compounds. Herein, the relationship between the structure of polybenzenoid hydrocarbons (PBHs) and their spin-density distribution and aromatic character in the first triplet excited state is studied. Although a direct link is not immediately visible, classifying the PBHs according to their annulation sequence reveals regularities. Based on these, a set of guidelines is defined to qualitatively predict the location of spin and paratropicity and the singlet-triplet energy gap in larger PBHs, using only their smaller tri- and tetracyclic components, and subsequently tested on larger systems.

Tuning Magnetic Interactions Between Triphenylene Radicals by Variation of Crystal Packing in Structures with Alkali Metal Counterions.

The monoanion of triphenylene (C18H12, 1) was generated in THF using several alkali metals (Na, K, Rb, and Cs) as reducing agents and crystallized with the corresponding cations in the presence of 18-crown-6 ether. The UV-vis spectroscopy points to the metal-dependent coordination environment of the triphenylene monoanion-radicals, 1·-, in solution. The X-ray diffraction characterization confirmed the formation of a solvent-separated ion pair (SSIP) with sodium ions, [{Na+(18-crown-6)(THF)2}(1·-)] (2), and three contact-ion pair (CIP) complexes formed by larger alkali metal ions, [{K+(18-crown-6)}(1·-)] (3), [{Rb+(18-crown-6)}(1·-)] (4), and [{Cs+(18-crown-6)}(1·-)] (5). Structural analysis of the series reveals a notable geometry perturbation of the triphenylene framework in 2 caused by one-electron acquisition, which is further enhanced by direct metal binding in 3-5. This has been correlated with the aromaticity changes and charge redistribution upon one-electron reduction of 1, as revealed by the computational studies. The EPR spectroscopy and magnetic susceptibility measurements confirm antiferromagnetic interactions corresponding to an S = 1/2 system in the solid state. The magnetic behavior of 3-5 correlates with the arrangement of triphenylene radicals in the crystal structures. All three compounds exhibit antiferromagnetic (AFM) interactions between S = 1/2 radicals in the solid state, but the exchange coupling in 4 and 5 is notably stronger than that in 3, which leads to AFM ordering at 3.8 K in 4 and at 2.0 K in 5. The magnetic phase transitions in 4 and 5 can be interpreted as originating from interactions between the chains of the AFM-coupled S = 1/2 radicals.

Extensive Redox Non-Innocence in Iron Bipyridine-Diimine Complexes: a Combined Spectroscopic and Computational Study.

Metal-ligand cooperation is an important aspect in earth-abundant metal catalysis. Utilizing ligands as electron reservoirs to supplement the redox chemistry of the metal has resulted in many new exciting discoveries. Here, we demonstrate that iron bipyridine-diimine (BDI) complexes exhibit an extensive electron-transfer series that spans a total of five oxidation states, ranging from the trication [Fe(BDI)]3+ to the monoanion [Fe(BDI]-1. Structural characterization by X-ray crystallography revealed the multifaceted redox noninnocence of the BDI ligand, while spectroscopic (e.g., 57Fe Mössbauer and EPR spectroscopy) and computational studies were employed to elucidate the electronic structure of the isolated complexes, which are further discussed in this report.

Predicting bond-currents in polybenzenoid hydrocarbons with an additivity scheme.

We report on the construction and application of a new bond-current additivity scheme for polybenzenoid hydrocarbons. The method is based on identification of the smaller substructures contained in the system, up to tricyclic subunits. Thus, it enables the prediction of any cata-condensed unbranched polybenzenoid hydrocarbon, using a library consisting of only four building blocks. The predicted bond-currents can then be used to generate Nucleus Independent Chemical Shift (NICS) values, the results of which validate previous observations of additivity with NICS-XY-Scans. The limitations of the method are probed, leading to clearly delineated and apparently constant error boundaries, which are independent of the molecular size. It is shown that there is a relationship between the accuracy of the predictions and the molecular structure and specific motifs that are especially challenging are identified. The results of the additivity method, combined with the transparent description of its strengths and weaknesses, ensure that this method can be used with well-defined reliability for characterization of polybenzenoid hydrocarbons. The resource-efficient and rapid nature of the method makes it a promising tool for screening and molecular design.

Supramolecular Modulation of Hybrid Perovskite Solar Cells via Bifunctional Halogen Bonding Revealed by Two-Dimensional 19F Solid-State NMR Spectroscopy.

There has been an ongoing effort to overcome the limitations associated with the stability of hybrid organic-inorganic perovskite solar cells by using different organic agents as additives to the perovskite formulations. The functionality of organic additives has been predominantly limited to exploiting hydrogen-bonding interactions, while the relevant atomic-level binding modes remain elusive. Herein, we introduce a bifunctional supramolecular modulator, 1,2,4,5-tetrafluoro-3,6-diiodobenzene, which interacts with the surface of the triple-cation double-halide perovskite material via halogen bonding. We elucidate its binding mode using two-dimensional solid-state 19F NMR spectroscopy in conjunction with density functional theory calculations. As a result, we demonstrate a stability enhancement of the perovskite solar cells upon supramolecular modulation, without compromising the photovoltaic performances.

Negative Charge as a Lens for Concentrating Antiaromaticity: Using a Pentagonal "Defect" and Helicene Strain for Cyclizations.

Incorporation of a five-membered ring into a helicene framework disrupts aromatic conjugation and provides a site for selective deprotonation. The deprotonation creates an anionic cyclopentadienyl unit, switches on conjugation, leads to a >200 nm red-shift in the absorbance spectrum and injects a charge into a helical conjugated π-system without injecting a spin. Structural consequences of deprotonation were revealed via analysis of a monoanionic helicene co-crystallized with {K+ (18-crown-6)(THF)} and {Cs+ 2 (18-crown-6)3 }. UV/Vis-monitoring of these systems shows a time-dependent formation of mono- and dianionic species, and the latter was isolated and crystallographically characterized. The ability of the twisted helicene frame to delocalize the negative charge was probed as a perturbation of aromaticity using NICS scans. Relief of strain, avoidance of antiaromaticity, and increase in charge delocalization assist in the additional dehydrogenative ring closures that yield a new planarized decacyclic dianion.

An Additivity Scheme for Aromaticity: The Heteroatom Case.

The nucleus-independent chemical shift (NICS)-XY-Scan is a simple and easy tool for the quantitative measurement of the aromaticity of polycyclic aromatic hydrocarbons and identification of the existence of local and global ring currents. We recently introduced an additivity scheme that uses the NICS-XY-Scans of smaller building blocks to predict the aromatic profiles of larger polycyclic aromatic hydrocarbon systems. We now report on an expansion of the methodology to include systems of varying aromatic natures containing the heteroatoms B, N, O, and S. The additivity approach allows for rapid and resource-efficient generation of NICS-XY-Scans of large, complex systems. Moreover, it reveals that the magnetic criterion of aromaticity behaves in an additive manner, and that the ring currents of multi-ring systems appear to be mostly localized within subunits of up to three rings.

Piecing it Together: An Additivity Scheme for Aromaticity using NICS-XY Scans.

Aromatic compounds are prevalent in both nature and man-made materials, yet their properties are still not fully understood and are therefore hard to predict. Herein, we introduce an additivity scheme for the prediction of the aromatic character of polycyclic aromatic hydrocarbons. Using a small set of building blocks and combination rules, we demonstrate the simple and intuitive construction of complete NICS-XY-scans for several test cases of one- and two-dimensional systems comprising six-membered rings. Partitioning the contribution of discrete building blocks provides insight into the aromatic character of these systems. The results obtained with this methodology provide a new perspective on the distribution of ring currents within polycyclic compounds and the effect of topology on the overall aromatic profile. The concept and strategy presented here are general and highly customizable. The scheme is easily applied to a wide range of interesting systems; it is especially beneficial for the investigation of large systems, as NICS-XY-scans are useful for their predictive utility with respect to optoelectronic and thermochemical properties.

The predictive power of aromaticity: quantitative correlation between aromaticity and ionization potentials and HOMO-LUMO gaps in oligomers of benzene, pyrrole, furan, and thiophene.

Aromaticity is a central and ubiquitous concept in organic chemistry, and is used extensively to explain various phenomena. Yet, aromaticity cannot be observed or measured as a property in its own right and, to date, only qualitative and semi-quantitative relationships have been described between aromaticity and an observable property. We now demonstrate for the first time a robust quantitative relationship between the HOMO-LUMO gap and adiabatic ionization potential of a polycyclic aromatic hydrocarbon oligomer - both measurable physical quantities - and its aromaticity, as quantified by the Nucleus Independent Chemical Shift (NICS) index. The agreement found for a range of structurally and electronically diverse oligomeric systems of varying lengths is so well-behaved as to enable accurate prediction of the properties of longer members of the respective oligomer family. The established correlation allows for preliminary screening of compounds geared towards functional use.

The Carbon-Nitrogen Bonds in Ammonium Compounds Are Charge Shift Bonds.

A comprehensive investigation, utilizing valence-bond- and molecular-orbital-based techniques, reveals that C-N bonds in protonated and methylated ammonium compounds are charge shift bonds. Moreover, the ammonium compounds are predominantly covalent at equilibrium distance, yet have two competing dissociation channels, which determine whether or not a state function crossing occurs during dissociation. The location of the crossing point can be predicted with relative ionization potentials, and may be tuned by altering the stabilization of the respective fragments. By closely examining the nature of the C-N bond as it is stretched, a correlation is observed between the extent of the charge shift nature and bond length. Identification of charge shift bonds in these ubiquitous organic species specifically affects understanding of their role in organocatalysis, in which the C-N bond is stretched.

Response to "Covalent Bonding and Charge Shift Bonds: Comment on 'The Carbon-Nitrogen Bonds in Ammonium Compounds Are Charge Shift Bonds'".

A response to the comment by Gernot Frenking, outlining common ground, as well as differences, with regard to a recent paper on charge-shift bonding in quaternary ammonium cations.

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.