Structures, Stabilities, and Spectral Properties of Endohedral Borospherenes M@B40 0/- (M = H2, HF, and H2O)

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.

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.

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.

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 MB12 0/- (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+B12 2-, resulting in symmetrically distributed chemical bonds and electrostatic interactions between cationic M+ and boron atoms. The calculated spectra indicate that MB12 0/- (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.

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.