Probing the structures and bonding of size-selected boron and doped-boron clusters

Formation of Delocalized Linear M-B-M Covalent Bonds: A Combined Experimental and Theoretical Study of BM2(CO)8+ (M = Co, Rh, Ir) Complexes

Investigations of transition-metal boride clusters not only lead to novel structures but also provide important information about the metal-boron bonds that are critical to understanding the properties of boride materials. The geometric structures and bonding features of heteronuclear boron-containing transition metal carbonyl cluster cations BM(CO)6+ and BM2(CO)8+ (M = Co, Rh, and Ir) are studied by a combination of the infrared photodissociation spectroscopy and density functional calculations at B3LYP/def2-TZVP level. The completely coordinated BM2(CO)8+ complexes are characterized as a sandwich structure composed of two staggered M(CO)4 fragments and a boron cation, featuring a D3d symmetry and 1Eg electronic ground state as well as metal-anchored carbonyls in an end-on manner. In conjunction with theoretical calculations, multifold metal-boron-metal bonding interactions in BM2(CO)8+ complexes involving the filled d orbitals of the metals and the empty p orbitals of the boron cation were unveiled, namely, one σ-type M-B-M bond and two π-type M-B-M bonds. Accordingly, the BM2(CO)8+ complexes can be described as a linear conjugated (OC)4M═B═M(CO)4 skeleton with a formal B-M bond index of 1.5. The three delocalized d-p-d covalent bonds render compensation for the electron deficiency of the cationic boron center and endow both metal centers with the favorable 18-electron structure, thus contributing much to the overall structural stability of the BM2(CO)8+ cations. As a comparison, the saturated BRh(CO)6+ and BIr(CO)6+ complexes are determined to be a doublet Cs-symmetry structure with an unbridged (OC)2B-M(CO)4 pattern, involving a two-center σ-type (OC)2B → M(CO)4+ dative single bond along with a weak covalent B-M half bond. This work offers important insight into the structure and bonding of late transition metal boride carbonyl cluster cations.

Formation of Supernarrow Borophene Nanoribbons

Borophenes have sparked considerable interest owing to their fascinating physical characteristics and diverse polymorphism. However, borophene nanoribbons (BNRs) with widths less than 2 nm have not been achieved. Herein, we report the experimental realization of supernarrow BNRs. Combining scanning tunneling microscopy imaging with density functional theory modeling and ab initio molecular dynamics simulations, we demonstrate that, under the applied growth conditions, boron atoms can penetrate the outermost layer of Au(111) and form BNRs composed of a pair of zigzag (2,2) boron rows. The BNRs have a width self-contained to ∼1 nm and dipoles at the edges to keep them separated. They are embedded in the outermost Au layer and shielded on top by the evacuated Au atoms, free of the need for post-passivation. Scanning tunneling spectroscopy reveals distinct edge states, primarily attributed to the localized spin at the BNRs' zigzag edges. This work adds a new member to the boron material family and introduces a new physical feature to borophenes.

Stability and chemical bonding in a series of inverse sandwich actinide boride clusters (An2B8) with δ bonding

An inverse sandwich structure has been computationally predicted for uranium boride and extended to the series of actinide elements (An) from Th to Cm. The electronic structure and chemical bonding of these novel compounds have been analyzed using density functional theory and multireference wave-function based methods. We report the trends in electronic structure and bonding for An2B8, and found that (d-π)π and (d-p)δ are the most important factors in the stability of An2B8. The (f-p)δ bond provides extra stabilization for Pa2B8 and U2B8, owing to the extensive interactions of An-B8-An, resulting in a short distance for the Pa-Pa and U-U bonds.

MB16 - (M=Sc, Y, La): Perfect Bowl-Like Boron Clusters

The introduction of transition-metal doping has engendered a remarkable array of unprecedented boron motifs characterized by distinctive geometries and bonding, particularly those heretofore unobserved in pure boron clusters. In this study, we present a perfect (no defects) boron framework manifesting an inherently high-symmetry, bowl-like architecture, denoted as MB16 - (M=Sc, Y, La). In MB16 -, the B16 is coordinated to M atoms along the C5v-symmetry axis. The bowl-shaped MB16 - structure is predicted to be the lowest-energy structure with superior stability, owing to its concentric (2 π+10 π) dual π aromaticity. Notably, the C5v-symmetry bowl-like B16 - is profoundly stabilized through the doping of an M atom, facilitated by strong d-pπ interactions between M and boron motifs, in conjunction with additional electrostatic stabilization by an electron transfer from M to the boron motifs. This concerted interplay of covalent and electrostatic interactions between M and bowl-like B16 renders MB16 - a species of exceptional thermodynamic stability, thus making it a viable candidate for gas-phase experimental detection.

An Extended Rudolph Diagram: B3H5 and B3H6+ Relate 3D-, 2D-, 1D-, and 0D-Boron Allotropes

The structural chemistry of boron goes beyond the sp, sp2, and sp3 hybridization paradigms of carbon chemistry. We relate the apparently unconnected polyhedral boranes and 3D allotropes on the one hand and 2D clusters, borophenes, and multilayer borophenes on the other hand, through an extended Rudolph diagram. All-boron equivalents of cyclopropenium cation viz the flat B3H5 and the nonplanar B3H6+ constitute the missing links. The nonplanar B3H6+ (C3v) is the starting point for construction of polyhedral boranes; e.g., fusion of two of them leads to octahedral B6H62-. On the other hand, planar B3H6+ and B3H5 relate to borophenes with hexagonal holes. These borophene sheets can be further stacked with diverse interlayer BB bonds, ranging from bilayers to infinite layers. The tendency to achieve electron sufficiency as in the parent C3H3+ dictates the preference for hexagonal holes in the constituent layers and the interlayer bonds between them in multilayer borophenes. The design principles and theoretical validations for the formation of multilayer borophenes are also presented, indicating the variety and complexities involved.

Searching for stable copper borozene complexes in CuB7- and CuB8. Phys Chem Chem Phys 2024;26:12928-12938. [PMID: 38456623 DOI: 10.1039/d4cp00296b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Copper has been shown to be an important substrate for the growth of borophenes. Copper-boron binary clusters are ideal platforms to study the interactions between copper and boron, which may provide insight about the underlying growth mechanisms of borophene on copper substrates. Here we report a joint photoelectron spectroscopy and theoretical study on two copper-doped boron clusters, CuB7- and CuB8-. Well resolved photoelectron spectra are obtained for the two clusters at different wavelengths and are used to understand the structures and bonding properties of the two CuBn- clusters. We find that CuB8- is a highly stable borozene complex, which possesses a half-sandwich structure with a Cu+ species interacting with the doubly aromatic η8-B82- borozene. The CuB7- cluster is found to consist of a terminal copper atom bonded to a double-chain B7 motif, but it has a low-lying isomer composed of a half-sandwich structure with a Cu+ species interacting with an open-shell η7-B72- borozene. Both ionic and covalent interactions are found to be possible in the binary Cu-B clusters, resulting in different structures.

Electronic Control of the Position of the Pb Atom on the Surface of B8 Borozene in the PbB8 Cluster

Spontaneous symmetry-breaking is common in chemical and physical systems. Here, we show that by adding an electron to the C7v PbB8 cluster, which consists of a planar B8 disk with the Pb atom situated along the C7 axis, the Pb atom spontaneously moves to the off-axis position in the PbB8- anion. Photoelectron spectroscopy of PbB8- reveals a broad ground-state transition and a large energy gap, suggesting a highly stable closed-shell PbB8 borozene complex and a significant geometry change upon electron detachment. Quantum chemistry calculations indicate that the lowest unoccupied molecular orbital of the C7v PbB8 cluster is a degenerate π orbital mainly consisting of the Pb 6px and 6py atomic orbitals. Occupation of one of the 6p orbitals spontaneously break the C7v symmetry in the anion due to the Jahn-Teller effect. The large amplitude of the position change of Pb in PbB8- relative to PbB8 is surprising owing to bonding interactions between the Pb 6p orbital with the π orbital of the B8 borozene.

B63: The Most Stable Bilayer Structure with Dual Aromaticity

The emergence of a bilayer B48 cluster, which has been both theoretically predicted and experimentally observed, as well as the recent experimental synthesis of bilayer borophene sheets on Ag and Cu surfaces, has generated tremendous curiosity in the bilayer structures of boron clusters. However, the connection between bilayer boron cluster and bilayer borophene remains unknown. By combining a genetic algorithm and density functional theory calculations, a global search for the low-energy structures of the B63 cluster was conducted, revealing that the Cs bilayer structure with three interlayer B-B bonds is the most stable bilayer structure. This structure was further examined in terms of its structural stability, chemical bonding, and aromaticity. Interestingly, the interlayer bonds induce strong electronegativity and robust aromaticity. Furthermore, the dual aromaticity stems from diatropic currents originating from virtual translational transitions for both σ and π electrons. This unprecedent bilayer boron cluster is anticipated to enrich the concept of dual aromaticity and serve as a potential precursor for bilayer borophene.

Metal-Centered Boron-Wheel Cluster of Y©B112- with Rare D11h Symmetry

Molecules with high point-group symmetry are interesting prototype species in the textbook. As transition metal-centered boron clusters tend to have highly symmetric structures to fulfill multicenter bonding and high stability, new boron clusters with rare point-group symmetry may be viable. Through in-depth scrutiny over the structures of experimentally already observed transition metal-centered boron-wheel complexes, geometric and electronic design principles are summarized, based on which we studied M©B11k- (M = Y, La; Zr, Hf; k = 1, 2) clusters and found that a Y©B112- boron-wheel complex has an unprecedented D11h point-group symmetry. The remarkable stability of the planar Y©B112- complex is illustrated via extensive global-minimum structural search as well as comprehensive chemical bonding analyses. Similar to other boron-wheel complexes, the Y©B112- complex is shown to possess σ and π double aromaticity, indeed following the electronic design principle previously summarized. This new compound is expected to be experimentally identified, which will extend the currently known largest possible planar molecular symmetry and enrich the metal-centered boron-wheel class.

Size effects and electronic properties of zinc-doped boron clusters Zn B n (n = 1-15)

CONTEXT

To gain a deeper understanding of zinc-doped boron clusters, theoretical calculations were performed to investigate the size effects and electronic properties of zinc-doped boron clusters. The study of the electronic properties, spectral characteristics, and geometric structures of Zn B n (n = 1-15) is of great significance in the fields of semiconductor materials science, material detection, and improving catalytic efficiency. The results indicate that Zn B n (n = 1-15) clusters predominantly exhibit planar or quasi-planar structures, with the Zn atom positioned in the outer regions of the B n framework. The second stable structure of Zn B 3 is a three-dimensional configuration, indicating that the structures of zinc-doped boron clusters begin to convert from the planar or quasi-planar structures to the 3D configurations. The second low-energy structure of Zn B 15 is a novel configuration. Relative stability analyses show that the Zn B 12 has better chemical stability than other clusters with a HOMO-LUMO gap of 2.79 eV. Electric charge analysis shows that part electrons on zinc atoms are transferred to boron atoms, and electrons prefer to cluster near the B n framework. According to the electron localization function, it gets harder to localize electrons as the equivalent face value drops, and it's challenging to see covalent bond formation between zinc and boron atoms. The spectrograms of Zn B n (n = 1-15) exhibit distinct properties and notable spectral features, which can be used as a theoretical basis for the identification and confirmation of boron clusters doped with single-atom transition metals.

METHODS

The calculations were performed using the ABCluster global search technique combined with density functional theory (DFT) methods. The selected low-energy structures were subjected to geometric optimization and frequency calculations at the PBE0/6-311 + G(d) level to ensure structural stability and eliminate any imaginary frequencies. To acquire more precise relative energies, we performed single-point energies calculations for the low-lying isomers of Zn B n (n = 1-15) at the CCSD(T)/6-311 + G(d)//PBE0/6-311 + G(d) level of theory. All calculations were performed using Gaussian 09 software. To facilitate analysis, we utilized software tools such as Multiwfn, and VMD.

Infrared Photodissociation Spectroscopy of Mass-Selected Dinuclear Transition Metal Boride Carbonyl Cluster Cations

The transition-metal-boron bonding interactions and geometric structures of heterodinuclear transition metal carbonyl cluster cations BM(CO)n+ (M = Co, Ni, and Cu) are studied by a combination of the infrared photodissociation spectroscopy and density functional theory calculations at the B3LYP/def2-TZVP level. The BCu(CO)5+ and BCo(CO)6+ cations are characterized as an (CO)2B-M(CO)3/4+ structure involving an σ-type (OC)2B → M(CO)3,4+ dative bonding with end-on carbonyls, while for BNi(CO)5,6+ complexes with a bridged carbonyl, a 3c-2e bond involving the 5σ electrons of the bridged carbonyl and an electron-sharing bond between the B(CO)2 fragment and the Ni(CO)2,3+ subunits were revealed. Moreover, the fundamental driving force of the exclusive existence of a bridged carbonyl group in the boron-nickel complexes has been demonstrated to stem from the desire of the B and Ni centers for the favorable 8- and 18-electron structures.

Build Borophite from Borophenes: A Boron Analogue Graphite

Unlike graphene derived from graphite, borophenes represent a distinct class of synthetic two-dimensional materials devoid of analogous bulk-layered allotropes, leading to covalent bonding within borophenes instead of van der Waals (vdW) stacking. Our investigation focuses on 665 vdW-stacking boron bilayers to uncover potential bulk-layered boron allotropes through vdW stacking. Systematic high-throughput screening and stability analysis reveal a prevailing inclination toward covalently bonded layers in the majority of boron bilayers. However, an intriguing outlier emerges in δ5 borophene, demonstrating potential as a vdW-stacking candidate. We delve into electronic and topological structural similarities between δ5 borophene and graphene, shedding light on the structural integrity and stability of vdW-stacked boron structures across bilayers, multilayers, and bulk-layered allotropes. The δ5 borophene analogues exhibit metallic properties and characteristics of phonon-mediated superconductors, boasting a critical temperature near 22 K. This study paves the way for the concept of "borophite", a long-awaited boron analogue of graphite.

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.

Investigation of Pb-B Bonding in PbB2(BO)n- (n = 0-2): Transformation from Aromatic PbB2- to Pb[B2(BO)2]-/0 Complexes with BB Triple Bonds

Boron has been found to be able to form multiple bonds with lead. To probe Pb-B bonding, here we report an investigation of three Pb-doped boron clusters, PbB2-, PbB3O-, and PbB4O2-, which are produced by a laser ablation cluster source and characterized by photoelectron spectroscopy and ab initio calculations. The most stable structures of PbB2-, PbB3O-, and PbB4O2- are found to follow the formula, [PbB2(BO)n]- (n = 0-2), with zero, one, and two boronyl ligands coordinated to a triangular and aromatic PbB2 core, respectively. The PbB2- cluster contains a BB double bond and two Pb-B single bonds. The coordination of BO is observed to weaken Pb-B bonding but strengthen the BB bond in [PbB2(BO)n]- (n = 1, 2). The anionic [PbB2(BO)2]- and its corresponding neutral closed-shell [PbB2(BO)2] contain a BB triple bond. A low-lying Y-shaped isomer is also observed for PbB4O2-, consisting of a central sp2 hybridized B atom bonded to two boronyl ligands and a PbB unit.

Transition Metal Behavior of Heavier Alkaline Earth Elements in Doped Monocyclic and Tubular Boron Clusters

Quantum chemical calculations are carried out to design highly symmetric-doped boron clusters by employing the transition metal behavior of heavier alkaline earth (Ae = Ca, Sr, and Ba) metals. Following an electron counting rule, a set of monocyclic and tubular boron clusters capped by two heavier Ae metals were tested, which leads to the highly symmetric Ae2B8, Ae2B18, and Ae2B30 clusters as true minima on the potential energy surface having a monocyclic ring, two-ring tubular, and three-ring tubular boron motifs, respectively. Then, a thorough global minimum (GM) structural search reveals that a monocyclic B8 ring capped with two Ae atoms is indeed a GM for Ca2B8 and Ba2B8, while for Sr2B8 it is a low-lying isomer. Similarly, the present search also unambiguously shows the most stable isomers of Ae2B18 and Ae2B30 to be highly symmetric two- and three-ring tubular boron motifs, respectively, capped with two Ae atoms on each side of the tube. In these Ae-doped boron clusters, in addition to the electrostatic interactions, a substantial covalent interaction, specifically the bonding occurring between (n - 1)d orbitals of Ae and delocalized orbitals of boron motifs, provides the essential driving force behind their highly symmetrical structures and overall stability.

Planar Hexacoordinate Beryllium: Covalent Bonding Between s-block Metals

Achieving a planar hypercoordinate arrangement of s-block metals through covalent bonding with ligands is challenging due to the strong ionicity involved. Herein, we report the first case of a neutral binary global minimum containing a planar hexacoordinate beryllium atom. The central Be atom is coordinated by six active Be atoms, the latter in turn are enclosed by an equal number of more electronegative chlorine atoms in the periphery, forming a star-like phBe cluster (Be©Be6 Cl6 ). Importantly, the cluster exhibits dynamically stabilized stemming geometrically from the appropriate matching of metal-ligand size and electronically from adherence to the octet rule as well as possessing a 6σ/2π double aromaticity. Remarkably, energy decomposition analysis-natural orbitals for chemical valence (EDA-NOCV) analysis reveals a significant covalent interaction between the ligand and the central metal beryllium atoms, a fact further supported by a large Wiberg bond index. This cluster is a promising synthetic as its excellent electronic, dynamic and thermodynamic stability.

PB12+ and P2B12+/0/-: The Novel B12 Cage Doped by Nonmetallic P Atoms

A new kind of nonmetallic atom-doped boron cluster is described herein theoretically. When a phosphorus atom is added to the B12 motif and loses an electron, a novel B12 cage is obtained, composed of two B3 rings at both ends and one B6 ring in the middle, forming a triangular bifrustum. Interestingly, this B12 cage is formed by three B7 units joined together from three directions at an angle of 120°. When two P atoms are added to the B12 motif, this novel B12 cage is also obtained, and two P atoms are attached to the B3 rings at both ends of the triangular bifrustum, forming a triangular bipyramid (Johnson solid). Amazingly, the global minimums of neutral, monocationic, and monoanionic P2B12+/0/- have the same cage structure with a D3h symmetry; this is the smallest boron cage with the same structure. The P atom has five valence electrons, according to adaptive natural density partitioning bonding analyses of cage PB12+ and P2B12, in addition to one lone pair, the other three electrons of the P atom combine with an electron of each B atom on the B3 ring to form three 2c-2e σ bonds and form three electron sharing bonds with B atoms through covalent interactions, stabilizing the B12 cage. The calculated photoelectron spectra can be compared with future experimental values and provide a theoretical basis for the identification and confirmation of PnB12- (n = 1-2).

Semiconducting Bilayer Borophene with High Carrier Mobility

Borophene has attracted much interest due to its rich configurations and novel properties such as Dirac fermions and superconductivity. The recently emerged bilayer borophene mitigates the oxidation problem when exposed to air, yet most studies ignore the influence of charge transfer induced by metal substrates on structural stability. Here we identified 31 monolayer borophene polymorphs that are stabilized on Au(111), Ag(111), or Cu(111) substrates through first-principle calculations. Interestingly, two novel semiconducting bilayer borophene polymorphs with band gaps of 0.37 and 0.42 eV were screened by integrating these monolayers. The formation of interlayer bonding contributed by the delocalized electrons is responsible for the semiconductivity. The predicted highest electron mobility reaches 2.01 × 104 cm2V-1 s-1, implying the possibility as a semiconductor device with a low power consumption. Moreover, light was also systemically thrown on the factors that may affect the electronic properties of bilayer borophenes and the positional preference of interlayer bonds.

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.

Geometric and electronic structures of medium-sized boron clusters doped with plutonium

Doping metal heteroatoms is an effective strategy to regulate the geometric and electronic structure of boron based nanoclusters. However, the exploration of the ground state structures of metal-boron-based nanoclusters is still a challenge duo to the complexity of the bonding interactions between heterogeneous atoms and boron cluster and the number of isomers on the potential energy surface increases exponentially with cluster size. Here, we use the CALYPSO cluster structural search method in combination with density functional theory calculations to study the geometries and electronic properties of anionic boron clusters doped with plutonium (PuBn-,n= 10-20). Our results show that the medium-sized PuB14-cluster exhibits excellent stability with highest occupied molecular orbital-lowest unoccupied molecular orbital energy gap of 2.30 eV. The remarkable stability of the anionic PuB14-cluster is due to the robust interactions between the Pu metal and the B14skeleton, along with the strong covalent interactions between the B atoms. These findings enrich the geometric structure database of metal doped clusters and provide valuable insights for the future synthesis of boron based nanomaterials.

Structural Evolution and Electronic Properties of Two Sulfur Atom-Doped Boron Clusters

We present a theoretical study of structural evolution, electronic properties, and photoelectron spectra of two sulfur atom-doped boron clusters S2Bn0/- (n = 2-13), which reveal that the global minima of the S2Bn0/- (n = 2-13) clusters show an evolution from a linear-chain structure to a planar or quasi-planar structure. Some S-doped boron clusters have the skeleton of corresponding pure boron clusters; however, the addition of two sulfur atoms modified and improved some of the pure boron cluster structures. Boron is electron-deficient and boron clusters do not form linear chains. Here, two sulfur atom doping can adjust the pure boron clusters to a linear-chain structure (S2B20/-, S2B30/-, and S2B4-), a quasi-linear-chain structure (S2B6-), single- and double-chain structures (S2B6 and S2B9-), and double-chain structures (S2B5, and S2B9). In particular, the smallest linear-chain boron clusters S2B20/- are shown with an S atom attached to each end of B2. The S2B2 cluster possesses the largest highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap of 5.57 eV and the S2B2- cluster possesses the largest average binding energy Eb of 5.63 eV, which shows the superior chemical stability and relative stability, respectively. Interestingly, two S-atom doping can adjust the quasi-planar pure boron clusters (B7-, B10-, and B120/-) to a perfect planar structure. AdNDP bonding analyses reveal that linear S2B3 and planar SeB11- have π aromaticity and σ antiaromaticity; however, S2B2, planar S2B6, and planar S2B7- clusters have π antiaromaticity and σ aromaticity. Furthermore, AdNDP bonding analyses reveal that planar S2B4, S2B10, and S2B12 clusters are doubly (π and σ) aromatic, whereas S2B5-, S2B8, S2B9-, and S2B13- clusters are doubly (π and σ) antiaromatic. The electron localization function (ELF) analysis shows that S2Bn0/- (n = 2-13) clusters have different electron delocalization characteristics, and the spin density analysis shows that the open-shell clusters have different characteristics of electron spin distribution. The calculated photoelectron spectra indicate that S2Bn- (n = 2-13) have different characteristic peaks that can be compared with future experimental values and provide a theoretical basis for the identification and confirmation of these doped boron clusters. Our work enriches the new database of geometrical structures of doped boron clusters, provides new examples of aromaticity for doped boron clusters, and is promising to offer new ideas for nanomaterials and nanodevices.

A DFT study to investigate the physical, electrical, optical properties and thermodynamic functions of boron nanoclusters (MxB2n0; x=1,2, n=3,4,5)

First Principle DFT calculations employing the B3LYP/LanL2DZ/SDD level of theory were used to analyze the various characteristics of boron nanoclusters (B6, B8, and B10). These pure structures were further doped with four transition metals (Ta, Ti, Tc, and V) to examine the enhancement of the pure structures' structural, electrical, and optical features. To study structural stability, we have estimated cohesion energy and imaginary frequencies. Cohesion energies were entirely negative, with a range of -3.37 eV to -8.07 eV, and most constructions had no imaginary frequencies, indicating their structural occurrences. The calculated adsorption energy suggests that the order of stability of the pristine boron nanoclusters is B10>B8>B6, and TcB10 and Tc2B10 are the more stable structures. Mulliken charge, DOS, HOMO-LUMO, and the HOMO-LUMO gap have all been examined in-depth to provide insight into electrical characteristics. UV-Vis and CD measurements show the doped boron nanoclusters have excellent optical properties. Aside from calculating thermodynamic functions, we have also calculated the global DFT parameters, which give us a deep quantum mechanical understanding of the optimized structure for further research and applications in the field of science and technology.

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, V₂B₇^-, through computational global-minimum structure searches and quantum chemical calculations. This alloy cluster has a heptatomic boron ring as well as a perpendicular V₂ 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-B₂-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 V₂ 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.

Photoelectron Spectroscopy and Theoretical Study of Di-Copper-Boron Clusters: Cu2B3- and Cu2B4. J Phys Chem A 2023. [PMID: 37235389 DOI: 10.1021/acs.jpca.3c02417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Copper has been found to be able to mediate the formation of bilayer borophenes. Copper-boron binary clusters are ideal model systems to probe the copper-boron interactions, which are essential to understand the growth mechanisms of borophenes on copper substrates. Here, we report a joint photoelectron spectroscopy and theoretical study on two di-copper-doped boron clusters: Cu₂B₃^- and Cu₂B₄^-. Well-resolved photoelectron spectra are obtained, revealing the presence of a low-lying isomer in both cases. Theoretical calculations show that the global minimum of Cu₂B₃^- (C_2v, ¹A₁) contains a doubly aromatic B₃^- unit weakly interacting with a Cu₂ dimer, while the low-lying isomer (C_2v, ¹A₁) consists of a B₃ triangle with the two Cu atoms covalently bonded to two B atoms at two vertexes. The global minimum of Cu₂B₄^- (D_2h, ²A_g) is found to consist of a rhombus B₄ unit covalently bonded to the two Cu atoms at two opposite vertexes, whereas in the low-lying isomer (C_s, ²A'), one of the two Cu atoms is bonded to two B atoms.

Sc@B28-, Ti@B28, V@B28+, and V@B292-: Spherically Aromatic Endohedral Seashell-like Metallo-Borospherenes

Transition-metal-doped boron nanoclusters exhibit unique structures and bonding in chemistry. Using the experimentally observed seashell-like borospherenes C₂ B₂₈^-/0 and C_s B₂₉^- as ligands and based on extensive first-principles theory calculations, we predict herein a series of novel transition-metal-centered endohedral seashell-like metallo-borospherenes C₂ Sc@B₂₈^- (1), C₂ Ti@B₂₈ (2), C₂ V@B₂₈⁺ (3), and C_s V@B₂₉^2- (4) which, as the global minima of the complex systems, turn out to be the boron analogues of dibenzenechromium D_6h Cr(C₆H₆)₂ with two B₁₂ ligands on the top and bottom interconnected by four or five corner boron atoms on the waist and one transition-metal "pearl" sandwiched at the center in between. Detailed molecular orbital, adaptive natural density partitioning (AdNDP), and iso-chemical shielding surface (ICSS) analyses indicate that, similar to Cr(C₆H₆)₂, these endohedral seashell-like complexes follow the 18-electron rule in bonding patterns (1S²1P⁶1D¹⁰), rendering spherical aromaticity and extra stability to the systems.

Geometric and electronic diversity of metal doped boron clusters

Being intermediate between small compounds and bulk materials, nanoparticles possess unique properties different from those of atoms, molecules, and bulk matter. In the past two decades, a combination of cluster structure prediction algorithms and experimental spectroscopy techniques was successfully used for exploration of the ground-state structures of pure and metal-doped boron clusters. The fruitfulness of this dual approach is well illustrated by the discovery of intriguing microstructures and unique physicochemical properties such as aromaticity and bond fluxionality for both boron and metal-doped boron clusters. Our review starts with an overview of geometrical configurations of pure boron clusters B_n, which are presented by planar, nanotube, bilayer, fullerene-like and core-shell structures, in a wide range ofnvalues. We consider next recent advances in studies of boron clusters doped with metal atoms paying close and thoughtful attention to modifications of geometric and electronic structures of pure boron clusters by heteroatoms. Finally, we discuss the possibility of constructing boron-based nanomaterials with specific functions from metal-boron clusters. Despite a variety of fruitful results obtained in numerous studies of boron clusters, the exploration of boron-based chemistry has not yet reached its peak. The intensive research continues in this area, and it should be expected that it brings exciting discoveries of intriguing new structures.

Heptacoordinate transition-metal-decorated metallo-borospherenes and multiple-helix metallo-boronanotubes

The recent discovery of lanthanide-metal-decorated metallo-borospherenes LM₃B₁₈^- (LM = La, Tb) marks the onset of a new class of boron-metal binary nanomaterials. Using the experimentally observed or theoretically predicted borospherenes as ligands and based on extensive first-principles theory calculations, we predict herein a series of novel chiral metallo-borospherenes C₂ Ni₆ ∈ B₃₉^- (1), C₁ Ni₆ ∈ B₄₁⁺ (3), C₂ Ni₆ ∈ B₄₂²⁺ (4), C₂ Ni₆ ∈ B₄₂ (5), and C₂ Ni₈ ∈ B₅₆ (6) as the global minima of the systems decorated with quasi-planar heptacoordinate Ni (phNi) centers in η⁷-B₇ heptagons on the cage surfaces, which are found to be obviously better favoured in coordination energies than hexacoordinate Ni centers in previously reported D_2d Ni₆ ∈ B₄₀ (2). Detailed bonding analyses indicate that these phNi-decorated metallo-borospherenes follow the σ + π double delocalization bonding pattern, with two effective (d-p)σ coordination bonds formed between each phNi and its η⁷-B₇ ligand, rendering spherical aromaticity and extra stability to the systems. The structural motif in elongated axially chiral Ni₆ ∈ B₄₂²⁺ (4), Ni₆ ∈ B₄₂ (5), and Ni₈ ∈ B₅₆ (6) can be extended to form the metallic phNi-decorated boron double chain (BDC) double-helix Ni₄ ∈ B₂₈ (2, 0) (P4̄m2) (8), triple-helix Ni₆ ∈ B₄₂ (3, 0) (P3̄m1) (9), and quadruple-helix Ni₈ ∈ B₅₆ (4, 0) (P4mm) (10) metallo-boronanotubes, which can be viewed as quasi-multiple-helix DNAs composed of interconnected BDCs decorated with phNi centers in η⁷-B₇ heptagons on the tube surfaces in the atomic ratio of Ni : B = 1 : 7.

Perfect Core-Shell Octahedral B@B38 + , Be@B38 , and Zn@B38 with an Octa-Coordinate Center as Superatoms Following the Octet Rule

Planar, tubular, cage-like, and bilayer boron clusters B_n ^+/0/- (n=3∼48) have been observed in joint experimental and theoretical investigations in the past two decades. Based on extensive global searches augmented with first-principles theory calculations, we predict herein the smallest perfect core-shell octahedral borospherene O_h B@B₃₈ ⁺ (1) and its endohedral metallo-borospherene analogs O_h Be@B₃₈ (2), and O_h Zn@B₃₈ (3) which, with an octa-coordinate B, Be or Zn atom located exactly at the center, turn out to be the well-defined global minima of the systems highly stable both thermodynamically and dynamically. B@B₃₈ ⁺ (1) represents the first boron-containing molecule reported to date which contains an octa-coordinate B center covalently coordinated by eight face-capping boron atoms at the corners of a perfect cube in the first coordination sphere. Detailed natural bonding orbital (NBO) and adaptive natural density partitioning (AdNDP) bonding analyses indicate that these high-symmetry core-shell complexes X@B₃₈ ^+/0/- (X=B, Be, Zn) as super-noble gas atoms follow the octet rule in coordination bonding patterns (1S² 1P⁶ ), with one delocalized 9c-2e S-type coordination bond and three delocalized 39c-2e P-type coordination bonds formed between the octa-coordinate X center and its octahedral O_h B₃₈ ligand to effectively stabilize the systems. Their IR, Raman, and UV-Vis spectra are computationally simulated to facilitate their spectroscopic characterizations.

σ-Aromatic MAl6S6 (M = Ni, Pd, Pt) Stars Containing Planar Hexacoordinate Transition Metals

Hypercoordinate transition-metal species are mainly dominated by the 18-valence-electron (18ve) counting. Herein, we report ternary MAl₆S₆ (M = Ni, Pd, Pt) clusters with the planar hexacoordinate metal (phM) centers, which feature 16ve counting instead of the classic 18ve rule. These global-minimum clusters are established via unbiased global searches, followed by PBE0 and single-point CCSD(T) calculations. The phM MAl₆ units are stabilized by six peripheral bridging S atoms in these star-like species. Chemical bonding analyses reveal that there are 10 delocalized electrons around the phM center, which can render the aromaticity according to the (4n + 2) Hückel rule. It is worth noting that adding an (or two) electron(s) to its π-type lowest unoccupied molecular orbital (LUMO) will make the system unstable.

A plier-shaped binary molecular wheel B7Mg4+ cluster: hybrid in-plane heptacoordination, double π/σ aromaticity, and electronic transmutation

A plier-shaped charge-transfer [Mg2]2+[Mg2B7]− complex cluster exhibits double 6π/6σ aromaticity, whose hybrid molecular wheel structure is rationalized using the concept of electronic transmutation.

Structural Evolution and Electronic Properties of Selenium-Doped Boron Clusters SeBn0/- (n = 3-16)

A theoretical research of structural evolution, electronic properties, and photoelectron spectra of selenium-doped boron clusters SeB_n^0/- (n = 3-16) is performed using particle swarm optimization (CALYPSO) software in combination with density functional theory calculations. The lowest energy structures of SeB_n^0/- (n = 3-16) clusters tend to form quasi-planar or planar structures. Some selenium-doped boron clusters keep a skeleton of the corresponding pure boron clusters; however, the addition of a Se atom modified and improved some of the pure boron cluster structures. In particular, the Se atoms of SeB₇^-, SeB₈^-, SeB₁₀^-, and SeB₁₂^- are connected to the pure quasi-planar B₇^-, B₈^-, B₁₀^-, and B₁₂^- clusters, which leads to planar SeB₇^-, SeB₈^-, SeB₁₀^-, and SeB₁₂^-, respectively. Interestingly, the lowest energy structure of SeB₉^- is a three-dimensional mushroom-shaped structure, and the SeB₉^- cluster displays the largest HOMO-LUMO gap of 5.08 eV, which shows the superior chemical stability. Adaptive natural density partitioning (AdNDP) bonding analysis reveals that SeB₈ is doubly aromatic, with 6 delocalized π electrons and 6 delocalized σ electrons, whereas SeB₉^- is doubly antiaromatic, with 4 delocalized π electrons and 12 delocalized σ electrons. Similarly, quasi-planar SeB₁₂ is doubly aromatic, with 6 delocalized π electrons and 14 delocalized σ electrons. The electron localization function (ELF) analysis shows that SeB_n^0/- (n = 3-16) clusters have different local electron delocalization and whole electron delocalization effects. The simulated photoelectron spectra of SeB_n^- (n = 3-16) have different characteristic bands that can identify and confirm SeB_n^- (n = 3-16) combined with future experimental photoelectron spectra. Our research enriches the geometrical structures of small doped boron clusters and can offer insight for boron-based nanomaterials.

Theoretical Investigation of Geometries and Bonding of Indium Hydrides in the In2Hx and In3Hy (x = 0-4,6; y = 0-5) Series

Boron hydrides have been an object of intensive theoretical and experimental investigation for many decades due to their unusual and somewhat unique bonding patterns. Despite boron being a neighboring element to carbon, boron hydrides almost always form non-classical structures with multi-center bonds. However, we expect indium to form its interesting molecules with non-classical patterns, though such molecules still need to be extensively studied theoretically. In this work, we investigated indium hydrides of In₂H_x (x = 0-4,6) and In₃H_y (y = 0-5) series via DFT and ab initio quantum chemistry methods, performing a global minimum search, chemical bonding analysis, and studies of their thermodynamical stability. We found that the bonding pattern of indium hydrides differs from the classical structures composed of 1c-2e lone pairs and 2c-2e bonds and the bonding pattern of earlier investigated boron hydrides of the B_nH_n+2 series. The studied stoichiometries are characterized by multi-center bonds, aromaticity, and the tendency for indium to preserve the 1c-2e lone pair.

Actinide-doped boron clusters: from borophenes to borospherenes

Similar to graphene and fullerene, metal-doping has been considered to be an effective approach to the construction of highly stable boron clusters. In this work, a series of actinide metal-doped boron clusters AnB₃₆ (An = Pa, Np, Pu, Am, Cm, Bk, and Cf) have been explored using extensive first-principles calculations. We found that the quasi-planar structure of B₃₆ transforms to an endohedral borospherene An@B₃₆ after actinide metal doping. Actinoborospherenes exhibit C_2h symmetry with Pa, Np, and Pu dopants and for Am, Cm, Bk and Cf dopants with larger atomic radii, the symmetry of An@B₃₆ is reduced to C_i. Bonding property analyses such as bond order, molecular orbital (MO) and quantum theory of atoms in molecules (QTAIM) analysis show that the covalency of the An-B bonds in C_2h An@B₃₆ (An = Pa, Np, and Pu) is higher than that in C_i An@B₃₆ (An = Am, Cm, Bk, and Cf). These endohedral borospherenes are robust according to thermodynamic and dynamic analyses. As expected, the C_i An@B₃₆ clusters are less stable compared to C_2h An@B₃₆, which is consistent with the stronger covalent bonds of the latter. These results indicate that the existence of the actinide-boron bonding is essential for the high stability of the An@B₃₆ clusters, confirming that the fullerene-like boron cages can be stabilized by actinide encapsulation. This work is expected to provide potential routes for the construction of robust borospherenes.

Unprecedented Prediction of a B160 Cluster Stuffed by Dual-Icosahedron B12

Serving as the premise to understand bulk allotropes, boron clusters have been intriguing experimentalists and theoreticians to study their geometries and chemical bonding. Here, we designed a complete core-shell B₁₆₀ cluster stuffed by two B₁₂ cores, which is energetically preferable over the bilayer structure of the same size. The unprecedented peanutlike structure with C_i symmetry has superior stability and exhibits superatomic electronic configuration and spherical aromaticity. Our theoretical work not only proposed the core-shell structure of dual icosahedrons for the first time but also indicated the multi-B₁₂ core-shell structural pattern in boron particles, bridging to boron crystalline structures.

Superatomic icosahedral-CnB12-n (n = 0, 1, 2) Stuffed mononuclear and binuclear borafullerene and borospherene nanoclusters with spherical aromaticity

Boron and boron-based nanoclusters exhibit unique structural and bonding patterns in chemistry. Extensive density functional theory calculations performed in this work predict the mononuclear walnut-like Ci C50B54 (1) (C2B10@C48B44), C1 C50B54 (2) (CB11@C49B43), and S10 C50B54 (3) (B12@C50B42) which contain one icosahedral-CnB12-n core (n = 0, 1, 2) at the center following the Wade’s skeletal electron counting rules and the approximately electron sufficient binuclear peanut-like Cs C88B78 (4) ((C2B10)2@C84B58), Cs C88B78 (5) ((CB11)2@C86B56), Cs C88B78 (6) ((B12)2@C88B54), Cs B180 (7) ((B12)2@B156), Cs B182 (8) ((B12)2@B158), and Cs B184 (9) ((B12)2@B160) which encapsulate two interconnected CnB12-n icosahedrons inside. These novel core–shell borafullerene and borospherene nanoclusters appear to be the most stable species in thermodynamics in the corresponding cluster size ranges reported to date. Detailed bonding analyses indicate that the icosahedral B122−, CB11−, and C2B10 cores in these core–shell structures possess the superatomic electronic configuration of 1S21P61D101F8, rendering spherical aromaticity and extra stability to the systems. Such superatomic icosahedral-CnB12-n stuffed borafullerenes and borospherenes with spherical aromaticity may serve as embryos to form bulk boron allotropes and their carbon-boron binary counterparts in bottom-up approaches.

Synthesis and isolation of dinuclear N,C-chelate organoboron compounds bridged by neutral, anionic, and dianionic 4,4'-bipyridine via reductive coupling of pyridines

The reaction of Bppy(Mes)2 (BN1; ppy = 2-phenylpyridine) and BCH₂ppy(Mes)₂ (BN3) with the reducing reagent KC₈ resulted in C-C bond formation via intermolecular radical coupling to generate the 4,4'-bipyridyl ligand compounds BN2 and BN4. Adding 1 equivalent of KC₈ to a THF solution of BN2 and BN4 generated the 4,4'-bipyridyl radical anions BN2K and BN4K. The dianion species BN2K2 and BN4K2 could be obtained by adding 2 equivalents of KC₈ to the THF solution of BN2 and BN4. In the presence of 2,2,2-cryptand or 18-crown-6, the radical anion salt BN2K(crypt) and the dianion salt BN2K2(18c6)2 were isolated for single-crystal X-ray diffraction analysis. Structural, spectroscopic, and computational studies were performed on the three species of BN2 derivatives (neutral, radical anion, and dianion species). BN2 and BN4 were stable and did not undergo photoisomerization or photoelimination under UV light irradiation.

Structural evolution and relative stability of vanadium-doped boron clusters

Cluster is the intermediate of individual atom and larger agglomeration. The structural evolutions of clusters are critically important to explore the physical properties of bulk solids. Here, we carry out systematic structure predictions of medium-sized vanadium-doped boron clusters by using crystal structure analysis by particle swarm optimization method combined with density function theory calculations. A great deal of low-lying isomers with attractive geometries are discovered, such as the crown-like VB₁₈^-cluster and the drum-like VB₂₀^-cluster. Interestingly, the VB₁₂^-cluster possesses excellently relative stability due to its higher second-order difference and larger highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) energy gap. The molecular orbitals (MOs) and adaptive natural density partitioning (AdNDP) analysis indicate that the 3dorbitals of V atom and the 2pand 2sorbitals of B atoms are the primary constituents of the MOs, and the interactions between V and B atoms are the main factor for the robust stabilization of the anionic VB₁₂^-cluster. The present findings advance the understanding of the structural evolution of transition metal doped boron clusters and offer crucial insights for future experiments.

Lanthanide/actinide boride nanoclusters and nanomaterials based on boron frameworks consisting of conjoined Bn rings (n = 7-9)

Extensive global minimum searches augmented with first-principles theory calculations performed in this work indicate that the experimentally observed perfect inverse sandwich lanthanide boride complexes D_7h La₂B₇^- (1), D_8h La₂B₈ (3), D_9h La₂B₉^- (7) can be extended to their actinide counterparts C_2v Ac₂B₇^- (1'), D_8h Ac₂B₈ (3'), D_9h Ac₂B₉^- (7') with a B_n monocyclic ring (n = 7-9) sandwiched by two Ac dopants. Such M₂B_n^-/0 inverse sandwiches (1/1', 3/3', 7/7') can be used as building blocks to generate the ground-state C₂ La₄B₁₃^- (2)/Ac₄B₁₃^- (2'), D₂ La₄B₁₅^- (4)/Ac₄B₁₅^- (4'), C_3v/C₃ La₄B₁₈ (5)/Ac₄B₁₈ (5'), O_h Ac₇B₂₄⁺ (6'), O_h Ac₇B₂₄, T_d Ac₄B₂₄ (8'), C₁ La₅B₂₄⁺ (9)/Ac₅B₂₄⁺ (9'), and T_d Ac₄B₂₉^- (10') which are based on boron frameworks consisting of multiple conjoined B_n rings (n = 7-9). Detailed bonding analyses show that effective (d-p)σ, (d-p)π and (d-p)δ coordination bonds are formed between the B_n rings and metal doping centers, conferring three-dimensional aromaticity and extra stability to the systems. In particular, the perfect body-centered cubic O_h Ac₇B₂₄⁺ (6') and O_h Ac₇B₂₄ with six conjoined B₈ rings can be extended in x, y, and z dimensions to form one-dimensional Ac₁₀B₃₂ (11'), two-dimensional Ac₃B₁₀ (12'), and three-dimensional AcB₆ (13') nanomaterials, presenting a B₈-based bottom-up approach from metal boride nanoclusters to their low-dimensional nanomaterials.

A bottom-up approach from medium-sized bilayer boron nanoclusters to bilayer borophene nanomaterials

Inspired by the experimentally observed bilayer B₄₈^-/0 and theoretically predicted bilayer B₅₀-B₇₂ and based on extensive density functional theory calculations, we report herein a series of novel medium-sized bilayer boron nanoclusters C₁ B₈₄ (I), C_2v B₈₆ (II), C₁ B₈₈ (III), C₁ B₉₀ (IV), C₁ B₉₂ (V), C₁ B₉₄ (VI), C_2v B₉₆ (VII), and C₁ B₉₈ (VIII) which are the most stable isomers of the systems reported to date effectively stabilized by optimum numbers of interlayer B-B σ bonds between the inward-buckled atoms on top and bottom layers. Detailed bonding analyses indicate that these bilayer species follow the universal bonding pattern of σ + π double delocalization, rendering three-dimensional aromaticity in the systems. More interestingly, the AA-stacked bilayer structural motif in B₉₆ (VII) with a B₇₂ bilayer hexagonal prism at the center can be extended to form bilayer C₂ B₁₂₈ (IX), D_2h B₂₁₄ (X), C_2v B₂₆₀ (XI), D_2h B₃₇₂ (XII), and D₂ B₈₂₈ (XIII) which contain one or multiple conjoined B₇₂ bilayer hexagonal prisms sharing interwoven zig-zag boron triple chains between them. Such bilayer species or their close-lying AB isomers can be viewed as embryos of the newly reported most stable freestanding BL-α⁺ bilayer borophenes and quasi-freestanding bilayer borophenes on Ag(111) which are composed of interwoven zig-zag boron triple chains shared by conjoined BL B₇₂ hexagonal prisms, presenting a bottom-up approach from medium-sized bilayer boron nanoclusters to two-dimensional bilayer borophene nanomaterials.

Construction of the Largest Metal-Centered Double-Ring Tubular Boron Clusters Based on Actinide Metal Doping

Metal doping has been considered to be an effective approach to stabilize various boron clusters. In this work, we constructed a series of largest metal-centered double-ring tubular boron clusters An@B₂₄ (An = Th, Pa, Pu, and Am). Extensive global minimum structural searches combined with density functional theory predicted that the global minima of An@B₂₄ (An = Th, Pu, and Am) are double-ring tubular structures. Formation energy analysis indicates that these boron clusters are highly stable, especially for Th@B₂₄ and Pa@B₂₄. Detailed bonding analysis shows that the significant stability of An@B₂₄ is determined by the covalent character of the An-B bonding, which stems from the interactions of An 5f and 6d orbitals and B 2p orbitals. These results show that actinide metal doping is a feasible route to construct stable large metal-centered double-ring tubular boron clusters, offering the possibility to design boron nanomaterials with special physiochemical properties.

B96: a complete core-shell structure with high symmetry

A complete core-shell and highly symmetric B₉₆ was designed. The core-shell B₉₆ of T_h symmetry is energetically favorable compared to the bilayer motif and the core-shell structure can be well maintained during first-principles molecular dynamics simulations at high temperatures (up to 1000 K). Moreover, it exhibits a superatomic electronic configuration and spherical aromaticity. Our theoretical work not only confirmed that the core-shell structural pattern is more energetically favorable for large-sized boron clusters, but also provided a strategy to design large boron clusters with a core-shell structure of high symmetry.

Structures, and electronic and spectral properties of single-atom transition metal-doped boron clusters MB24 - (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni)

A theoretical study of geometrical structures, electronic properties, and spectral properties of single-atom transition metal-doped boron clusters MB₂₄ ^- (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) is performed using the CALYPSO approach for the global minimum search, followed by density functional theory calculations. The global minima obtained for the MB₂₄ ^- (M = Sc, Ti, V, and Cr) clusters correspond to cage structures, and the MB₂₄ ^- (M = Mn, Fe, and Co) clusters have similar distorted four-ring tubes with six boron atoms each. Interestingly, the global minima obtained for the NiB₂₄ ^- cluster tend to a quasi-planar structure. Charge population analyses and valence electron density analyses reveal that almost one electron on the transition-metal atoms transfers to the boron atoms. The electron localization function (ELF) of MB₂₄ ^- (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) indicates that the local delocalization of MB₂₄ ^- (M = Sc, Ti, V, Cr, and Ni) is weaker than that of MB₂₄ ^- (M = Mn, Fe, and Co), and there is no obvious covalent bond between doped metal and B atoms. The spin density and spin population analyses reveal that open-shell MB₂₄ ^- (M = Ti, Cr, Fe, and Ni) has different spin characteristics which are expected to lead to interesting magnetic properties and potential applications in molecular devices. The polarizability of MB₂₄ ^- (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) shows that MB₂₄ ^- (M = Mn, Fe, and Co) has larger first hyperpolarizability, indicating that MB₂₄ ^- (M = Mn, Fe, and Co) has a strong nonlinear optical response. Hence, MB₂₄ ^- (M = Mn, Fe, and Co) might be considered as a promising nonlinear optical boron-based nanomaterial. The calculated spectra indicate that MB₂₄ ^- (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) has different and meaningful characteristic peaks that can be compared with future experimental values and provide a theoretical basis for the identification and confirmation of these single-atom transition metal-doped boron clusters. Our work enriches the database of geometrical structures of doped boron clusters and can provide an insight into new doped boron clusters.