Doped Boron Clusters Research Articles

Geometric and electronic properties of anionic precious metal gold doped boron clusters AuB n − (n = 10–20)

Recently, the Au–B covalent bonds in gold doped boron clusters has attracted great attention. However, there are fewer theoretical reports on exploration their ground state structures and stabilities, especially for the medium sizes. Here, we study the structural evolution and electronic properties of the anionic Au doped boron clusters with medium sizes of n from 10 to 20 using the unbiased cluster structural searches combined with density functional theory (DFT) calculations. The results reveal that the quasi-planar AuB18 − (1A, C 1) cluster shows excellent stability and a large vertical separation energy (VDE) of 4.25 eV. The good consistency between the computationally simulated photoelectron spectra and the experimental spectra strongly supports the correctness of our low-lying structures. Further bonding analyses show that the well-stabilized aromatic AuB18 − cluster is due to the active σ interactions between Au atom (6s orbitals) and B units (2p orbitals), as well as the large number of σ–bonds in the B18 − moiety with π-aromaticity. These findings enriched the family of Au-B alloy clusters and metal-doped boron-based aromatic clusters, which provide valuable information for the experimental characterization and preparation of boron-rich alloy nanoclusters in the future.

Structural Evolution of Small-Sized Phosphorus-Doped Boron Clusters: A Half-Sandwich-Structured PB15 Cluster.

The present study is a theoretical investigation into the structural evolution, electronic properties, and photoelectron spectra of phosphorus-doped boron clusters PBn0/- (n = 3-17). The results of this study revealed that the lowest energy structures of PBn- (n = 3-17) clusters, except for PB17-, exhibit planar or quasi-planar structures. The lowest energy structures of PBn (n = 3-17), with the exceptions of PB7, PB9, and PB15, are planar or quasi-planar. The ground state of PB7 has an umbrella-shaped structure, with C6V symmetry. Interestingly, the neutral cluster PB15 has a half-sandwich-like structure, in which the P atom is attached to three B atoms at one end of the sandwich, exhibiting excellent relative and chemical stability due to its higher second-order energy difference and larger HOMO-LUMO energy gap of 4.31 eV. Subsequently, adaptive natural density partitioning (AdNDP) and electron localization function (ELF) analyses demonstrate the bonding characteristics of PB7 and PB15, providing support for the validity of their stability. The calculated photoelectron spectra show distinct characteristic peaks of PBn- (n = 3-17) clusters, thus providing theoretical evidence for the future identification of doped boron clusters. In summary, our work has significant implications for understanding the structural evolution of doped boron clusters PBn0/- (n = 3-17), motivating further experiments regarding doped boron clusters.

Geometry and electronic properties of alkali metal (rubidium) doped boron clusters

AbstractAlkali metal‐doped boron clusters have captured much attention because of their novel electronic properties and structural evolution. In the study of RbBn0/− (n = 2–12) clusters, the minimum global search of the potential energy surface and structure optimization at the level of PBE1PBE by using the CALYPSO method and Gaussian package coupled with DFT calculation; the geometrical structures and electronic properties are systematically investigated. At n = 8, the ground‐state structures are composed of an Rb atom above B atoms, forming a structurally stable pagoda cone. By stability analysis and charge transfer calculation, the RbB8− cluster shows more stability. It found that s‐p hybridization between Rb atom and B atoms as well as s‐p hybridization between B atoms is one of the reasons for the outstanding stability exhibited in the RbB80/− clusters by using DOS and HOMO–LUMO orbital contour maps. The chemical bonding of the RbB80/− groups was analyzed by using the AdNDP method, and B atoms with larger numbers readily form multi‐center chemical bonds with the Rb atom. From the results of the bonding analysis, the interaction between the Rb atom and B atoms strengthens the stability of the RbB80/− clusters. It is hoped that this work provides a direction for experimental manipulation.

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.

Reversible hydrogen storage in Ti decorated small boron clusters: Insights from molecular dynamics simulations

We report the hydrogen adsorption capacity of Ti doped boron clusters (BnTi2, n = 6–10) using dispersion corrected density functional theory. The stability and the reactivity of the studied clusters were confirmed by their reactivity indices. It was shown that each titanium atom can accommodate up to 5H2 molecules with average adsorption energy falling in the range of 0.15–0.21 eV/H2 via Niu-Rao-Jena kind of interaction. An elongation in H–H bond length was found to be 0.76–0.80 Å with average Ti–H2 bond distance being 2.754 Å which confirmed molecular physisorption of H2. The molecular dynamics simulation demonstrated that, these clusters could adsorb all H2 molecules at low temperature (77K) while desorbing almost all the H2 molecules at 300K and 373K. The occupation number calculation revealed that these clusters could hold up to 10.20 wt % of H2 molecules at storage condition (100K-60bar).

Inquiring into geometric structures and electronic properties of sodium doped boron clusters: DFT study of NaB n (n = 1–12) clusters

The systems of alkali metal atom doped boron have been investigative topics over the years because of their diverse configurations and unique properties. Based on this, with the assistance of the CALYPSO software facilitates the structure prediction process, our group went into the geometric structures and bonding characteristics of the NaB n (n = 1–12) clusters. From the structural analysis, the participation of the Na atom does not dramatically change the configurations of the clusters, and the NaB n (n = 1–12) clusters maintain the parallel planar configuration evolution model with the pure boron clusters. The local most stable structures emerge at n = 3, 5, 8, the NEC result shows electrons transfer from Na-3s and B-2s orbitals to B-2p orbitals, indicating the s-p hybridisation exists between Na-B, B-B atoms, and the DOS analysis also prove this point. In the subsequent bonding analysis of the NaB3 and NaB8 clusters, the high proportions of B-2s, B-2p, or Na-3s orbitals at HOMO, LUMO, demonstrate the formation of strong B-B, Na-B covalent bonds, which is the benefit to forming high stability. Last, the calculations of ICSS show, perpendicular to the boron ring plane direction, the NaB3 exhibits strong aromaticity while the NaB8 shows anti-aromaticity.

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 SeBn0/- (n = 3-16) is performed using particle swarm optimization (CALYPSO) software in combination with density functional theory calculations. The lowest energy structures of SeBn0/- (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 SeB7-, SeB8-, SeB10-, and SeB12- are connected to the pure quasi-planar B7-, B8-, B10-, and B12- clusters, which leads to planar SeB7-, SeB8-, SeB10-, and SeB12-, respectively. Interestingly, the lowest energy structure of SeB9- is a three-dimensional mushroom-shaped structure, and the SeB9- 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 SeB8 is doubly aromatic, with 6 delocalized π electrons and 6 delocalized σ electrons, whereas SeB9- is doubly antiaromatic, with 4 delocalized π electrons and 12 delocalized σ electrons. Similarly, quasi-planar SeB12 is doubly aromatic, with 6 delocalized π electrons and 14 delocalized σ electrons. The electron localization function (ELF) analysis shows that SeBn0/- (n = 3-16) clusters have different local electron delocalization and whole electron delocalization effects. The simulated photoelectron spectra of SeBn- (n = 3-16) have different characteristic bands that can identify and confirm SeBn- (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.

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 VB18 − cluster and the drum-like VB20 − cluster. Interestingly, the VB12 − 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 3d orbitals of V atom and the 2p and 2s orbitals 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 VB12 − cluster. The present findings advance the understanding of the structural evolution of transition metal doped boron clusters and offer crucial insights for future experiments.

Two-layer molecular rotors: A zinc dimer rotating over planar hypercoordinate motifs.

Multi-layer molecular rotors represent a class of unique combination of topology and bonding, featuring a barrier-free rotation of one layer with respect to other layers. This emerging fluxional behavior has been found in a few doped boron clusters. Herein, we strongly enrich this intriguing family followed by an effective design strategy, summarized as essential factors: i) considerable electrostatic interactions originated from a strong charge transfer between layers; ii) the absence of strong covalent bonds between layers; and iii) fully delocalized σ/π electrons from at least one layer. We found that planar hypercoordinate motifs consisting of monocyclic boron rings and metals with σ + π dual aromaticity can be regarded as one promising layer, which can support the suspended X2 (X=Zn, Cd, Hg) dimers. By detailed investigations of thermodynamic and kinetic stabilities of 60 species, eventually, MB7 X2 - and MB8 X2 (X=Zn, Cd; M=Be, Ru, Os; Be works only for Zn-based cases) clusters were verified to be the global-minimum two-layer molecular rotors. Especially, their electronic structure analyses vividly confirm the practicability of the electronic structure requirements mentioned above for designing multi-layer molecular rotors.

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 MB24− (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 MB24− (M = Sc, Ti, V, and Cr) clusters correspond to cage structures, and the MB24− (M = Mn, Fe, and Co) clusters have similar distorted four-ring tubes with six boron atoms each. Interestingly, the global minima obtained for the NiB24− 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 MB24− (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) indicates that the local delocalization of MB24− (M = Sc, Ti, V, Cr, and Ni) is weaker than that of MB24− (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 MB24− (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 MB24− (M = Sc, Ti, V, Cr, Mn, Fe, Co, and Ni) shows that MB24− (M = Mn, Fe, and Co) has larger first hyperpolarizability, indicating that MB24− (M = Mn, Fe, and Co) has a strong nonlinear optical response. Hence, MB24− (M = Mn, Fe, and Co) might be considered as a promising nonlinear optical boron-based nanomaterial. The calculated spectra indicate that MB24− (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.

Exploring the structure and electronic properties of germanium doped boron clusters using density functional theory based global optimization method

Density functional theory calculations to investigate the effect of single and double germanium atom doping on the geometric structure and electronic properties of boron clusters with 10 to 20 atoms.

Structural and Electronic Properties of Single-Atom Transition Metal-Doped Boron Clusters MB24 (M = Sc, V, and Mn).

A theoretical study of geometrical structures, electronic properties, and spectral properties of single-atom transition metal-doped boron clusters MB24 (M = Sc, V, and Mn) is performed using the CALYPSO approach for the global minimum search, followed by density functional theory calculations. The global minima obtained for the VB24 and MnB24 clusters correspond to cage structures. Interestingly, the global minima obtained for the ScB24 cluster tend to a three-ring tubular structure. Population analyses and valence electron density analyses reveal that partial electrons on transition-metal atoms transfer to boron atoms. The localized orbital locator of MB24 (M = Sc, V, and Mn) indicates that the electron delocalization of ScB24 is stronger than that of VB24 and MnB24, and there is no obvious covalent bond between doped metals and B atoms. The spin density and spin population analyses reveal that MB24 (M = Sc, V, and Mn) have different spin characteristics which are expected to lead to interesting magnetic properties and potential applications in molecular devices. The calculated spectra indicate that MB24 (M = Sc, V, and Mn) has 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.

Structures, Electronic, and Spectral Properties of Doped Boron Clusters MB12 0/- (M = Li, Na, and K).

Structures and electronic properties of alkali metal atom-doped boron clusters MB120/– (M = Li, Na, K) are determined using the CALYPSO method for the global minimum search followed by density functional theory. It is found that the global minima obtained for the neutral clusters correspond to the half-sandwich structure and those of the monoanionic clusters correspond to the boat-shaped structure. The neutral MB12 (M = Li, Na, K) can be considered as a member of the half-sandwich doped B12 clusters, and the geometrical pattern of anion MB12– (M = Li, Na, K) is a new structure that is different from other doped B12 clusters. Natural population and chemical bonding analyses reveal that the alkali metal atom-doped boron clusters MB12– are characterized as charge transfer complexes, M+B122–, resulting in symmetrically distributed chemical bonds and electrostatic interactions between cationic M+ and boron atoms. The calculated spectra indicate that MB120/– (M = Li, Na, K) has meaningful spectral features that can be compared with future experimental data. Our work enriches the varieties of geometrical structures of doped boron clusters and can provide much insight into boron nanomaterials.

Theoretical investigation on anti-sandwich beryllium-boron clusters Be2mBn (m = 1–3): Fluxionality and multi-aromaticity

Recently, a kind of anti-sandwich beryllium doped boron clusters with middle boron ring and beryllium units on both sides have been experimentally obtained. These clusters attract our interests due to their special three layers geometry with structural fluxionality and aromaticity. Therefore, the anti-sandwich Be2mBn (m = 1–3) clusters are theoretically investigated in this paper to get a general view of this series of beryllium-boron binary systems. The results indicate that all these clusters possessed three-layer geometries after optimization, whose stability is demonstrated by diverse multi-aromaticity. Among these structures, Be4B102+ could be viewed as an analogue of naphthalene depending on its pattern of delocalized π electrons, and Be6B12 has a distorted boron plane and high fluxionality with potential applications as a 3D Wankel motor. Our study further reveals structural rules of anti-sandwich beryllium-boron binary clusters and enriches their structural diversity.

LiB13: A New Member of Tetrahedral-Typed B13 Ligand Half-Surround Cluster

It will get entirely unusual derivatives with gratifying chemical bonding schemes for boron clusters by doping with lithium, the lightest alkalis. The geometric structures and electronic properties of the LiBn0/− (n = 10−20) clusters have been studied through Crystal structure AnaLYsis by Particle Swarm Optimization (CALYPSO) structural search approach along with the density functional theory (DFT) calculations. The low-lying candidates of LiBn0/− (n = 10–20) are reoptimized at the B3LYP functional in conjunction with 6–311 + G(d) basis set. Three forms of geometric configurations are identified for the ground-state structures of LiBn0/− clusters: half-sandwich-type, quasi-planar and drum-type structures. The photoelectron spectra (PES) of the LiBn− clusters have been calculated through time-dependent density functional theory (TD-DFT). A promising LiB13 with tetrahedral-typed B13 ligand half-surround cluster and robust stability is uncovered. The molecular orbital and adaptive natural density partitioning (AdNDP) analysis show that B-B bonds in the B13 moiety combined with the interaction between the B13 shell and Li atom stabilize the C2v LiB13 cluster. Our results advance the fundamental understanding about the alkali metal doped boron clusters.

Structures, stabilities and aromatic properties of endohedrally transition metal doped boron clusters M@B22, M = Sc and Ti: a theoretical study.

A genetic search algorithm in conjunction with density functional theory calculations was used to determine the lowest-energy minima of the pure B22 cluster and thereby to evaluate the capacity of its isomers to form endohedrally doped cages with two transition metal atoms M (M = Sc and Ti). An important charge transfer from metal atoms M to the boron cage takes place, stabilizing the endohedral compounds, as predicted with the genetic algorithm implemented. High-level coupled-cluster theory CCSD(T) calculations were carried out to confirm that the structures found are the lowest-energy isomers. For a deeper understanding of the doping effects and related charge transfer, the best structural motif of the B22 isomers was also determined when the bare cages are in anionic states, such as B222- and B224-. It was found that B22 has an appropriate size, geometric shape and electronic state to host the chosen metal atoms and, consequently, to form stable endohedrally doped compounds Ti@B22 (C2v, 4-Ti) and Sc@B22 (C2v, 5-Sc). The chemical bonding was analyzed in order to understand the molecular orbitals that these novel systems form. The cage aromaticity was evaluated by means of the nuclear independent chemical shift (NICS(0)iso) indices, the isochemical shielding surface (ICSSzz), the anisotropy of the current induced density (ACID) maps, and the magnetic ring current Gauge-Including Magnetically Induced Current (GIMIC) method, indicating that aromaticity plays a crucial role in the stabilization of endohedrally doped boron clusters. Finally, the thermodynamic stability of the latter, using parameters derived from density functional theory (DFT), was evaluated. Ab initio molecular dynamics (AIMD) simulations were performed to elucidate the stability, at high temperature, of the most stable endohedrally doped boron clusters 4-Ti and 5-Sc.

Probing the structure and electronic properties of beryllium doped boron clusters: A planar BeB16\u2212 cluster motif for metallo-borophene

Beryllium-doped boron clusters display essential similarities to borophene (boron sheet) with a molecular structure characterized by remarkable properties, such as anisotropy, metallicity and high conductivity. Here we have determined low-energy structures of BeBn0/− (n = 10–20) clusters by utilizing CALYPSO searching program and DFT optimization. The results indicated that most ground states of clusters prefer plane or quasi-plane structures by doped Be atom. A novel unexpected fascinating planar BeB16− cluster with C2v symmetry is uncovered which possesses robust relative stability. Furthermore, planar BeB16− offers a possibility to construct metallo-borophene nano-materials. Molecular orbital and chemical bonding analysis reveal the peculiarities of BeB16− cluster brings forth the aromaticity and the strong interaction of B-B σ-bonds in boron network.

Structural evolution and electronic properties of medium-sized boron clusters doped with scandium

Doping of boron-based materials with transition metal atoms allows one to tune or modify the properties and structure of the materials. In this work, an extensive search for the global minima on potential energy surfaces of ScBn and ScB clusters has been performed using the CALYPSO method. The structural evolution of scandium doped boron clusters of this range is found to proceed in three steps; namely, the formation of half-sandwich type structures is followed by the formation of drum-like structures with the Sc atom located at the center and terminates with the cage-like structures. It is also found that highly symmetrical geometric structures are more common for the smaller size range of . The neutral ScB13 cluster is identified as magic on the basis of an analysis of relative stabilities in the ScBn series. Our analysis of chemical bonding has shown that the stability of this cluster is mainly due to the formation of several delocalized -bonding molecular orbitals composed of Sc 3d and B 2s atomic orbitals. These bonds appear to be responsible for the enhanced stability of ScB13 with respect to other Sc-doped boron clusters.

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

Because of their interesting structures and bonding and potentials as motifs for new nanomaterials, size-selected boron clusters have received tremendous interest in recent years. In particular, boron cluster anions (Bn-) have allowed systematic joint photoelectron spectroscopy and theoretical studies, revealing predominantly two-dimensional structures. The discovery of the planar B36 cluster with a central hexagonal vacancy provided the first experimental evidence of the viability of 2D borons, giving rise to the concept of borophene. The finding of the B40 cage cluster unveiled the existence of fullerene-like boron clusters (borospherenes). Metal-doping can significantly extend the structural and bonding repertoire of boron clusters. Main-group metals interact with boron through s/p orbitals, resulting in either half-sandwich-type structures or substitutional structures. Transition metals are more versatile in bonding with boron, forming a variety of structures including half-sandwich structures, metal-centered boron rings, and metal-centered boron drums. Transition metal atoms have also been found to be able to be doped into the plane of 2D boron clusters, suggesting the possibility of metalloborophenes. Early studies of di-metal-doped boron clusters focused on gold, revealing ladder-like boron structures with terminal gold atoms. Recent observations of highly symmetric Ta2B6- and Ln2Bn- (n = 7-9) clusters have established a family of inverse sandwich structures with monocyclic boron rings stabilized by two metal atoms. The study of size-selected boron and doped-boron clusters is a burgeoning field of research. Further investigations will continue to reveal more interesting structures and novel chemical bonding, paving the foundation for new boron-based chemical compounds and nanomaterials.

NbB12-: a new member of half-sandwich type doped boron clusters with high stability.

A theoretical study of geometrical structures and electronic properties of niobium-doped boron clusters is performed using the CALYPSO approach for the global minimum search followed by density functional theory calculations. It is found that the global minima obtained for the neutral clusters correspond to the half-sandwich structures at n = 10-17 and the tubular-type structures at n = 18-20. The geometrical patterns in the anion series are more complex. The geometries undergo a transformation from the wheel-like structure of NbB10- to the half-sandwich ones beginning at n = 11 and finally to the drum-shaped structures at n ≥ 18. A fascinating NbB12- cluster is uncovered by our structural search, which shows robust stability and can be considered as a new member of the half-sandwich transition metal doped boron clusters. The chemical bond analysis indicates that the high stability is due to the strong interactions between the Nb atom and the B12 host as well as to the strong B-B covalent bonds. Our study will enrich the database of geometrical structures of transition metal doped boron clusters and will stimulate future synthesis of boron-based nanomaterials.