Axially Modified Square-Pyramidal CoN4-F1 Sites Enabling High-Performance Zn-Air Batteries.

Cobalt-nitrogen-carbon (Co-N-C) catalysts with a CoN4 structure exhibit great potential for oxygen reduction reaction (ORR), but the imperfect adsorption energy toward oxygen species greatly limits their reduction efficiency and practical application potential. Here, F-coordinated Co-N-C catalysts with square-pyramidal CoN4-F1 configuration are successfully synthesized using F atoms to regulate the axial coordination of Co centers via hydrothermal and chemical vapor deposition methods. During the synthesis process, the geometry structure of the Co atom converts from six-coordinated Co-F6 to square-pyramidal CoN4-F1 in the coordinatively unsaturated state, which provides an open binding site for the O2. The introduction of axial F atoms into the CoN4 plane alters the local atomic environment around Co, significantly improving the ORR activity and Zn-air batteries performance. In situ spectroscopy proves that CoN4-F1 sites strongly combine with the OOH* intermediate and facilitate the splitting of O-O bond, making OOH* readily decompose into O* and OH* via a dissociative pathway. Theoretical calculations confirm that the axial F atom effectively reduces the electronic density of the Co centers and facilitates the desorption of the OH* intermediate, efficiently accelerating the overall ORR kinetics. This work advances a feasible synthesis mechanism of axial ligands and provides a route to construct efficient high-coordination catalysts.

Boron-based ternary MgTa2B6 cluster: a turning nanoclock with dynamic structural fluxionality.

Boron-based complex clusters are a fertile ground for the exploration of exotic chemical bonding and dynamic structural fluxionality. Here we report on the computational design of a ternary MgTa2B6 cluster via global structural searches and quantum chemical calculations. The cluster turns out to be a new member of the molecular rotor family, closely mimicking a turning clock at the subnanoscale. It is composed of a hexagonal B6 ring with a capping Ta atom at the top and bottom, whereas the Mg atom is linked to one Ta site as a radial Ta-Mg dimer. These components serve as the dial, axis, and hand of a nanoclock, respectively. Chemical bonding analyses reveal that the inverse sandwich Ta2B6 motif in the cluster features 6π/6σ double aromaticity, whose electron counting conforms to the (4n + 2) Hückel rule. The Ta-Mg dimer has a Lewis-type σ bond, and the Mg site has negligible bonding with B6 ring. The ternary cluster can be formulated as an [Mg]0[Ta2B6]0 complex. Molecular dynamics simulations suggest that the cluster is structurally fluxional analogous to a nanoclock, even at a low temperature of 100 K. The Ta-Mg hand turns almost freely around the Ta2 axis and along the B6 dial. The tiny intramolecular rotation barrier is less than 0.3 kcal mol-1, being dictated by the bonding nature of double 6π/6σ aromaticity. The present system offers a new type of molecular rotor in physical chemistry.

Boron-based Pd3B26 alloy cluster as a nanoscale antifriction bearing system: tubular core-shell structure, double π/σ aromaticity, and dynamic structural fluxionality.

Boron-based nanoclusters show unique geometric structures, nonclassical chemical bonding, and dynamic structural fluxionality. We report here on the theoretical prediction of a binary Pd3B26 cluster, which is composed of a triangular Pd3 core and a tubular double-ring B26 unit in a coaxial fashion, as identified through global structural searches and electronic structure calculations. Molecular dynamics simulations indicate that in the core-shell alloy cluster, the B26 double-ring unit can rotate freely around its Pd3 core at room temperature and beyond. The intramolecular rotation is virtually barrier free, thus giving rise to an antifriction bearing system (or ball bearing) at the nanoscale. The dimension of the dynamic system is only 0.66 nm. Chemical bonding analysis reveals that Pd3B26 cluster possesses double 14π/14σ aromaticity, following the (4n + 2) Hückel rule. Among 54 pairs of valence electrons in the cluster, the overwhelming majority are spatially isolated from each other and situated on either the B26 tube or the Pd3 core. Only one pair of electrons are primarily responsible for chemical bonding between the tube and the core, which greatly weaken the bonding within the Pd3 core and offers structural flexibility. This is a key mechanism that effectively diminishes the intramolecular rotation barrier and facilitates dynamic structural fluxionality of the system. The current work enriches the field of nanorotors and nanomachines.

Boron-Based Inverse Sandwich V2B7- Cluster: Double π/σ Aromaticity, Metal-Metal Bonding, and Chemical Analogy to Planar Hypercoordinate Molecular Wheels.

Inverse sandwich clusters composed of a monocyclic boron ring and two capping transition metal atoms are interesting alloy cluster systems, yet their chemical bonding nature has not been sufficiently elucidated to date. We report herein on the theoretical prediction of a new example of boron-based inverse sandwich alloy clusters, V2B7-, through computational global-minimum structure searches and quantum chemical calculations. This alloy cluster has a heptatomic boron ring as well as a perpendicular V2 dimer unit that penetrates through the ring. Chemical bonding analysis suggests that the inverse sandwich cluster is governed by globally delocalized 6π and 6σ frameworks, that is, double 6π/6σ aromaticity following the (4n + 2) Hückel rule. The skeleton B-B σ bonding in the cluster is shown not to be strictly Lewis-type two-center two-electron (2c-2e) σ bonds. Rather, these are quasi-Lewis-type, roof-like 4c-2e V-B2-V σ bonds, which amount to seven in total and cover the whole surface of inverse sandwich in a truly three-dimensional manner. Theoretical evidence is revealed for a 2c-2e Lewis σ single bond within the V2 dimer. Direct metal-metal bonding is scarce in inverse sandwich alloy clusters. The present inverse sandwich alloy cluster also offers a new type of electronic transmutation in physical chemistry, which helps establish an intriguing chemical analogy between inverse sandwich clusters and planar hypercoordinate molecular wheels.

Chemical Bonding and Dynamic Structural Fluxionality of a Boron-Based Na5B7 Sandwich Cluster.

Doping alkali metals into boron clusters can effectively compensate for the intrinsic electron deficiency of boron and lead to interesting boron-based binary clusters, owing to the small electronegativity of the former elements. We report on the computational design of a three-layered sandwich cluster, Na5B7, on the basis of global-minimum (GM) searches and electronic structure calculations. It is shown that the Na5B7 cluster can be described as a charge-transfer complex: [Na4]2+[B7]3-[Na]+. In this sandwich cluster, the [B7]3- core assumes a molecular wheel in shape and features in-plane hexagonal coordination. The magic 6π/6σ double aromaticity underlies the stability of the [B7]3- molecular wheel, following the (4n + 2) Hückel rule. The tetrahedral Na4 ligand in the sandwich has a [Na4]2+ charge-state, which is the simplest example of three-dimensional aromaticity, spherical aromaticity, or superatom. Its 2σ electron counting renders σ aromaticity for the ligand. Overall, the sandwich cluster has three-fold 6π/6σ/2σ aromaticity. Molecular dynamics simulation shows that the sandwich cluster is dynamically fluxional even at room temperature, with a negligible energy barrier for intramolecular twisting between the B7 wheel and the Na4 ligand. The Na5B7 cluster offers a new example for dynamic structural fluxionality in molecular systems.

NAl4X4 + (X = S, Se, Te): Clusters with a planar tetracoordinate nitrogen and significantly improved stability.

The design of clusters featuring non-classical planar hypercoordinate atoms (phAs) often depends on the delocalized multicenter bonds involving reactive electron-deficient elements, which both destabilize the clusters and lead to difficulty in achieving the phA arrangement for electronegative elements such as nitrogen due to their preference for localized bonds. In this work, we computationally designed a series of aluminum chalcogenide clusters NAl4X4 + (X = S, Se, Te) with a desired planar tetracoordinate nitrogen and meaningfully improved chemical stability, as evidenced by the wide gaps (6.51-7.23 eV) between their highest occupied molecular orbitals and lowest unoccupied molecular orbitals, high molecular rigidity (dynamically stable up to 1500 K), and exclusively low global energy minima nature (their isomers locate at least 51.2 kcal/mol higher). Remarkably, these clusters are stabilized by peripheral chalcogen atoms, which not only sterically protect the NAl4 core moiety but also electronically compensate for the electron-deficient aluminum atoms via X → Al π back bonds, meeting the description of our recently proposed "electron-compensation" strategy.

Chemical bonding and dynamic structural fluxionality of a boron-based Al2B8 binary cluster: the robustness of a doubly 6π/6σ aromatic [B8]2- molecular wheel.

Despite the isovalency between Al and B elements, Al-doping in boron clusters can deviate substantially from an isoelectronic substitution process. We report herein on a unique sandwich di-Al-doped boron cluster, Al2B8, using global structural searches and quantum chemical calculations. The cluster features a perfectly planar B8 molecular wheel, with two isolated Al atoms symmetrically floating above and below it. The two Al atoms are offset from the center of the molecular wheel, resulting in a C 2v symmetry for the cluster. The Al2B8 cluster is shown to be dynamically fluxional even at far below room temperature (100 K), in which a vertical Al2 rod slides or rotates freely within a circular rail on the B8 plate, although there is no direct Al-Al interaction. The energy barrier for intramolecular rotation is only 0.01 kcal mol-1 at the single-point CCSD(T) level. Chemical bonding analysis shows that the cluster is a charge-transfer complex and can be formulated as [Al]+[B8]2-[Al]+. The [B8]2- molecular wheel in sandwich cluster has magic 6π/6σ double aromaticity, which underlies the dynamic fluxionality, despite strong electrostatic interactions between the [Al]+, [B8]2-, and [Al]+ layers.

Be3B11- cluster: a dynamically fluxional beryllo-borospherene.

The beryllium-doped Be3B11- cluster has two nearly isoenergetic isomers, adopting the smallest trihedral spherical geometries with a boron single-chain skeleton. The B11 skeleton in the global minimum (C2v, 1A1) comprises three conjoined boron rings (one B8/two B7) on the waist, sharing two B3 equilateral triangles at the top and bottom, respectively. However, the local minimum (Cs, 1A') has one deformed B4 pyramid at the top. The drastic structural transformation of B11 skeletons from perfectly planar B11 clusters mainly profited from robust electrostatic interaction between Be atoms and B11 skeletons. The dynamic simulations suggest that two species can interconvert via a novel mechanism, that is "triangle-pyramid-triangle", which facilitates the free migration of boron atoms in the B11 skeleton, thereby showing the fascinating dynamic fluxionality. The chemical bonding analyses reveal that the B11 skeleton is covered by two types of delocalized π bonds in an orthogonal direction, which leads to its spherical aromaticity.

Y©B8C4 cluster: a boron-carbon molecular wheel with dodeca-coordination number in plane.

Computational evidence is reported for the largest planar molecular wheel of the Y©B8C4 cluster, featuring an yttrium atom enclosed by a highly symmetric B8C4 ring. The B8C4 ring is viable in the -(BCB)4- form with double 9π/10σ aromaticity. The centered yttrium atom is dodeca-coordinated with the peripheral B8C4 ring, which sets a record coordination number for a planar structure in chemistry heretofore.

Bare and ligand protected planar hexacoordinate silicon in SiSb3M3 + (M = Ca, Sr, Ba) clusters.

The occurrence of planar hexacoordination is very rare in main group elements. We report here a class of clusters containing a planar hexacoordinate silicon (phSi) atom with the formula SiSb3M3 + (M = Ca, Sr, Ba), which have D 3h (1A1') symmetry in their global minimum structure. The unique ability of heavier alkaline-earth atoms to use their vacant d atomic orbitals in bonding effectively stabilizes the peripheral ring and is responsible for covalent interaction with the Si center. Although the interaction between Si and Sb is significantly stronger than the Si-M one, sizable stabilization energies (-27.4 to -35.4 kcal mol-1) also originated from the combined electrostatic and covalent attraction between Si and M centers. The lighter homologues, SiE3M3 + (E = N, P, As; M = Ca, Sr, Ba) clusters, also possess similar D 3h symmetric structures as the global minima. However, the repulsive electrostatic interaction between Si and M dominates over covalent attraction making the Si-M contacts repulsive in nature. Most interestingly, the planarity of the phSi core and the attractive nature of all the six contacts of phSi are maintained in N-heterocyclic carbene (NHC) and benzene (Bz) bound SiSb3M3(NHC)6 + and SiSb3M3(Bz)6 + (M = Ca, Sr, Ba) complexes. Therefore, bare and ligand-protected SiSb3M3 + clusters are suitable candidates for gas-phase detection and large-scale synthesis, respectively.

Boron Oxide B5O6 - Cluster as a Boronyl-Based Inorganic Analog of Phenolate Anion.

Boron oxide clusters have structural richness and exotic chemical bonding. We report a quantum chemical study on the binary B5O6 - cluster, which is relatively oxygen-rich. A global structural search reveals planar C 2v (1A1) geometry as the global minimum structure, featuring a heteroatomic hexagonal B3O3 ring as its core. The three unsaturated B sites are terminated by two boronyl (BO) groups and an O- ligand. The B5O6 - cluster can be faithfully formulated as B3O3(BO)2O-. This structure is in stark contrast to that of its predecessors, C s B5O5 - and T d B5O4 -, both of which have a tetrahedral B center. Thus, there exists a major structural transformation in B5O n - series upon oxidation, indicating intriguing competition between tetrahedral and heterocyclic structures. The chemical bonding analyses show weak 6π aromaticity in the B5O6 - cluster, rendering it a boronyl analog of phenolate anion (C6H5O-) or boronyl boroxine. The calculated vertical detachment energy of B5O6 - cluster is 5.26 eV at PBE0, which greatly surpasses the electron affinities of halogens (Cl: 3.61 eV), suggesting that the cluster belongs to superhalogen anions.

Structures and chemical bonding of boron-based B12O and B11Au clusters. A counterexample in boronyl chemistry.

Boron oxide clusters have structural diversity and unique chemical bonding, and recent literature has shown that boronyl complexes dominate boron-rich oxide clusters. A counterexample in boronyl chemistry is presented in this work. Using global structural searches, electronic structure calculations, and chemical bonding analyses, we shall report on the computational design of two boron-based quasi-planar or planar clusters: B12O and B11Au. Contrary to expectation, the B12O cluster has a circular quasi-planar shape with a peripheral B-O-B bridge, which resembles bare B12 cluster. It does not contain a boronyl ligand. The isomeric boronyl complex turns out to be 10.32 kcal mol-1 higher in energy at the single-point CCSD(T) level. In contrast, B11Au cluster behaves normally with an elongated B11 moiety and a terminal Au ligand. Chemical bonding analyses reveal three-fold π/σ aromaticity in circular B12O cluster, including global 6π aromaticity, as well as spatially isolated inner 2σ aromaticity and outer 10σ aromaticity. The three-fold 6π/2σ/10σ aromaticity underlies the stability of B12O cluster. This bonding picture is unknown for bare B12 cluster and its derivatives. The elongated B11Au cluster has conflicting π/σ aromaticity (with 6π versus 8σ electron-counting). The B12O cluster is actually isoelectronic with bare B12 cluster in terms of delocalized π/σ bonding, which inherits the structural and electronic robustness of the latter.

The unique sandwich K6Be2B6H6 cluster with a real borozene B6H6 core.

Theoretical evidence is reported for a boron-based K6Be2B6H6 sandwich cluster, showing a perfectly D 6h B6H6 ring, being capped by two tetrahedral K3Be ligands. Due to the comfortable charge transfer, the sandwich is viable in [K3Be]3+[B6H6]6-[BeK3]3+ ionic complex in nature. The [B6H6]6- core with 6π aromaticity vividly imitates the benzene (C6H6), occurring as a real borozene. In contrast, the tetrahedral [K3Be]3+ ligand is 2σ three-dimensional aromatic, acting as the simple superatom. Thus, this complex possesses a collectively three-fold 2σ/6π/2σ aromaticity. The interlaminar interaction is governed by the robust electrostatic attraction. The unique chemical bonding gives rise to interesting dynamic fluxionality.

Ternary 14-electron XB2Be2 (X = Si, Ge, Sn, Pb) clusters: a planar tetracoordinate silicon (ptSi) system and its ptGe/Sn/Pb congeners.

A class of ternary 14-electron clusters, XB2Be2 (X = Si, Ge, Sn, Pb), have been computationally predicted with a planar tetracoordinate silicon (ptSi) unit, as well as its heavier ptGe/Sn/Pb congeners. These pentaatomic ptSi/Ge/Sn/Pb species are established as global-minimum structures via computer global searches, followed by electronic structure calculations at the PBE0-D3, B3LYP-D3, and single-point CCSD(T) levels. Molecular dynamics simulations indicate that they are also kinetically stable against isomerization or decomposition. Chemical bonding analyses show that the clusters have double 2π/2σ aromaticity. The latter concept underlies the stability of ptSi/Ge/Sn/Pb clusters, overriding the 14-electron count or its variants, such as the 18-electron rule. No sp3 hybridization occurs in these species, which naturally explains why they are ptSi/Ge/Sn/Pb (rather than traditional tetrahedral) systems.

Planar pentacoordinate carbon in a sulphur-surrounded boron wheel: the global minimum of CB5S5.

We report a σ + π double aromatic CB5S5+ cluster, the first global minimum unusually having a planar hypercoordinate carbon inside a boron wheel. Five peripheral sulfur atoms stabilize the carbon-centered boron wheel by weakening the electron deficiency of the boron atoms through strong S → B π back-bonding.

Transition-metal-like bonding behaviors of a boron atom in a boron-cluster boronyl complex [(η7-B7)-B-BO].

Boron displays many unusual structural and bonding properties due to its electron deficiency. Here we show that a boron atom in a boron monoxide cluster (B9O-) exhibits transition-metal-like properties. Temperature-dependent photoelectron spectroscopy provided evidence of the existence of two isomers for B9O-: the main isomer has an adiabatic detachment energy (ADE) of 4.19 eV and a higher energy isomer with an ADE of 3.59 eV. The global minimum of B9O- is found surprisingly to be an umbrella-like structure (C 6v, 1A1) and its simulated spectrum agrees well with that of the main isomer observed. A low-lying isomer (C s, 1A') consisting of a BO unit bonded to a disk-like B8 cluster agrees well with the 3.59 eV ADE species. The unexpected umbrella-like global minimum of B9O- can be viewed as a central boron atom coordinated by a η7-B7 ligand on one side and a BO ligand on the other side, [(η7-B7)-B-BO]-. The central B atom is found to share its valence electrons with the B7 unit to fulfill double aromaticity, similar to that in half-sandwich [(η7-B7)-Zn-CO]- or [(η7-B7)-Fe(CO)3]- transition-metal complexes. The ability of boron to form a half-sandwich complex with an aromatic ligand, a prototypical property of transition metals, brings out new metallomimetic properties of boron.

Concentric Inner 2π/6σ and Outer 10π/14σ Aromaticity Underlies the Dynamic Structural Fluxionality of Planar B19- Wankel Motor Cluster.

Planar C2v B19- global-minimum (GM) cluster is known as a molecular Wankel motor, featuring unique chemical bonding and structural fluxionality. While the geometry, bonding, and molecular dynamics of the cluster are documented in the literature, it remains warranted to fully understand its bonding nature and unravel the mechanism behind the structural dynamics. We shall offer herein an updated bonding model on the bases of canonical molecular orbital (CMO) analysis and adaptive natural density partitioning (AdNDP), further aided by natural bond orbital (NBO) analysis and orbital composition calculations. The computational data indicate that the B19- cluster has inner 2π/6σ and outer 10π/14σ concentric 4-fold π/σ aromaticity. Being spatially isolated from each other, the inner B6 disk supports 2π and 6σ subsystems, whereas the outer B18 double-ring ribbon has 10π and 14σ subsystems. All 4-fold π/σ subsystems are intrinsically delocalized and conform to the (4n + 2) Hückel rule for aromaticity. The change of Wiberg bond index (WBI) from GM to transition-state (TS) for radial B-B links is minimal and uniform, which offers a semiquantitative measure of structural dynamics and underlies the low energy barrier.

Anchoring a bow-shaped boron single chain in binary Be6B7- cluster: hybrid octagonal ring, multifold π/σ aromaticity, and dual electronic transmutation.

Elemental boron clusters do not form linear chain or monocyclic ring structures, which is in contrast to carbon. Based on computer global searches and quantum chemical calculations, we report on the viability of a curved boron single chain in binary Be6B7- cluster. The boron motif assumes a bow shape, being anchored on a Be6 prism. Such a motif, which appears to be highly strained in its free-standing form, is exotic in boron-based clusters and nanostructures. Chemically, the cluster is analogous to a "clam-and-pearl-chain" system at the nanoscale (about 1 nm in size), in which a Be6 clam moderately opens its mouth, except that a B7 pearl chain is too large to be encapsulated inside. The picture differs from a three-layered sandwich. This cluster features a hybrid Be2B7 monocyclic ring, which is octagonal in nature and supports double 10π/6σ aromaticity. The number of π bonds substantially surpasses that in bare boron clusters of similar sizes. Two Be3 rings in the prism are also σ aromatic, albeit with effective 1σ/1σ electron-counting only. The unique multifold 1σ/10π/6σ/1σ aromaticity governs the geometry of the Be6B7- cluster, which can also be rationalized using the concept of dual electronic transmutation.

Are all planar and quasi-planar boron clusters aromatic? Counter examples of island or global π antiaromaticity from chemical bonding analysis.

Boron is an electron-deficient element. The flatland of planar or quasi-planar (2D) boron clusters is believed to possess aromaticity for all members, which remains a fundamental issue in debate in boron chemistry. Using a selected set of D2h B62-, C2h B282-, and C2v B29- clusters as counter examples, we shall present computational evidence for global or island π antiaromaticity in 2D boron clusters. The latter two are flattened for the purpose of clarity, which model their quasi-planar C2 or Cs monoanion clusters observed in prior gas-phase experiments. Chemical bonding in the clusters is elucidated collectively on the basis of canonical molecular orbital (CMO) analysis, adaptive natural density partitioning (AdNDP), electron localization functions (ELFs), and localized molecular orbital (LMO) analysis. These results are complementary to each other and yet highly coherent. As a quantitative indicator, nucleus-independent chemical shifts (NICSs) are calculated at selected specific points in the clusters, which help differentiate between π aromaticity and antiaromaticity. Intriguingly, triangular sites in the same boron cluster can be aromatic, antiaromatic, or nonaromatic, despite the fact that they are physically indistinguishable. The phenomenon is understood in analogy to hydrocarbons and polycyclic aromatic hydrocarbons (PAHs). Even perfect sheet-like boron clusters are convertible to the PAH analogous systems. This work provides compelling examples for global and island π antiaromaticity in the 2D boron clusters.

Boron-based ternary Rb6Be2B6 cluster featuring unique sandwich geometry and a naked hexagonal boron ring.

Computational evidence is reported on a boron-based ternary Rb6Be2B6 cluster as the "Big Mac" sandwich on a subnanoscale with thickness of 0.58 nm. The core hexagonal B6 ring, occurring in the naked form due to double 6π/6σ aromaticity, is capped by two tetrahedral BeRb3 ligands. Such a B6 motif is scarce in boron clusters. The sandwich cluster has four-fold 2σ/6π/6σ/2σ aromaticity and its tetrahedral BeRb3 ligand is the simplest case of three-dimensional aromaticity (or spherical aromaticity). The sandwich can be formulated as a charge-transfer complex, [Rb3Be]3+[B6]6-[BeRb3]3+, whose components are held together by robust electrostatics, facilitating dual-mode dynamic fluxionality.