Magnetically Induced Ring-Current Strengths of Planar and Nonplanar Molecules: New Insights from the Pseudo-π Model

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.

Aromaticity in Semi-Condensed Figure-Eight Molecules

Electron delocalization and aromaticity was comparatively evaluated in recently synthesized figure-eight molecules made of two condensed U-shaped polycyclic aromatic hydrocarbon moieties connected either by two single bonds or by two para-phenylene groups. The selected examples include molecules that incorporate eight-membered and sixteen-membered rings, as well as a doubly [5]helicene-bridged (1,4)cyclophane. We probe whether some electron delocalization could occur through the stereogenic single bonds in these molecules: Is aromaticity purely (semi-)local, or possibly also global in these molecules? It was concluded that the situation can go from a purely (semi-)local character when the dihedral angle at the connecting single bonds is large, such as in biphenyl, to a predominantly (semi-)local character with a minor global contribution when the dihedral angle is small, such as in the para-phenylene connectors of the [5] helicene-bridged cyclophane.

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.

Understanding Aromaticity in [5]Helicene-Bridged Cyclophanes: A Comprehensive Study

This study explores the aromaticity of doubly [5]helicene-bridged (1,4)cyclophane and triply [5]helicene-bridged (1,3,5)cyclophane via calculations of the magnetic response and of electronic aromaticity indices. The primary objective is to assess the π-electron delocalization to determine whether they sustain global ring currents associated with π aromaticity. The molecules show local ring currents in the presence of an external magnetic field. The ring currents flow diatropically in the stacked six-membered rings and in the helicene arms. However, these π currents are not interconnected due to the discontinuity of the π delocalization at the C-C single bonds connecting the central six-membered rings to the helicene arms. Electronic indices suggest that the helicene-arm systems have significantly smaller electron delocalization than benzene. The reduction in the delocalization does not compromise their ability to exhibit ring currents in the presence of an external magnetic field. The analysis provides further evidence that the magnetic criteria yield a different degree of aromaticity for the helicene arms than obtained in the calculation of the electronic aromaticity indices. However, both approaches confirm that the studied molecules are not globally aromatic.

Isomerism tunes the diradical character of difluorenopyrroles at constant Hückel-level anti-aromaticity

A new diradical based on diindenocarbazole or difluorenopyrrole was synthesized and experimentally characterized by optical, electrochemical, and magnetic techniques, as well as quantum chemical calculations. The isomerism of these structures tunes the diradical character and the associated properties, representing a unique case of such important modulation. A full study of the electronic structure was carried out considering the perturbative interactions between different canonical forms as well as the anti-aromatic character of the molecular cores. Such a study reveals how we can tune diradical character simply by reorganizing the bonding patterns at constant chemical costs (composition).

Current-Density Calculations on Zn-Porphyrin40 Nanorings

Two porphyrinoid nanorings have been studied computationally. They were built by linking 40 Zn-porphyrin units with butadiyne bridges. The molecular structures belonging to the D40h point group were fully optimized with the Turbomole program at the density functional theory (DFT) level using the B3LYP functional and the def2-SVP basis sets. The aromatic character was studied at the DFT level by calculating the magnetically induced current-density (MICD) susceptibility using the GIMIC program. The neutral molecules are globally non-aromatic with aromatic Zn-porphyrin units. Charged nanorings could not be studied because almost degenerate frontier orbitals led to vanishing optical gaps for the cations. Since DFT calculations of the MICD are computationally expensive, we also calculated the MICD using three pseudo-π models. Appropriate pseudo-π models were constructed by removing the outer hydrogen atoms and replacing all carbon and nitrogen atoms with hydrogen atoms. The central Zn atom was either replaced with a beryllium atom or with two inner hydrogen atoms. Calculations with the computationally inexpensive pseudo-π models yielded qualitatively the same magnetic response as obtained in the all-electron calculations.

Charge delocalization and aromaticity of doubly reduced double-walled carbon nanohoops

Cycloparaphenylenes (CPPs) exhibit selective host capabilities, featuring the ability to incorporate smaller CPPs to form double-walled host-guest complexes. Moreover, CPPs can also be stabilized by global aromaticity under twofold oxidation or reduction, involving electronic conjugation along with the overall structural backbone. Herein we explore the structural modifications, bonding, electron delocalization and magnetic properties of doubly reduced double-walled CPP complexes with DFT methods, in the isolated and aggregate [n + 5]CPP⊃[n]CPP^2- (n = 5-8) species. Our results show that the hosts undergo structural, bonding and delocalization deformations towards quinoidal configurations and exhibit global long-ranged shielding cones similar to global aromatic free dianionic CPPs, accounting for charge delocalization on the outer nanohoops, whereas the guests preserve local aromatic benzenoid configurations, resulting in global and local aromatic circuits within the host-guest aggregate. This observation suggests that in multi-layered related species electronic delocalization will be retained at the outer structural surface. The aromaticity of the hosts is manifested in the strong upfield shifts of the guests ¹H-NMR signals. Hence, CPP complexes can be extended to doubly reduced species stabilized by global host aromaticity expanding our understanding of doubled-walled nanotubes at the nanoscale regime.

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.

Local and global aromaticity under rotation: analysis of two- and three-dimensional representative carbon nanostructures

Nanoscaled 2D and 3D carbon structures with closed curved π-surfaces are of relevance in the development of desirable building units for materials science. Such species are able to sustain local and global aromatic circuits involving isolated regions or the overall structural backbone, respectively. Here we account for local and global aromaticity under rotation of representative two- and three-dimensional species involving para-connected and fused edge-sharing phenyl rings ([8]CPP, [10]CPP, CNB), and C₆₀ fullerene at different charge states. Our results denote that nanoscaled 2D global aromatics mimic the behaviour of the most prototypical aromatic 6π-circuit, given by benzene, where the shielding cone properties vary along the rotation motion. In contrast, 3D spherical aromatics remain almost invariant under rotation, given the distinctive characteristics of such species, differing from 2D global aromatics. Dissection of orbital contributions reveals that π-orbitals are determinants for shifting from non-aromatic to spherical aromatic species. Under rotation, the variation of the anisotropic effect inherent to such nanoscaled structures is accounted for, which is relevant to rationalize variation in NMR signal shifts upon the formation of host-guest aggregates.

On the antiaromatic-aromatic-antiaromatic transition of the stacked cyclobutadiene dimer

We have studied the changes in the aromatic nature of two cyclobutadiene (C₄H₄) molecules on decreasing the intermolecular distance and approaching the cubane structure in a face-to-face fashion. The analysis based on the calculations of the magnetically induced current density and the induced magnetic field shows that the aromaticity of the two C₄H₄ rings changes from a strongly antiaromatic character at long distances to an aromatic transition state of stacked C₄H₄ rings at intermediate internuclear distances when approaching the antiaromatic state of cubane.

Magnetic response properties of carbon nano-onions

The magnetic response of a number of double- and triple-layer carbon nano-onions (CNOs) is analyzed by calculating the magnetically induced current density and the induced magnetic field using the pseudo-π model. Qualitatively the same magnetic response was obtained in calculations at the all-electron level. The calculations show that the CNOs exhibit strong net diatropic (paratropic) ring currents when the external magnetic field points perpendicularly to one of the six-membered (five-membered) rings. They are deshielded inside and shielded outside the CNO; the latter dominates for larger CNOs. The magnetic response originates from a combination of spherical paratropic current densities on the inside of each carbon layer and diatropic current densities on the outside of them. The quantitative differences in the aromaticity of the CNOs as compared to single fullerenes are discussed in terms of ring-current strengths. The magnetic response of some of the CNOs is approximately the sum of the magnetic response of the individual layers, whereas deviations are significant for CNOs containing C₈₀. For the largest CNOs, the deviation from the sum of the fullerene contributions is larger, especially when the external magnetic field is perpendicular to a six-membered ring.

Core-electron contributions to the molecular magnetic response

Orbital contributions to the magnetic response depend on the method used to compute them. Here, we show that dissecting nuclear magnetic shielding tensors using natural localized molecular orbitals (NLMOs) leads to anomalous core contributions. The arbitrariness of the assignment might significantly affect the interpretation of the magnetic response of nonplanar molecules such as C₆₀ or [14]helicene and the assessment of their aromatic character. We solve this problem by computing the core- and σ-components of the induced magnetic field (and NICS) and the magnetically induced current density by removing the valence electrons (RVE). We estimate the core contributions to the magnetic response by performing calculations on the corresponding highly charged molecules, such as C₆H₆³⁰⁺ for benzene, using gauge-including atomic orbitals and canonical molecular orbitals (CMOs). The orbital contributions to nuclear magnetic shielding tensors are usually estimated by employing a natural chemical shielding (NCS) analysis in NLMO or CMO bases. The RVE approach shows that the core contribution to the magnetic response is small and localized at the nuclei, contrary to what NCS calculations suggest. This may lead to a completely incorrect interpretation of the magnetic σ-orbital response of nonplanar structures, which may play a major role in the overall magnetic shielding of the system. The RVE approach is thus a simple and inexpensive way to determine the magnetic response of the core- and σ-electrons.

Non-intersecting ring currents in [12]infinitene

The aromaticity of the newly synthesized [12]infinitene is addressed via analysis of the magnetically induced current density and the induced magnetic field. Our calculations reveal that [12]infinitene responds to an external magnetic field by creating two current-density pathways that flow diatropically along the edges of the molecule. The current-density pathways do not intersect. The entire structure is completely shielded suggesting that the infinitene molecule is aromatic, contrary to what the Möbius rule for twisted annulene structures predicts. We also studied the dication of [12]infinitene, which sustains two paratropic ring currents flowing along the edges. The space between the stacked rings at the crossing point is shorter for the dication as compared to the neutral molecule. Hence, a strong through-space current density appears at the crossing point of π-π stacked rings.

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.

Integration of global ring currents using the Ampère-Maxwell law

Magnetically induced ring currents are calculated from the magnetic shielding tensor by employing the Ampère-Maxwell law. The feasibility of the method is demonstrated by integrating the zz component of the shielding tensor along the symmetry axis of highly symmetric ring-shaped aromatic, antiaromatic and nonaromatic molecules. The calculated ring-current strengths agree perfectly with the ones obtained by integrating the current-density flux passing through a plane cutting half the molecular ring. The method can be used in combination with all electronic structure codes capable of calculating nuclear magnetic resonance (NMR) shielding tensors in general points in space. We also show that nucleus independent chemical shifts (NICS) along the symmetry axis are related to the spatial derivative of the strength of the global ring-current along the z axis.

The Magnetic Response of Starphenes

The aromaticity of [n]starphenes (n = 1, 4, 7, 10, 13, 16), as well as starphene-based [19]dendriphene, is addressed by calculating the magnetically induced current density and the induced magnetic field, using the pseudo-π model. When an external magnetic field is applied, these systems create diatropic currents that split into a global peripheral current surrounding the starphene skeleton and several local currents in the acene-based arms, resulting in large shielding cones above the arms. In particular, the arm currents are smaller than their linear analogs, and in general, the strengths of the ring currents tend to weaken as the starphene get larger.

Visualizing electron delocalization in contorted polycyclic aromatic hydrocarbons

Electron delocalization in contorted polycyclic aromatic hydrocarbon (PAH) molecules was examined through 3D isotropic magnetic shielding (IMS) contour maps built around the molecules using pseudo-van der Waals surfaces. The resulting maps of electron delocalization provided an intuitive, yet detailed and quantitative evaluation of the aromatic, non aromatic, and antiaromatic character of the local and global conjugated cyclic circuits distributed over the molecules. An attractive pictural feature of the 3D IMS contour maps is that they are reminiscent of the Clar π-sextet model of aromaticity. The difference in delocalization patterns between the two faces of the electron circuits in contorted PAHs was clearly visualized. For π-extended contorted PAHs, some splits of the π system resulted in recognizable patterns typical of smaller PAHs. The differences between the delocalization patterns of diastereomeric chiral PAHs could also be visualized. Mapping IMS on pseudo-van der Waals surfaces around contorted PAHs allowed visualization of their superimposed preferred circuits for electron delocalization and hence their local and global aromaticity patterns.

Electronic Currents and Anapolar Response Induced in Molecules by Monochromatic Light

It is shown that the electric dipole- and electric quadrupole–anapole polarizabilities, denoted respectively by fαβ′ and gα,βγ′, and the anapole magnetizability aαβ, are intrinsic properties of the electron cloud of molecules responding to optical fields. fαβ′ is a nonvanishing property of chiral and achiral compounds, whereas aαβ is suitable for enantiomer discrimination of chiral species. They can conveniently be evaluated by numerical integration, employing a formulation complementary to that provided by perturbation theory and relying on the preliminary computation of electronic current density tensors all over the molecular domain. The origin dependence of the dynamic anapolar response is rationalized via related computational techniques employing numerical integration, as well as definitions of molecular property tensors, for example, electric dipole and electric quadrupole polarizabilties and magnetizability. A preliminary application of the theory is reported for the Ra enantiomer of the hydrogen peroxide molecule, evaluating tensor components of electric dipole-anapole polarizability and anapole magnetizability as functions of the dihedral angle ϕ≡∠ H-O-O-H in the range 0∘≤ϕ≤180∘.