Structural Evolution and Electronic Properties of Selenium-Doped Boron Clusters SeB_n^0/- (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.