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.

InnTl4-nH+ (n = 0∼4): Tetracoordinate Hydrogen in a Planar Fashion?

The recent report of planar tetracoordinate hydrogen (ptH) in In4H+ is very intriguing in planar hypercoordinate chemistry. Our high-level CCSD(T) calculations revealed that the proposed D4h-symmetric ptH In4H+ is a first-order saddle point with an imaginary frequency in the out-of-plane mode of the hydrogen atom. In fact, at the CCSD(T)/aug-cc-pV5Z/aug-cc-pV5Z-PP level, the C4v isomer, with the H atom located 0.70 Å above the In4 plane, is 0.5 kcal/mol more stable than the D4h isomer. However, given the small perturbation from planarity and essentially barrierless C4v ↔ D4h ↔ C4v transition, the vibrationally averaged structure can still be considered as a planar. Extending our exploration to the InnTl4-nH+ (n = 0-3) systems, we found all these ptH structures, except for In2Tl2H+, to be the putative global minimum. The single σ-delocalized interaction between the central hydrogen atom and InnTl4-n ligand rings proves pivotal in establishing planarity and aromaticity and conferring substantial stability upon these rule-breaking ptH species.

Electronic Structure and Aromaticity of an Unusual Cyclo[18]carbon Precursor, C18 Br6.

Herein, the electronic structure and bonding character of the stable cyclo[18]carbon (C18 ) precursor, C18 Br6 , are thoroughly characterized by molecular orbital (MO), density of states (DOS), bond order (BO), and interaction region indicator (IRI) analyses. The delocalization characters of out-of-plane and in-plane π-electrons (labeled as πout - and πin -electrons, respectively) in bonding regions were examined using localized orbital locator (LOL) and electron localization function (ELF). The aromaticity was investigated, studying the molecular magnetic response to external magnetic field by computing the magnetically induced current density (Jind ), iso-chemical shielding surface (ICSS), anisotropy of the induced current density (AICD), and the induced magnetic field (Bind ). All these analyses indicate that C18 Br6 is a globally aromatic species with lower aromaticity than C18 , and the blocking of in-plane π-conjugation (labeled as πin -conjugation) by the presence of -Br substituents in it is the underlying cause for the weakening of molecular aromaticity.

Exploring the Aromaticity Differences of Isoelectronic Species of Cyclo[18]carbon (C18), B6C6N6, and B9N9: The Role of Carbon Atoms as Connecting Bridges.

The cyclo[18]carbon (C18) has piqued widespread interest in recent years for its geometrical aesthetic and unique electronic structure. Inspired by it, theoretical investigations of its isoelectronic B9N9 have been published occasionally; however, few studies considered their other companion B6C6N6. In this work, we study the geometric structure, charge distribution, bonding characteristic, aromaticity, and electron delocalization of B6C6N6 and B9N9 for the first time and compare the relevant results with those of C18. Based on the comprehensive analysis of aromaticity indicators such as AV1245, nucleus-independent chemical shifts, anisotropy of the induced current density, magnetically induced current density, iso-chemical shielding surface, and induced magnetic field (Bind), we found that B6C6N6 has definitely a double aromatic character similar to C18 and the aromaticities of the two are very close, while B9N9 is a nonaromatic species. In response to this novel finding, we delved into its nature from an electron delocalization perspective through a localized orbital locator, electron localization function, Fermi hole, and atomic remote delocalization index analyses. The C atom between B and N as an interconnecting bridge strengthens the electron delocalization of the conjugate path, which is the essence of the significant enhancement of the molecular aromaticity from B9N9 to B6C6N6. This work elucidates that within the framework of the isoelectronicity of C18, different methods of atomic doping can achieve molecules with completely different properties.

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

We have studied the changes in the aromatic nature of two cyclobutadiene (C4H4) 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 C4H4 rings changes from a strongly antiaromatic character at long distances to an aromatic transition state of stacked C4H4 rings at intermediate internuclear distances when approaching the antiaromatic state of cubane.

The magnetically induced current density of the [12]infinitene dianion.

The magnetically induced current-density susceptibility (MICD) of the [12]infinitene dianion and the induced magnetic field around it have been calculated at the density functional theory level. Separation of the MICD into diatropic and paratropic contributions shows that the MICD is dominated by the diatropic contribution in contrast to the notion that it is antiaromatic, which was reported in a recently published article. The MICD of the [12]infinitene dianion exhibits several through-space MICD pathways, whereas it sustains only weak local paratropic current-density contributions. We identified four main current-density pathways of which two are similar to the ones for neutral [12]infinitene. It is difficult to judge from calculations of the nucleus independent shielding constants and the induced magnetic field around the [12]infinitene dianion whether it sustains diatropic or paratropic ring currents.

Correction: Structure and bonding of molecular stirrers with formula B7M2- and B8M2 (M = Zn, Cd, Hg).

Correction for 'Structure and bonding of molecular stirrers with formula B7M2- and B8M2 (M = Zn, Cd, Hg)' by Rui Yu et al., Phys. Chem. Chem. Phys., 2020, 22, 12312-12320, https://doi.org/10.1039/D0CP01603A.

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 C80. 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.

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.

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 C60 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 C6H630+ 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.

B3Al4+: A Three-Dimensional Molecular Reuleaux Triangle.

We systematically explore the potential energy surface of the B3Al4+ combination of atoms. The putative global minimum corresponds to a structure formed by an Al4 square facing a B3 triangle. Interestingly, the dynamical behavior can be described as a Reuleaux molecular triangle since it involves the rotation of the B3 triangle at the top of the Al4 square. The molecular dynamics simulations, corroborating with the very small rotational barriers of the B3 triangle, show its nearly free rotation on the Al4 ring, confirming the fluxional character of the cluster. Moreover, while the chemical bonding analysis suggests that the multicenter interaction between the two fragments determines its fluxionality, the magnetic response analysis reveals this cluster as a true and fully three-dimensional aromatic system.

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

The pseudo-π model yields current densities and induced magnetic fields that mimic the π-component, allowing investigations of large molecular structures, whether they are planar or not, at a low computational cost but with high accuracy. Herein the π-contribution to the magnetically induced current densities and induced magnetic fields of large planar molecules and nonplanar molecules (such as [10]cyclophenacene and chiral toroidal nanotubes C2016 and C2196) were computed using the pseudo-π model with the gauge-including magnetically induced currents method. Additionally, we provide a way to determine the π-component of the ring-current strengths, which can be used for assessing the aromatic character of large carbon molecules.

Planar Tetracoordinate Carbons in Allene-Type Structures.

The exhaustive exploration of the potential energy surfaces of CE2M2 (E = Si-Pb; M = Li and Na) revealed seven global minima containing a planar tetracoordinate carbon (ptC). The design, based on a π-localization strategy, resulted in a ptC with two double bonds forming a linear or a bent allene-type E═C═E motif. The magnetic response of the bent E═C═E fragments support a σ-aromaticity. The bonding analysis indicated that the ptCs form C-E covalent bonds and C-M electrostatic interactions.

Consequences of Curvature on Induced Magnetic Field: The Case of Helicenes.

Helicenes consist of several fused rings twisted around an axis, forming a cylindrical helix, with π-delocalized electrons in the non-planar rings. Induced magnetic fields dissecting the orbital contributions of [6]-, [7]-, and [14]helicene are discussed. Computations show a deshielding cone produced by the π-electrons along the helical axis. Unexpectedly, the response of the core electrons produces a shielding cone, which is cumulative and sensitive to the curvature of the systems owing to the overlap of the other ring responses. A warning is provided regarding the evaluation of the delocalization in curved systems in which the x- and y-components of the induced magnetic field become relevant.

Structure and bonding of molecular stirrers with formula B7M2- and B8M2 (M = Zn, Cd, Hg).

In this work, we systematically explored clusters with formula B7M2- and B8M2 (M = Zn, Cd, Hg). The putative global minima are formed by an M2 dimer and a disk-shaped boron wheel. Moreover, the chemical bonding analysis revealed that charge transfer from the metal atoms to the boron motifs resulted in (B7)3-(M2)2+ and (B8)2-(M2)2+ complexes with double (σ + π) aromatic boron wheels and a single bond for the metallic dimer. Above all, the computed rotational barriers of the M-M fragment with respect to the boron disk and molecular dynamics simulations indicate a virtually barrierless spin, resembling a magnetic stirrer on a baseplate.

Molecular Helmholtz coils.

How to build a molecular Helmholtz coil? The possibility to create a Helmholtz coil at the molecular level is studied via the induced magnetic responses of several small cyclic hydrocarbon dimers with formula (CnHn)2 and Dnh symmetry (n = 6-10). Our results reveal that for n ≥ 8, π-electrons give rise to a uniform magnetic field within the central region between rings, satisfying the Helmholtz coil condition. This uniformity is independent of the intensity of the induced magnetic field.

Planar pentacoordinate carbon in CGa5+ derivatives.

We report a family of systems having a planar pentacoordinate carbon (ppC) based on the next heavier analogue of CAl5+, the ppC system par excellence. Although because of the larger size of Ga, the ppC isomer is not even a local minimum in CGa5+, a single isoelectronic substitution of Ga by smaller sized Be maximizes the bonding in the ppC form. Retaining the 18 valence electron rule, the global minimum structures of CGa4Be, CGa3Be2-, CGa2Be32-, and CGaBe43- clusters and their corresponding lithium salts have a ppC.

E3 M3+ (E=C-Pb, M=Li-Cs) Clusters: The Smallest Molecular Stars.

Extensive potential energy surface explorations of twenty-five clusters with the formula E3 M3+ (E=Group 14 element and M=Group 1 element) through density functional theory and high-level ab initio computations reveal that the lowest-energy isomer for all these systems corresponds to a non-classical D3h star-like structure in the singlet state, where three M atoms interact electrostatically with the triangular E3 core, occupying three bridging positions around it. More than 18 200 calculations were done in the search for the minima structures, starting with a first phase at the PBE0/LANL2DZ level and ending with an analysis of the most representative clusters at the CCSD(T)/def2-TZVP//PBE0/def2-TZVP level. The title clusters represent the smallest molecular stars with three planar tetracoordinate E atoms (E=Group 14 element). All these E3 M3+ clusters behave like superalkali cations with small vertical electron affinities (smaller than Cs), large vertical electron detachment energies, and HOMO-LUMO energy gaps. Their energetics, bonding, and electron delocalization are discussed in detail. The high stability of these clusters is reflected from the large dissociation energy needed for different dissociation channels. The electron delocalization is confirmed by the presence of two delocalized π electrons over the E3 core and strong diatropic responses.

CAl11-: a molecular rotor with a quasi-planar tetracoordinate carbon.

In this work, we analyzed the bonding and fluxional character of the global minimum of CAl11-. Its structure is formed by two stacked layers, one of them resembles the well-known planar tetracoordinate carbon CAl4 on top of a hexagonal Al@Al6 wheel. Our results show that the CAl4 fragment rotates freely around the central axis. The exceptional stability and fluxionality of CAl11- derive from its particular electron distribution.

Planar pentacoordinate s-block metals.

The presence of a delocalized π-bond is often considered an essential criterion for achieving planar hypercoordination. Herein, we show that σ-delocalization could be sufficient to make the planar configuration the most stable isomer in a series of planar pentacoordinate s-block metals. High-level ab initio computations reveal that the global minimum of a series of interalkali and interalkali-alkaline earth clusters (LiNa5, Li5Mg+, Na5Mg+, K5Ca+, CaRb5+, Rb5Sr+, and SrCs5+) adopts a singlet D5h structure with a planar pentacoordinate lithium or alkaline earth metal (AE = Mg, Ca, Sr). These clusters are unusual combinations to stabilize a planar pentacoordinate atom, as all their constituents are electropositive. Despite the absence of π-electrons, Hückel's rule is fulfilled by the six σ-electrons. Furthermore, the systems exhibit a diatropic ring current in response to an external magnetic field and a strong magnetic shielding, so they might be classified as σ-aromatic. Therefore, multicenter σ-bonds and the resulting σ-delocalization stabilize these clusters, even though they lack π-aromaticity.