Triply Bonded Pancake π-Dimers Stabilized by Tetravalent Actinides

Aromatic π-stacking is a weakly attractive, noncovalent interaction often found in biological macromolecules and synthetic supramolecular chemistry. The weak nondirectional nature of π-stacking can present challenges in the design of materials owing to their weak, nondirectional nature. However, when aromatic π-systems contain an unpaired electron, stronger attraction involving face-to-face π-orbital overlap is possible, resulting in covalent so-called "pancake" bonds. Two-electron, multicenter single pancake bonds are well known, whereas four-electron double pancake bonds are rare. Higher-order pancake bonds have been predicted, but experimental systems are unknown. Here, we show that six-electron triple pancake bonds can be synthesized by a 3-fold reduction of hexaazatrinaphthylene (HAN) and subsequent stacking of the [HAN]3- triradicals. Our analysis reveals a multicenter covalent triple pancake bond consisting of a σ-orbital and two equivalent π-orbitals. An electrostatic stabilizing role is established for the tetravalent thorium and uranium ions in these systems. We also show that the electronic absorption spectrum of the triple pancake bonds closely matches computational predictions, providing experimental verification of these unique interactions. The discovery of conductivity in thin films of triply bonded π-dimers presents new opportunities for the discovery of single-component molecular conductors and other spin-based molecular materials.

Reimagining the Wave Function in Density Functional Theory: Exploring Strongly Correlated States in Pancake-Bonded Radical Dimers

Pancake bonding between π-conjugated radicals challenges conventional electronic structure approximations, due to the presence of both dispersion (van der Waals) interactions and "strong" electron correlation. Here we use a reimagined wave function-in-density functional theory (DFT) approach to model pancake bonds. Our generalized self-interaction correction extends DFT's reference system of noninteracting electrons, by introducing electron-electron interactions within an active space. We show that a small variation on our previous derivation recovers a DFT-corrected complete active space method proposed by Pijeau and Hohenstein. Comparison of the two approaches shows that the latter provides reasonable dissociation curves for single bonds and pancake bonds, including excited states inaccessible to conventional linear response time-dependent DFT. The results motivate broader adoption of wavefunction-in-DFT approaches for modeling pancake bonds.

Doublet Open-Shell Graphene Fragments

The recent advances on neutral delocalized radical species based on polycyclic aromatic hydrocarbons with fused hexagonal rings, herein defined as doublet open-shell graphene fragments, are summarized in this review. A few simple yet useful theoretical approaches for structural analysis and molecular design were introduced at first. Then, based on the number of fused hexagonal rings, molecular systems with different size, symmetry and edge structure were discussed with emphasis on those isolated in the crystalline form. Their unique self-association behavior, chemical reactivity and physical properties were summarized and discussed, and insights on their functions and potential applications were provided.

Bonding and uneven charge distribution in infinite pyrene π-stacks

Unusual intermolecular π-stacking in a new charge transfer salt of pyrene (Py), (Py)2+(Ga2Cl7)−, has been observed.

Rational design of organic semiconductors with low internal reorganization energies for hole and electron transport: position effect of aza-substitution in phenalenyl derivatives

Amphoteric-redox phenalenyl radical (PLY) is a suitable candidate used to design ambipolar organic materials. Because the singly occupied nonbonding molecular orbital (NBMO) of PLY has a perfect local nonbonding character, its internal reorganization energy (λ) for transporting holes (λ⁺) or electrons (λ^-) is known to be small. Herein, PLY is employed to study the position effect of the aza group on the λ. By adding or extracting an electron from the NBMO, the bond length alterations can be minute. Therefore, the PLY derivatives are also an excellent candidate to study the contributions from the bond angle alterations to the λ. Substituting the aza groups at the β- or α-positions of PLY shows two different trends. When consecutively substituting the aza group at the three β-positions of PLY, the λs are consistently decreased. Contrarily, a series of double functionalization of aza groups at the four α-positions of PLY, the λs are increased. It is because the local bonding or antibonding character in frontier orbitals (FMO) is observed in α2N-PLY and α4N-PLY. As the FMOs of the three β-substituted PLYs and α6N-PLY have perfect local nonbonding character, we found the bond angle alterations are the main contributors of λ. The λs for most aza-PLYs were smaller than 100 meV. Thus, we propose a design rule for substituting aza groups on the parent molecules with strong local nonbonding character in their FMOs. Based on the adiabatic ionization potential and electron affinity, two π-extended PLY derivatives with small λ were recommended for fabricating air-stable ambipolar OFET.

Acene-Linked Zethrenes and Bisphenalenyls: A DFT Search for Organic Tetraradicals

Polycyclic aromatic hydrocarbons are of special interest due to their promising nonlinear optical and magnetic properties. A series of acene-linked zethrenes and bisphenalenyls comprising from five to nine benzene rings in the linker group have been computationally studied by the DFT UB3LYP/6-311++G(d,p) quantum-chemical modeling of their electronic structure, possible spin states, and exchange interactions. The zethrenes with octacene and nonacene linkers as well as bisphenalenyls comprising heptacene, octacene, and nonacene linker groups have been revealed to possess tetraradicaloid nature, which makes them promising building blocks for organic optoelectronic and spintronic devices. The results obtained open a way of constructing tetraradicaloid organic molecules characterized by the presence of two types of paramagnetic centers.

Charge transport properties of open-shell graphene fragments: a computational study of the phenalenyl tilings

Thinking outside the box of the phenalenyl radical: a systematic structure design strategy, phenalenyl tiling, is found to benefit the electron transport properties of open-shell graphene fragments with one free radical. Compared with the closed-shell species, phenalenyl-based π-radicals exhibit smaller intramolecular reorganization energies and larger intermolecular electronic couplings. However, the on-site Coulomb repulsion can be too strong and impedes the charge transport efficiency of such materials. The repulsion can be weakened in radical species by spin delocalization. In this paper, the extended π-radicals we studied are categorized into three types of open-shell structures: the zigzag, the armchair and the discotic odd alternant hydrocarbons. The latter two belong to phenalenyl tilings. We found that the phenalenyl tilings fully inherit the desirable features of the singly occupied molecular orbital of the phenalenyl radical in a predictable and delocalized fashion, and their on-site Coulomb repulsion is effectively reduced. The zigzag π-radicals are less satisfactory. Therefore, the phenalenyl tilings are favorable candidates for charge transporting materials.

Metallic bands in chevron-type polyacenes

We present electronic structure calculations based on a single-parameter plane wave expansion method for basic graphene building blocks, namely n-oligophenylenes and n-oligoacenes, revealing excellent agreement with density-functional theory. When oligophenylene molecules are joined through meta (zigzag) or ortho (chevron) junctions, the resulting molecular dimers and polymers exhibit a semiconducting character. While zigzag dimers of oligoacenes also exhibit gapped electronic structures, their chevron-phase features a sharp metallic band at the Fermi energy. This zero-point-energy state, which transforms into Dirac-like band in chevron polymers, survives at the outer elbows of the dimer irrespective of the molecular length, and has the same origin as reported for the polyacetylene and topologically induced edge states at edge-decorated graphene nanoribbons. These findings assist the engineering of topological electronic states at the molecular level and complement the toolbox of quantum phases in carbon-based nanostructures.

We present electronic structure calculations based on a single-parameter plane wave expansion method for molecular nanostructures revealing excellent agreement with density functional theory and predicting metallic bands for chevron molecular dimers.

The Importance of Electronic Dimensionality in Multiorbital Radical Conductors

The exceptional performance of oxobenzene-bridged bis-1,2,3-dithiazolyls 6 as single-component neutral radical conductors arises from the presence of a low-lying π-lowest unoccupied molecular orbital, which reduces the potential barrier to charge transport and increases the kinetic stabilization energy of the metallic state. As part of ongoing efforts to modify the solid-state structures and transport properties of these so-called multiorbital materials, we report the preparation and characterization of the acetoxy, methoxy, and thiomethyl derivatives 6 (R = OAc, OMe, SMe). The crystal structures are based on ribbonlike arrays of radicals laced together by S···N' and S···O' secondary bonding interactions. The steric and electronic effects of the exocyclic ligands varies, affording one-dimensional (1D) π-stacked radicals for R = OAc, 1D cofacial dimer π-stacks for R = SMe, and a pseudo two-dimensional (2D) brick-wall arrangement for R = OMe. Variable-temperature magnetic and conductivity measurements reveal strong antiferromagnetic interactions and Mott insulating behavior for the two radical-based structures (R = OAc, OMe), with lower room-temperature conductivities (σ_RT ≈ 1 × 10^-4 and ∼1 × 10^-3 S cm^-1, respectively) and higher thermal activation energies ( E_act = 0.24 and 0.21 eV, respectively) than found for the ideal 2D brick-wall structure of 6 (R = F), where σ_RT ≈ 1 × 10^-2 S cm^-1 and E_act = 0.10 eV. The performance of R = OMe, OAc relative to that of R = F, is consistent with the results of density functional theory band electronic structure calculations, which indicate a lower kinetic stabilization energy of the putative metallic state arising from their reduced electronic dimensionality.

Triangular graphene nanofragments: open-shell character and doping

In this work we study the intricacies of the electronic structure properties of triangular graphene nanofragments (TGNFs) in their ground and low-lying excited states by means ofab initioquantum chemistry calculations.

Influence of organic cations on the stacking of semiquinone radical anions

A series of salts of tetrachloro- and tetrabromosemiquinone radical anions reveal four types of stacks: 1) pancake bonded dimers, 2) pancake-bonded trimers, 3) equidistant radicals and 4) a novel type of equidistant stacks of partially charged radicals.

Experimental and Computational Evidence for "Double Pancake Bonds": The Role of Dispersion-Corrected DFT Methods in Strongly Dimerized 5-Aryl-1λ2,3λ2-dithia-2,4,6-triazines

Crystal structures are reported for bicyclic 3-CF₃C₆H₄CN₅S₃ and monocyclic 3-CF₃C₆H₄CN₃S₂, the latter of which is strongly dimerized in a cis-cofacial geometry [3-CF₃C₆H₄CN₃S₂]₂. The title compounds have previously been characterized in solution by NMR, displaying spectra that are consistent with the structure of [3-CF₃C₆H₄CN₃S₂]₂ in the crystal with anti-oriented CF₃ substituents. The interannular binding was investigated using density functional theory (DFT) methods. However, the DFT-optimized geometry spreads the aryl rings too far apart (centroid-centroid distances of ≥4.353 Å versus experimental distance of 3.850 Å). Significant improvements are obtained with dispersion-corrected DFT functionals B3LYP-D3, B3LYP-D3BJ, M062X, and APFD using the 6-311+G(2d,p) basis set. However, all of these overbind the aryl rings with centroid-centroid distances of 3.612, 3.570, 3.526, and 3.511 Å, respectively. After selecting B3LYP-D3BJ/6-311+G(2d,p) as the best method, five alternative dimer geometries were tested, and all were found to be binding; however, anti cofacial-4 (matching the structure in the solid state) is the most stable. Computed energies of the remainder are as follows: +7.0 kJ mol^-1 (syn-cofacial-5), +26.7 kJ mol^-1 (anti-cofacial-64), +27.0 kJ mol^-1 (syn-cofacial-150), +102.0 kJ mol^-1 (S,S-antarafacial), and +103.7 kJ mol^-1 (S,N-antarafacial), where the suffixes are torsional angles around the CN₃S₂ thiazyl ring centroids. The binding in the four most stable cofacial dimers may be described by "double pancake bonding".

Energetics and Electronic Structure of Triangular Hexagonal Boron Nitride Nanoflakes

We studied the energetics and electronic structures of hexagonal boron nitrogen (h-BN) nanoflakes with hydrogenated edges and triangular shapes with respect to the edge atom species. Our calculations clarified that the hydrogenated h-BN nanoflakes with a triangular shape prefer the N edges rather than B edges irrespective of the flake size. The electronic structure of hydrogenated h-BN nanoflakes depends on the edge atom species and their flake size. The energy gap between the lowest unoccupied (LU) and the highest occupied (HO) states of the nanoflakes with N edges is narrower than that of the nanoflakes with B edges and the band gap of h-BN. The nanoflakes possess peculiar non-bonding states around their HO and LU states for the N and B edges, respectively, which cause spin polarization under hole or electron doping, depending on the edge atom species.

Interactions of Metal-Based and Ligand-Based Electronic Spins in Neutral Tripyrrindione π Dimers

The ability of tetrapyrrolic macrocycles to stabilize unpaired electrons and engage in π-π interactions is essential for many electron-transfer processes in biology and materials engineering. Herein, we demonstrate that the formation of π dimers is recapitulated in complexes of a linear tripyrrolic analogue of naturally occurring pigments derived from heme decomposition. Hexaethyltripyrrindione (H₃TD1) coordinates divalent transition metals (i.e., Pd, Cu, Ni) as a stable dianionic radical and was recently described as a robust redox-active ligand. The resulting planar complexes, which feature a delocalized ligand-based electronic spin, are stable at room temperature in air and support ligand-based one-electron processes. We detail the dimerization of neutral tripyrrindione complexes in solution through electron paramagnetic resonance (EPR) and visible absorption spectroscopic methods. Variable-temperature measurements using both EPR and absorption techniques allowed determination of the thermodynamic parameters of π dimerization, which resemble those previously reported for porphyrin radical cations. The inferred electronic structure, featuring coupling of ligand-based electronic spins in the π dimers, is supported by density functional theory (DFT) calculations.

Conformational behavior and stacking interactions of contorted polycyclic aromatics

Computational studies of non-covalent dimers of saddle-shaped molecules unveil widely varying conformations and stacking configurations.

Validation of density functionals for pancake-bonded π-dimers; dispersion is not enough

π-Stacking pancake bonding between radicals poses special challenges to density functional theories (DFTs) due to their shorter than van der Waals contact distances, their multireference singlet ground states and the concurrently important dispersion interactions.

The nature of the multicenter bonding in π-[TCNE]22− dimer: 4c/2e, 12c/2e, or 20c/2e?

Composition of the 20c–2e bonding orbital in the π-[TCNE]₂²⁻ dimer, and the partial occupancy numbers C1, C2 and N in the 20c–2e bond.

Stacked homodimers of substituted contorted hexabenzocoronenes and their complexes with C60 fullerene

Computations show that the tendency of contorted hexabenzocoronene (c-HBC) to form either homodimers or complexes with C₆₀ can be tuned by changing the curvature of the c-HBC via the addition of substituents.

Polyradical Character of Triangular Non-Kekulé Structures, Zethrenes, p-Quinodimethane-Linked Bisphenalenyl, and the Clar Goblet in Comparison: An Extended Multireference Study

In this work, two different classes of polyaromatic hydrocarbon (PAH) systems have been investigated in order to characterize the amount of polyradical character and to localize the specific regions of chemical reactivity: (a) the non-Kekulé triangular structures phenalenyl, triangulene and a π-extended triangulene system with high-spin ground state and (b) PAHs based on zethrenes, p-quinodimethane-linked bisphenalenyl, and the Clar goblet containing varying polyradical character in their singlet ground state. The first class of structures already have open-shell character because of their high-spin ground state, which follows from the bonding pattern, whereas for the second class the open-shell character is generated either because of the competition between the closed-shell quinoid Kekulé and the open-shell singlet biradical resonance structures or the topology of the π-electron arrangement of the non-Kekulé form. High-level ab initio calculations based on multireference theory have been carried out to compute singlet–triplet splitting for the above-listed compounds and to provide insight into their chemical reactivity based on the polyradical character by means of unpaired densities. Unrestricted density functional theory and Hartree–Fock calculations have been performed for comparison also in order to obtain better insight into their applicability to these types of complicated radical systems.

π–π stacking between polyaromatic hydrocarbon sheets beyond dispersion interactions

A slipped parallel structure of a stacked graphene flake showing a biconcave curvature.