Revisiting Kekulene: Synthesis and Single-Molecule Imaging

Collective Magnetism of Spin Coronoid via On-Surface Synthesis

Polyradicals obtained from open-shell coronoids hold promise for applications in spintronics and quantum technologies due to the strong interactions between spins in fully fused cyclic systems. Coronoid synthesis has long been considered difficult due to the cyclization of nanographene. It becomes an immense challenge to synthesize open-shell coronoids since radicals appear only when the macrocycle size exceeds a critical value. Here we present an open-shell coronoid with six radicals achieved through an on-surface synthesis. This spin coronoid displays a collective spin state arising from both the nearest-neighbor exchange interaction and the next-nearest-neighbor exchange interaction of six unpaired π electrons along the conjugation pathways. The characterization of the spin excitation from the ground state to the excited state was carried out by using inelastic electron tunneling spectroscopy. Additionally, we show that the spin coronoid can be utilized as a nanoscale platform to achieve short antiferromagnetic spin-1/2 Heisenberg chains through tip manipulation. Our findings present a design strategy for creating coronoids with polyradicals, which could provide inspiration for fabrication of open-shell coronoid or cyclic spintronic systems.

Helical Quintulene: Synthesis, Chirality, and Supramolecular Assembly

Quintulene is a quintuply symmetrical cycloarene with a positively curved molecular geometry. First described by Staab and Sauer in 1984, its successful synthesis was not achieved until 2020. Due to the challenges posed by its positive curvature, structural extensions of quintulene have been studied rarely. In this work, we report the synthesis of a chiral quintulene (1) featuring five [5]helicene units at its rim constructed through the introduction of steric hindrance. Single-crystal X-ray diffraction analysis confirmed the structure of 1 that showcases a unique chirality combining helicity and rotational symmetry. The enantiomers of 1 were separated and their chiroptical properties were examined, revealing the dissymmetric absorption and luminescence. Additionally, the bowl shape facilitates supramolecular assembly between 1 and fullerenes, with a binding constant to C60 10-fold higher than to C70. This work expands the structural diversity of quintulene and bridges the structures of cycloarenes and multi-helicenes, paving the way for the structural design of chiral cycloarenes.

Two Different Clar's Sextets: Clar's Rule Does Not Necessarily Contradict the Global Pathway of Conjugated Macrocycles

Two contradictory theories, Clar's sextet and resonance theory, explain the stability of conjugated macrocycles using localized and delocalized models, respectively. To reconcile these theories and gain better insights, PAH-containing porphyrinoids were chosen as a representative class of conjugated macrocycles. Two types of Clar's sextets were identified and proposed for the first time based on their interaction with global conjugated pathways: (i) shared sextets that integrate with and perturb the global resonance pathway, and (ii) independent sextets that are separate from the pathway and potentially stabilize or do not perturb it. To test our hypothesis, we synthesized precise regioisomers with and without an independent sextet to exclude variables such as steric hindrance, different PAH characteristics, and other elements. This task is challenging due to the extensive synthetic efforts needed for regioisomeric PAH-containing porphyrinoids. We focused on 2,3-vinylnaphthiporphyrin, which features shared sextets, because its o-regioisomers are readily accessible. Two distinct regioisomers (1,2-vinylnaphthiporphyrins) with an independent sextet were synthesized along with unexpected nonaromatic N-fused porphyrinoids. Our analysis shows that 1,2-vinylnaphthiporphyrins with an independent sextet are more aromatic than 2,3-vinylnaphthiporphyrin (nonaromatic). This sextet analysis was also applied to other reported and yet-to-be-synthesized PAH-containing porphyrinoids, and these results are consistent with our current study.

Recent Progress of Imaging Chemical Bonds by Scanning Probe Microscopy: A Review

In the past decades, the invention of scanning probe microscopy (SPM) as the versatile surface-based characterization of organic molecules has triggered significant interest throughout multidisciplinary fields. In particular, the bond-resolved imaging acquired by SPM techniques has extended its fundamental function of not only unraveling the chemical structure but also allowing us to resolve the structure-property relationship. Here, we present a systematical review on the history of chemical bonds imaged by means of noncontact atomic force microscopy (nc-AFM) and bond-resolved scanning tunneling microscopy (BR-STM) techniques. We first summarize the advancement of real-space imaging of covalent bonds and the investigation of intermolecular noncovalent bonds. Beyond the bond imaging, we also highlight the applications of the bond-resolved SPM techniques such as on-surface synthesis, the determination of the reaction pathway, the identification of molecular configurations and unknown products, and the generation of artificial molecules created via tip manipulation. Lastly, we discuss the current status of SPM techniques and highlight several key technical challenges that must be solved in the coming years. In comparison to the existing reviews, this work invokes researchers from surface science, chemistry, condensed matter physics, and theoretical physics to uncover the bond-resolved SPM technique as an emerging tool in exploiting the molecule/surface system and their future applications.

In-Depth Theoretical Investigations of Borazine's Aromaticity: Tailoring Electron Delocalization through Substituent Effects

The current study investigates the influence of several R substituents (e.g., Me, SiH3, F, Cl, Br, OH, NH2, etc.) on the aromaticity of borazine, also known as the "inorganic benzene". By performing hybrid DFT methods, blended with several computational techniques, e.g., Natural Bond Orbital (NBO), Quantum Theory of Atoms in Molecules (QTAIM), Gauge-Including Magnetically Induced Current (GIMIC), Nucleus-Independent Chemical Shift (NICS), and following a simultaneous evaluation of four different aromaticity indices (para-delocalization index (PDI), multi-centre bond order (MCBO), ring current strength (RCS), and NICS parameters), it is emphasized that the aromatic character of B-substituted (B3R3N3H3) and N-substituted (B3H3N3R3) borazine derivatives can be tailored by modulating the electronic effects of R groups. It is also highlighted that the position of R substituents on the ring structure is crucial in tuning the aromaticity. Systematic comparisons of calculated aromaticity index values (i.e., via regression analyses and correlation matrices) ensure that the reported trends in aromaticity variation are accurately described, while the influence of different R groups on electron delocalization and related aromaticity phenomena is quantitatively assessed based on NBO analyses. The most relevant interactions impacting the aromatic character of investigated systems are (i) the electron conjugations occurring between the p lone pair electrons (LP) on the F, Cl, Br, O or N atoms, of R groups, and the π*(B=N) orbitals on the borazine ring (i.e., LP(R)→π*(B=N) donations), and (ii) the steric-exchange (Pauli) interactions between the same LP and the π(B=N) bonds (i.e., LP(R)↔π(B=N) repulsions), while inductive/field effects influence the aromaticity of the investigated trisubstituted borazine systems to a much lesser extent. This work highlights that although the aromatic character of borazine can be enhanced by grafting electron-donor substituents (F, OH, NH2, O-, NH-) on the N atoms, the stabilization due to aromaticity has only a moderate impact on these systems. By replacing the H substituents on the B atoms with similar R groups, the aromatic character of borazine is decreased due to strong exocyclic LP(R)→π*(B=N) donations affecting the delocalization of π-electrons on the borazine ring.

Generalized kekulenes and clarenes as novel families of cycloarenes: structures, stability, and spectroscopic properties

Cycloarenes constitute a captivating class of polycyclic aromatic hydrocarbons with unique structures and properties, but their synthesis represents a challenging task in organic chemistry. Kekulenes and edge-extended kekulenes as classic types of cycloarenes play an important role in the comprehension of π electron distribution, but their sparse molecular diversity considerably limits their further development and application. In this work, we propose two novel classes of cycloarenes, the generalized kekulenes and the clarenes. Using density functional theory, we carry out a comprehensive study of all possible isomers of the generalized kekulenes and clarenes with different sizes. By applying a simple Hückel model, we show that π delocalization plays a crucial role in determining the relative stability of isomers. We also discover that π-π stacking is commonly present in certain larger clarenes and provides a considerable additional stabilization effect, making the corresponding isomers the lowest-energy ones. Among all considered typical looped polyarenes, generalized kekulenes and/or clarenes are revealed to be the energetically most stable forms, suggesting that these novel cycloarenes proposed here would be viable targets for future synthetic work. The simulated 1H NMR spectra and UV-vis absorption spectra provide valuable information about the electronic and optoelectronic properties for the most stable generalized kekulene and clarene species and may support their identification in future synthesis and experimental characterization.

Butterfly-Shaped Dibenz[a,j]anthracenes: Synthesis and Photophysical Properties

A strategy for the synthesis of dibenz[a,j]anthracenes (DBAs) from cyclohexa-2,5-diene-1-carboxylic acids is presented. Our approach involves sequential C-H olefination, cycloaddition, and decarboxylative aromatization. In the key step for DBA skeleton construction, the bis-C-H olefination products, 1,3-dienes, are utilized as substrates for [4 + 2] cycloaddition with benzyne. This concise synthetic route allows for regioselective ring formation and functional group introduction. The structural features and photophysical properties of the resulting DBA molecules are discussed.

Electronic properties and carrier mobilities of nanocarbons formed by non-benzoidal building blocks

Exotic 1D and 2D carbon nanostructures have been grown in the laboratory in the last few years by means of surface-assisted chemical routes. In these processes, the strategical choice of a molecular precursor plays a dominant role in the determination of the synthesized nanocarbon. Further variations of these techniques are able to produce non-benzoidal carbon quantum-dots (QDs). Considering this experimental scenario as motivation, we propose a series of nanoribbon systems based on concatenating recently synthesized carbon QDs containing pentagonal, hexagonal, and heptagonal rings. We use density functional theory (DFT) simulations to reveal their properties can range from metallic to semiconducting depending on the concatenation hierarchy used to form the nanoribbons. This DFT implementation is based on a LCAO approach to describe valence wavefunctions and most of the simulations employ the PBE-GGA functional. Since this functional is known to underestimate band gaps, we also use the B3LYP functional in a plane-wave DFT approach for a selected case for comparison purposes. These systems show a different gap versus width relationship compared to conventional graphene nanoribbons setups and a particular set of carrier mobility values. We further discuss the interplay between the QD's frontier states and the electronic properties of the nanoribbons in light of their structural details.

Computational Analysis of Some More Rectangular Tessellations of Kekulenes and Their Molecular Characterizations

Cycloarene molecules are benzene-ring-based polycyclic aromatic hydrocarbons that have been fused in a circular manner and are surrounded by carbon-hydrogen bonds that point inward. Due to their magnetic, geometric, and electronic characteristics and superaromaticity, these polycyclic aromatics have received attention in a number of studies. The kekulene molecule is a cyclically organized benzene ring in the shape of a doughnut and is the very first example of such a conjugated macrocyclic compound. Due to its structural characteristics and molecular characterizations, it serves as a great model for theoretical research involving the investigation of π electron conjugation circuits. Therefore, in order to unravel their novel electrical and molecular characteristics and foresee potential applications, the characterization of such components is crucial. In our current research, we describe two unique series of enormous polycyclic molecules made from the extensively studied base kekulene molecule, utilizing the essential graph-theoretical tools to identify their structural characterization via topological quantities. Rectangular kekulene Type-I and rectangular kekulene Type-II structures were obtained from base kekulene molecules arranged in a rectangular fashion. We also employ two subcases for each Type and, for all of these, we derived ten topological indices. We can investigate the physiochemical characteristics of rectangular kekulenes using these topological indices.

Bottom-up Solution Synthesis of Graphene Nanoribbons with Precisely Engineered Nanopores

The incorporation of nanopores into graphene nanostructures has been demonstrated as an efficient tool in tuning their band gaps and electronic structures. However, precisely embedding the uniform nanopores into graphene nanoribbons (GNRs) at the atomic level remains underdeveloped especially for in-solution synthesis due to the lack of efficient synthetic strategies. Herein we report the first case of solution-synthesized porous GNR (pGNR) with a fully conjugated backbone via the efficient Scholl reaction of tailor-made polyphenylene precursor (P1) bearing pre-installed hexagonal nanopores. The resultant pGNR features periodic subnanometer pores with a uniform diameter of 0.6 nm and an adjacent-pores-distance of 1.7 nm. To solidify our design strategy, two porous model compounds (1 a, 1 b) containing the same pore size as the shortcuts of pGNR, are successfully synthesized. The chemical structure and photophysical properties of pGNR are investigated by various spectroscopic analyses. Notably, the embedded periodic nanopores largely reduce the π-conjugation degree and alleviate the inter-ribbon π-π interactions, compared to the nonporous GNRs with similar widths, affording pGNR with a notably enlarged band gap and enhanced liquid-phase processability.

The Topology of Molecules with Twelve Fused Phenyl Rings ([12]Circulenes): Rings, Infinitenes, and Möbius Infinitenes

Following the recent preparation of infinitene (J. Am. Chem. Soc. 2022, 144, 862-871), a computational (ωB97XD/6-311G(d)) exploration of 42 isomeric compounds with 12 fused phenyl rings identified structures with linking number of zero (ring, saddle, and ribbon shapes), two (infinitene-like shape), and one (Möbius infinitene shape) is reported. An infinitene isomer composed of two [5]helicene fragments connected to two stacked phenyl rings and a Möbius infinitene isomer are identified that are more stable than the known infinitene. The energies of the structures are examined by assessing their macrocyclization (strain) energies, π-stacking, and possible aromaticity. Examples of fused phenyl molecules with linking numbers of 3, 4, 5, and 6 are shown, indicating the potential topological range that these molecules can possess.

Infinitenes as the Most Stable Form of Cycloarenes: The Interplay among π Delocalization, Strain, and π-π Stacking

The recent successful preparation of infinitene has sparked widespread attention due to its aesthetic appeal and synthetic challenge. Spectroscopic measurements and follow-up computational investigations suggest that infinitene holds fundamental significance and potential applications in chiroptics, optoelectronics, asymmetric synthesis, and supramolecular chemistry. However, unlike other looped polyarenes enriched with sizes and shapes, the infinitene molecule seems, so far, the only known example of this fascinating new form of nanocarbons, whose further exploitation would be considerably limited because of the lack of molecular diversity. Here, we introduce a whole new family of generalized infinitenes with different sizes and topologies. Three types of infinitene structures are rationally designed by joining two units of coronene, kekulene, or their extended analogs. The constructed molecules of varying sizes, each with a large number of possible topoisomers, are systematically studied by DFT calculations. Comprehensive analysis using a simple energy decomposition model uncovers that the stability of infinitenes is governed by the interplay among π delocalization, steric strain, and π-π stacking. While the first two factors are crucial to the stability of smaller infinitenes, the latter is the primary stabilizing interaction for larger infinitenes. Most importantly, we show that larger-sized infinitenes are actually the energetically most favorable form among all known looped polyarenes; their substantial thermodynamic stability surpassing that of circulenes, various carbon nanobelts, and kekulene-like macrocycles renders them promising targets for synthesis. The simulated ¹H NMR, UV-vis, and circular dichroism spectra along with optical rotations for the most stable infinitene species may help their identification in future synthetic efforts.

Length-dependent symmetry in narrow chevron-like graphene nanoribbons

We report the structural and electronic properties of narrow chevron-like graphene nanoribbons (GNRs), which depending on their length are either mirror or inversion symmetric. Additionally, GNRs of different length can form molecular heterojunctions based on an unusual binding motif.

On-Surface Synthesis of C144 Hexagonal Coronoid with Zigzag Edges

Coronoids as polycyclic aromatic macrocycles enclosing a cavity have attracted a lot of attention due to their distinctive molecular and electronic structures. They can be also regarded as nanoporous graphene molecules whose electronic properties are critically dependent on the size and topology of their outer and inner peripheries. However, because of their synthetic challenges, the extended hexagonal coronoids with zigzag outer edges have not been reported yet. Here, we report the on-surface synthesis of C144 hexagonal coronoid with outer zigzag edges on a designed precursor undergoing hierarchical Ullmann coupling and cyclodehydrogenation on the Au(111) surface. The molecular structure is unambiguously characterized by bond-resolved noncontact atomic force microscopy imaging. The electronic properties are further investigated by scanning tunneling spectroscopy measurements, in combination with the density functional theory calculations. Moreover, the values of the harmonic oscillator model of aromaticity are derived from calculations that suggest that the molecular structure is ideally represented by Clar's model. Our results provide approaches toward realizing a hexagonal coronoid with zigzag edges, potentially inspiring fabrication of hexagonal zigzag coronoids with multiple radical characters in the future.

Quantitative Resonance Theory Based on the Clar Sextet Model

Despite its great explanatory power in understanding the chemistry of polycyclic aromatic hydrocarbons (PAHs) and related systems, the Clar sextet rule still remains an intuitive and qualitative model with notable exceptions in some cases. Here we develop a quantitative theory of chemical resonance based on semilocalized Clar-type resonance structures (named the Clar resonators) consisting of variable numbers of Clar sextets and C═C bonds. The constructed wave functions of the Clar resonators are used to expand the actual wave function of the π-conjugated system obtained from a DFT or Hartree-Fock calculation. The resultant weights and one-electron energies of the Clar resonators can serve as a quantitative measure of the importance of these resonators. Implementing the theory in our open-source python code EzReson and applying it to over a thousand PAH molecules of different sizes and shapes, we show that the weight of the Clar resonators increases exponentially with increasing number of sextets and that their energy decreases linearly with the latter, thus confirming the general validity of the Clar rule. On the basis of such a large-scale resonance analysis, we propose three extended Clar rules, along with a unified quantitative model, that are able to evaluate the importance of all Clar resonators and the ring aromaticity for PAHs. Using the present theories, we uncover the essential role that the minor Clar resonators may play in correctly understanding the resonance stabilization and local aromaticity of rings, which was totally overlooked in the original Clar model.

Chemisorption-Induced Formation of Biphenylene Dimer on Ag(111)

We report an example that demonstrates the clear interdependence between surface-supported reactions and molecular-adsorption configurations. Two biphenyl-based molecules with two and four bromine substituents, i.e., 2,2'-dibromobiphenyl (DBBP) and 2,2',6,6'-tetrabromo-1,1'-biphenyl (TBBP), show completely different reaction pathways on a Ag(111) surface, leading to the selective formation of dibenzo[e,l]pyrene and biphenylene dimer, respectively. By combining low-temperature scanning tunneling microscopy, synchrotron radiation photoemission spectroscopy, and density functional theory calculations, we unravel the underlying reaction mechanism. After debromination, a biradical biphenyl can be stabilized by surface Ag adatoms, while a four-radical biphenyl undergoes spontaneous intramolecular annulation due to its extreme instability on Ag(111). Such different chemisorption-induced precursor states between DBBP and TBBP consequently lead to different reaction pathways after further annealing. In addition, using bond-resolving scanning tunneling microscopy and scanning tunneling spectroscopy, we determine with atomic precision the bond-length alternation of the biphenylene dimer product, which contains 4-, 6-, and 8-membered rings. The 4-membered ring units turn out to be radialene structures.

Symmetry and Combinatorial Concepts for Cyclopolyarenes, Nanotubes and 2D-Sheets: Enumerations, Isomers, Structures Spectra & Properties

This review article highlights recent developments in symmetry, combinatorics, topology, entropy, chirality, spectroscopy and thermochemistry pertinent to 2D and 1D nanomaterials such as circumscribed-cyclopolyarenes and their heterocyclic analogs, carbon and heteronanotubes and heteronano wires, as well as tessellations of cyclopolyarenes, for example, kekulenes, septulenes and octulenes. We establish that the generalization of Sheehan’s modification of Pólya’s theorem to all irreducible representations of point groups yields robust generating functions for the enumeration of chiral, achiral, position isomers, NMR, multiple quantum NMR and ESR hyperfine patterns. We also show distance, degree and graph entropy based topological measures combined with techniques for distance degree vector sequences, edge and vertex partitions of nanomaterials yield robust and powerful techniques for thermochemistry, bond energies and spectroscopic computations of these species. We have demonstrated the existence of isentropic tessellations of kekulenes which were further studied using combinatorial, topological and spectral techniques. The combinatorial generating functions obtained not only enumerate the chiral and achiral isomers but also aid in the machine construction of various spectroscopic and ESR hyperfine patterns of the nanomaterials that were considered in this review. Combinatorial and topological tools can become an integral part of robust machine learning techniques for rapid computation of the combinatorial library of isomers and their properties of nanomaterials. Future applications to metal organic frameworks and fullerene polymers are pointed out.

Infinitene: A Helically Twisted Figure-Eight [12]Circulene Topoisomer

New forms of molecular nanocarbon particularly looped polyarenes adopting various topologies contribute to the fundamental science and practical applications. Here we report the synthesis of an infinity-shaped polyarene, infinitene (1) (cyclo[c.c.c.c.c.c.e.e.e.e.e.e]dodecakisbenzene), comprising consecutively fused 12-benzene rings forming an enclosed loop with a strain energy of 60.2 kcal·mol^-1. Infinitene (1) represents a topoisomer of still-hypothetical [12]circulene, and its scaffold can be formally visualized as the outcome of the "stitching" of two homochiral [6]helicene subunits by both their ends. The synthetic strategy encompasses transformation of a rationally designed dithiacyclophane to cyclophadiene through the Stevens rearrangement and pyrolysis of the corresponding S,S'-bis(oxide) followed by the photocyclization. The structure of 1 is a unique hybrid of helicene and circulene with a molecular formula of C₄₈H₂₄, which can be regarded as an isomer for kekulene, [6,6]carbon nanobelt ([6,6]CNB), and [12]cyclacene. Infinitene (1) is a bench-stable yellow solid with green fluorescence and soluble to common organic solvents. Its figure-eight molecular structure was unambiguously confirmed by X-ray crystallography. The scaffold of 1 is significantly compressed as manifested by a remarkably shortened distance (3.152-3.192 Å) between the centroids of two π-π stacked central benzene rings and the closest C···C distance of 2.920 Å. Fundamental photophysical properties of 1 were thoroughly elucidated by UV-vis absorption and fluorescence spectroscopic studies and density functional theory calculations. Its configurational stability enabled separation of the corresponding enantiomers (P,P) and (M,M) by a chiral HPLC. Circular dichroism (CD) and circularly polarized luminescence (CPL) measurements revealed that 1 has moderate |g_CD| and |g_CPL| values.

Dianion and Dication of Tetracyclopentatetraphenylene as Decoupled Annulene-within-an-Annulene Models

The dianion and dication of tetramesityl-substituted tetracyclopentatetraphenylene, a circulene consisting of alternating five- and six-membered rings, have been generated by reduction with alkali metals and oxidation with antimony(V) halides, respectively. They are theoretically predicted to adopt double annulenoid structures called annulene-within-an-annulene models in which the outer and inner conjugation circuits are significantly decoupled. The theoretical structures were experimentally proven by X-ray crystallographic analyses and the electronic configurations were supported by MCD spectra. Based on the ¹³ C NMR chemical shifts, negative and positive charges are shown to be located mainly at the outer periphery, indicating that the dianion and dication have delocalized 22-π and 18-π electron outer perimeters, respectively, and 8-π electron structure at the inner ring. Notably, the dianion has an open-shell character, whereas the dication has a closed-shell ground state.

Synthesis and Chiral Resolution of Twisted Carbon Nanobelts

Twisted carbon nanobelts could display persistent chirality, which is desirable for material applications, but their synthesis is very challenging. Herein, we report the successful synthesis and chiral resolution of such a kind of molecules (1-H and 1) with a figure-eight configuration. 1-H was synthesized first by macrocyclization through Suzuki coupling reaction followed by benzannulation via Bi(OTf)₃-mediated cyclization reaction of vinyl ether. Oxidative dehydrogenation of 1-H gave the fully π-conjugated 1. Their twisted structures were confirmed by X-ray crystallographic analysis and calculations, and they can be resolved by chiral high-performance liquid chromatography. The isolated enantiomers showed persistent chiroptical properties according to the circular dichroism measurements, with moderate |g_abs| values (0.0016 for 1-H and 0.005-0.007 for 1). Their photophysical properties were also briefly studied.

Expanded Kekulenes

The synthesis of kekulene and its higher homologues is a challenging task in organic chemistry. The first successful synthesis and characterization of the parent kekulene were reported by Diederich and Staab in 1978. Herein, we report the facile preparation of a series of edge-extended kekulenes by bismuth(III) triflate-catalyzed cyclization of vinyl ethers from the properly designed macrocyclic precursors. Their molecular structures were confirmed by X-ray crystallographic analysis and NMR spectroscopy. Their size- and symmetry-dependent electronic structures (frontier molecular orbitals, aromaticity) and physical properties (optical and electrochemical) were investigated by various spectroscopic measurements, assisted by theoretical calculations. Particularly, the acene-like units along each zigzag edge demonstrate a dominant local aromatic character. Our studies provide an easy synthetic strategy toward various fully fused carbon nanostructures and give some insights into the electronic properties of cycloarenes.

Cluster many-body expansion: A many-body expansion of the electron correlation energy about a cluster mean field reference

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. In the past, applications of MBE to correlation energy have been unfeasible for large systems, but recent improvements to computing resources have sparked renewed interest in capturing the correlation energy using the generalized nth order Bethe-Goldstone equation. In this work, we extend this approach, originally proposed for a Slater determinant, to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean field reference TPS configuration and includes the correlation contribution of the excited TPSs using the MBE. This method, named cluster MBE (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from a RHF reference. We present numerical results for strongly correlated systems, such as the one- and two-dimensional Hubbard models and the chromium dimer. The performance of the cMBE method is also tested by partitioning the extended π space of several large π-conjugated systems, including a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 10⁶⁶ determinants for the exact FCI solution.

On-Surface Synthesis and Molecular Engineering of Carbon-Based Nanoarchitectures

On-surface synthesis via covalent coupling of adsorbed precursor molecules on metal surfaces has emerged as a promising strategy for the design and fabrication of novel organic nanoarchitectures with unique properties and potential applications in nanoelectronics, optoelectronics, spintronics, catalysis, etc. Surface-chemistry-driven molecular engineering (i.e., bond cleavage, linkage, and rearrangement) by means of thermal activation, light irradiation, and tip manipulation plays critical roles in various on-surface synthetic processes, as exemplified by the work from the Ernst group in a prior issue of ACS Nano. In this Perspective, we highlight recent advances in and discuss the outlook for on-surface syntheses and molecular engineering of carbon-based nanoarchitectures.

The effect of spin polarization on the electron transport of molecular wires with diradical character

Some of the most promising materials for application in molecular electronics and spintronics are based on diradical chains. Herein, the proposed relation between increasing conductance with length and diradical character is revisited using ab initio methods that account for the static electron correlation effects. Electron transmission was previously obtained from restricted single determinant wavefuntions or tight-binding approximations, which are unable to account for static correlation. Broken Symmetry Unrestricted Kohn-Sham Density Functional Theory (BS-UKS-DFT) in combination with electron transport analysis based on electron deformation orbitals (EDOs) reflects an exponential decay of the electrical conductance with length. Also, other important effects such as quantum interference are correctly accounted for, leading to a decrease of the conductance as the diradical character increases. As a proof-of-concept, the electrical conductance obtained from BS-UKS-DFT and CASSCF(2,2) wavefunctions were compared in diradical graphene strips in the frame of the pseudo-π approach, obtaining very similar results.

A simple Hückel model-driven strategy to overcome electronic barriers to retro-Brook silylation relevant to aryne and bisaryne precursor synthesis

ortho-Silylaryl triflate precursors (oSATs) have been responsible for many recent advances in aryne chemistry and are most commonly accessed from the corresponding 2-bromophenol. A retro-Brook O- to C-silyl transfer is a key step in this synthesis but not all aromatic species are amenable to the transformation, with no functionalized bisbenzyne oSATs reported. Simple Hückel models are presented which show that the calculated aromaticity at the brominated position is an accurate predictor of successful retro-Brook reaction, validated synthetically by a new success and a predicted failure. From this, the synthesis of a novel difunctionalized bisaryne precursor has been tested, requiring different approaches to install the two C-silyl groups. The first successful use of a disubstituted o-silylaryl sulfonate bisbenzyne precursor in Diels-Alder reactions is then shown.

On-Surface Synthesis and Characterization of [7]Triangulene Quantum Ring

The ability to engineer geometrically well-defined antidots in large triangulene homologues allows for creating an entire family of triangulene quantum rings (TQRs) with tunable high-spin ground state, crucial for next-generation molecular spintronic devices. Herein, we report the synthesis of an open-shell [7]triangulene quantum ring ([7]TQR) molecule on Au(111) through the surface-assisted cyclodehydrogenation of a rationally designed kekulene derivative. Bond-resolved scanning tunneling microscopy (BR-STM) unambiguously imaged the molecular backbone of a single [7]TQR with a triangular zigzag edge topology, which can be viewed as [7]triangulene decorated with a coronene-like antidot in the center. Additionally, dI/dV mapping reveals that both inner and outer zigzag edges contribute to the edge-localized and spin-polarized electronic states of [7]TQR. Both experimental results and spin-polarized density functional theory calculations indicate that [7]TQR retains its open-shell septuple ground state (S = 3) on Au(111). This work demonstrates a new route for the design of high-spin graphene quantum rings for future quantum devices.

Ultrahigh-yield on-surface synthesis and assembly of circumcoronene into a chiral electronic Kagome-honeycomb lattice

On-surface synthesis has revealed remarkable potential in the fabrication of atomically precise nanographenes. However, surface-assisted synthesis often involves multiple-step cascade reactions with competing pathways, leading to a limited yield of target nanographene products. Here, we devise a strategy for the ultrahigh-yield synthesis of circumcoronene molecules on Cu(111) via surface-assisted intramolecular dehydrogenation of the rationally designed precursor, followed by methyl radical-radical coupling and aromatization. An elegant electrostatic interaction between circumcoronenes and metallic surface drives their self-organization into an extended superlattice, as revealed by bond-resolved scanning probe microscopy measurements. Density functional theory and tight-binding calculations reveal that unique hexagonal zigzag topology of circumcoronenes, along with their periodic electrostatic landscape, confines two-dimensional electron gas in Cu(111) into a chiral electronic Kagome-honeycomb lattice with two emergent electronic flat bands. Our findings open up a new route for the high-yield fabrication of elusive nanographenes with zigzag topologies and their superlattices with possible nontrivial electronic properties.

Phenanthrenes/dihydrophenanthrenes: the selectivity controlled by different benzynes and allenes

A method for the intermolecular annulation of benzynes with allenes is disclosed. This protocol utilized allenes as an unconventional diene component for the selective synthesis of phenanthrenes and dihydrophenanthrenes under the control of different benzyne precursors, featuring high atom-economy and good functional group compatibility. Density functional theory (DFT) calculations reveal that different migratory routes of the aromatic C-H bond are crucial for the observed selectivity.

Clar Goblet and Aromaticity Driven Multiradical Nanographenes

The Clar Goblet, the first radical bowtie nanographene proposed by Erich Clar nearly 50 years ago, was recently synthesized. Bowtie nanographenes present quasi-degenerate magnetic ground states, which make them so elusive as unique. A thorough analysis is presented of the spin-state energetics of Clar Goblet and bowtie nanographenes by a battery of existing and novel ab initio procedures ranging from density functional theory to complete active space Hamiltonians. With this, it was proven that π radicals of bowtie nanographenes sit on BP (Benzo[cd]Pyrene) moieties driven by their local aromaticity, a purely chemical concept, which confers global stability to the whole structure. Besides, a novel Pauli energy densities analysis provided a visual intuitive explanation for this preference. These findings allow envisioning that analogous bowtie nanographenes with arbitrary polyradical character are not only feasible at the molecular scale but will share Clar Goblet's peculiar properties.

Kekulene: On-Surface Synthesis, Orbital Structure, and Aromatic Stabilization

We revisit the question of kekulene's aromaticity by focusing on the electronic structure of its frontier orbitals as determined by angle-resolved photoemission spectroscopy. To this end, we have developed a specially designed precursor, 1,4,7(2,7)-triphenanthrenacyclononaphane-2,5,8-triene, which allows us to prepare sufficient quantities of kekulene of high purity directly on a Cu(111) surface, as confirmed by scanning tunneling microscopy. Supported by density functional calculations, we determine the orbital structure of kekulene's highest occupied molecular orbital by photoemission tomography. In agreement with a recent aromaticity assessment of kekulene based solely on C-C bond lengths, we conclude that the π-conjugation of kekulene is better described by the Clar model rather than a superaromatic model. Thus, by exploiting the capabilities of photoemission tomography, we shed light on the question which consequences aromaticity holds for the frontier electronic structure of a π-conjugated molecule.

Kekulene: On-Surface Synthesis, Orbital Structure, and Aromatic Stabilization

We revisit the question of kekulene's aromaticity by focusing on the electronic structure of its frontier orbitals as determined by angle-resolved photoemission spectroscopy. To this end, we have developed a specially designed precursor, 1,4,7(2,7)-triphenanthrenacyclononaphane-2,5,8-triene, which allows us to prepare sufficient quantities of kekulene of high purity directly on a Cu(111) surface, as confirmed by scanning tunneling microscopy. Supported by density functional calculations, we determine the orbital structure of kekulene's highest occupied molecular orbital by photoemission tomography. In agreement with a recent aromaticity assessment of kekulene based solely on C-C bond lengths, we conclude that the π-conjugation of kekulene is better described by the Clar model rather than a superaromatic model. Thus, by exploiting the capabilities of photoemission tomography, we shed light on the question which consequences aromaticity holds for the frontier electronic structure of a π-conjugated molecule.

Synthesis and assembly of extended quintulene

Quintulene, a non-graphitic cycloarene with fivefold symmetry, has remained synthetically elusive due to its high molecular strain originating from its curved structure. Here we report the construction of extended quintulene, which was unambiguously characterized by mass and NMR spectroscopy. The extended quintulene represents a naturally curved nanocarbon based on its conical molecular geometry. It undergoes dimerization in solution via π-π stacking to form a metastable, but isolable bilayer complex. Thermodynamic and kinetic characterization reveals the dimerization process as entropy-driven and following second-order kinetics with a high activation energy. These findings provide a deeper understanding of the assembly of conical nanocarbons. Comparison of optical properties of monomer and dimer points toward a H-type interlayer coupling in the dimer.

Large-Cavity Coronoids with Different Inner and Outer Edge Structures

Coronoids, polycyclic aromatic hydrocarbons with geometrically defined cavities, are promising model structures of porous graphene. Here, we report the on-surface synthesis of C168 and C140 coronoids, referred to as [6]- and [5]coronoid, respectively, using 5,9-dibromo-14-phenylbenzo[m]tetraphene as the precursor. These coronoids entail large cavities (>1 nm) with inner zigzag edges, distinct from their outer armchair edges. While [6]coronoid is planar, [5]coronoid is not. Low-temperature scanning tunneling microscopy/spectroscopy and noncontact atomic force microscopy unveil structural and electronic properties in accordance with those obtained from density functional theory calculations. Detailed analysis of ring current effects identifies the rings with the highest aromaticity of these coronoids, whose pattern matches their Clar structure. The pores of the obtained coronoids offer intriguing possibilities of further functionalization toward advanced host–guest applications.