Structures and bonding of auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3): comparison with dihydride polyboroenes

Peng Shao*a, Li-Ping Dinga, Cheng Lub, Jiang-Tao Caia, Bo Liua and Chang-Bo Suna
aCollege of Science, Shaanxi University of Science & Technology, Xi'an 710021, China. E-mail: scu_sp@163.com
bDepartment of Physics, Nanyang Normal University, Nanyang 473061, China

Received 9th August 2015 , Accepted 12th October 2015

First published on 12th October 2015


Abstract

Equilibrium structures of auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3) are obtained from density functional theory-based calculations. Results show that the ground states of Au2B9+, Au2B13+ and Au2B142− can be obtained by adding two Au atoms to the corresponding ground-state pure boron clusters. For Au2B102−, Au2B11, Au2B142− and Au2B15, the ladder structures are proven to be the ground states at TPSS, OVGF, and CCSD(T) levels, which is similar to that of dihydride polyboroenes. AdNDP analysis indicates that the two rows of boron atoms in these auropolyboroenes are bonded by delocalized three-, four-, or five-center σ and π bonds. Especially, the dominant bonding patterns in Au2B11 and Au2B142− bear similarities to those of dihydride polyboroenes. The photoelectron spectroscopy (PES) spectra for anionic clusters were simulated to facilitate the experimental PES spectra. In addition, the fragmentation energies and products against different decay channels are estimated and discussed.


1. Introduction

Over the past decades, pure boron clusters have attracted considerable attention because boron possesses interesting chemical bonding properties and plays essential roles in advancing chemical bonding models.1,2 A series of joint experimental and theoretical studies have presented that all the boron clusters possess planar or quasi-planar geometries in their ground states up to very large sizes.3–16 For the cations, the transitions from planar or quasi-planar to three-dimensional (3D) structures occur at B16+,15 at B20 for the neutrals12 and at least up to B23 for the anions,8–13 where the exact size for the anions is still unclear. Furthermore, boron clusters have been tested to possess rich chemical bonding properties, which can be understood based on the π and σ aromaticity and antiaromaticity. Boron-based materials, especially boron hydrides (boranes), have been popular due to their ability to store a good amount of hydrogen and release it upon chemical treatment when needed.17 To explain the structure of all known boron hydrides, Longuet-Higgins, Lipscomb and co-workers18,19 first put forward the concept of three-center two-electron (3c-2e) bonding. The 3c-2e bond represents a milestone in establishing the validity of the molecular orbital theory. From then on, a surge of studies were devoted to the structure and chemical bonding of boranes. Additionally, the electron delocalization and aromaticity have also been proved to be key bonding features for the electron-deficient boranes.

In previous studies, Au/H similarity is well supported by surprising experimental discovery of H/AuPR3 analogy.20 Further analogy between a bare Au atom and H has also been discovered in gas-phase binary Au clusters.21–23 The analogy extended the structures and bonding in variety of boron-hydrides to B–Au compounds. Such as the heptaboron aurides B7Au2 anion22 was first studied by Zhai et al. using photoelectron spectroscopy and ab initio calculations in 2006. The results showed that B7Au2 possesses an extremely stable planar structure, which is identical to that of B7H2. What's more, the B–Au bonding was shown to be covalent, similar to the B–H bonding, indicating the Au mimics H in its bonding to boron. A subsequent study23 via density functional theory calculations showed that the closo-auroboranes BxAux2− (x = 7–14) dianions possess structure and bonding analogous to the famous deltahedral closo-borane cages, BxHx2−, demonstrating the Au indeed similar to H. In 2010, Zhai groups24 found that the highly covalent B–Au bond in B10Au is also similar to the B–H bond in B10H. The BAun0/− (n = 1–4)25 and B2Aun0/− (n = 1, 3, 5)26 auroboranes have also been proved to possess similar geometrical structures and bonding patterns with corresponding boron hydrides. More recently, Chen et al.27 studied B6Aun0/− (n = 1–3), which provide new examples for the Au/H analogy in Au alloy clusters. These studies on B–Au clusters all indicate that Au serve as monovalent σ ligand apes hydrogen.

The investigations mentioned above unquestionably provide precious information on Au/H analogy, as well as the structure and chemical bonding of boron-gold alloy clusters. However, the previous studies only limited to neutral and anionic B–Au compounds. In the current study, we report a new specie cations with the stable closed structures, such as [Au2(B4)xB]+ in combination with [Au2(B4)xB3] and [Au2(B4)xB2]2− (x = 2, 3) auropolyboroenes. In previous investigation, Li et al.28 have conjectured that [H2(B4)xB3], [H2(B4)xB2]2− and [H2(B4)xB]+ should possess similar characteristics with the stable closed-shell H2B7, H2B102−, H2B9+ species, respectively. Taking Au/H analogy into account, we think the corresponding auropolyboroenes should also show interesting properties. In this case, we performed quantum chemical calculations on the auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2−, [Au2(B4)xB]+ (x = 2, 3) using density functional theory (DFT) method to address several issues. (1) What patterns of structure and chemical bonding exist in them? (2) Whether they are high stable molecules? (3) Whether their bonding pattern bear similarities to those in dihydride polyboroenes especially for cations. Additionally, in order to facilitate the synthesis of auropolyboroenes in the bulk or deposited on surfaces, we have calculated vertical detachment energy (VDE) for the concerned anions and adiabatic ionization potentials (AIP) for neutrals. Meanwhile, the most preferred fragmentation channels and products were also predicted.

2. Computational methods

Using density functional theory (DFT) method TPSS with the τ-dependent gradient-corrected functional,29 as implemented in GAUSSIAN09 codes,30 we performed quantum chemical calculations for auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3). The Stuttgart/Bonn relativistic effective core (SDD)31 was adopted for gold atom, and all-electron 6-311+G* basis set32 with polarization and diffuse functions for boron and hydrogen atoms. The reliability and accuracy of functional form was first guided by extensive tests performed on two-atom clusters (Au2, B2, AuB and AuB), using a variety of methods (the hybrid DFT methods: B3LYP,33 B3PW91;33–35 the DFT methods: TPSS,29 PBE,36 BPW91,34,37 and PW91;33,34 M06 (ref. 38) and CCSD39–41). The tested results combined with the corresponding experimental values and reliable theoretical values are summarized in Table 1. According to the results we found that the bond length and frequency of Au2 and B2 based on the TPSS/Au/SDD/B/6-311+G* level are in good agreement with the experimental values.42–44 The TPSS results for AuB and AuB dimers were also comparable with other DFT functions and ab initio calculations25 which further validate our selected methods.
Table 1 The bond length (r) and vibrational frequency (ω) of the binary clusters B2, Au2, AuB and AuB at different levels
  B3LYP B3PW91 TPSS PBE BPW91 PW91 M06 CCSD Exp./theo.
a Experiment ref. 42–44.b Theoretical studied bond length ref. 25.
Au2 (1gD∞h)
r 2.58 2.55 2.55 2.56 2.56 2.56 2.60 2.57 2.47a
ω 163 169 171 166 166 167 160 172 191a
[thin space (1/6-em)]
B2 (5uD∞h)
r 1.52 1.52 1.53 1.53 1.53 1.53 1.52 1.53 1.59a
ω 1280 1285 1240 1244 1246 1245 1298 1263 1052a
[thin space (1/6-em)]
AuB (1+C∞v)
r 1.95 1.95 1.95 1.94 1.94 1.94 1.98 1.93 1.95b
ω 634 646 655 660 653 660 589 661
[thin space (1/6-em)]
AuB (2+C∞v)
r 2.01 2.00 2.00 1.99 1.99 1.99 2.04 1.97 1.99b
ω 525 543 544 568 558 567 497 574


Our recipe for finding the equilibrium structures of auropolyboroenes is as follows. The geometries of pure B9+, B102−, B11, B13+, B142− and B15 clusters were first optimized referring to various previous reported structures.45–47 To search for the lowest energy structure of auropolyboroenes, a large number of initial structures were obtained by placing the gold atoms on each possible site of corresponding pure boron clusters as well as by substituting two B atoms using Au in Bn+2+/− clusters. Then, the simulated photoelectron spectroscopy spectra for obtained anionic isomers were performed and compared with available experimental spectra.28 In the geometry optimization procedure, different possible spin multiplicities were also taken into account for each of these isomers. All the structure optimizations were carried out without any symmetry constraints. Finally, we performed the vibrational frequency computations to ensure that the optimized geometry corresponds to a local minimum in potential energy. Afterwards, the first two lowest-energy isomers obtained at TPSS theory method were recalculated by high-quality methods such as CCSD(T)48 and the outer valence Green's function (OVGF).49–51 The vertical detachment energy and adiabatic ionization potentials (VDE and AIP) were calculated by the TPSS function. The chemical bonding analysis was also performed using the adaptive natural density partitioning (AdNDP) method52,53 as implemented in the Multiwfn 3.1 program.54 Molecular orbital visualization was done using the Gaussview 5.0.8 program.55

3. Results and discussion

3.1. Geometrical structures

3.1.1. Pure multicharged B9+, B102−, B11, B13+, B142− and B15 clusters. Based on the method mentioned above, we try our best to search for the equilibrium structures of auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3). The first step is to optimize the geometric structures of pure multicharged B9+, B102−, B11, B13+, B142− and B15 clusters. These pure boron clusters have been well studied by many DFT investigations in previous literatures45–47 except for B142− cluster. In this work, although all the possible isomers of bare boron clusters for each cluster size were extensively explored, only the first two lowest energy geometric structures were shown in Fig. 1. The first and second lowest energy isomers are designated by a lowercase letter a, b. It is worth pointing out that the ground state structures for each cluster are in good agreement with the results given by ref. 45–47. For B142− dianion, there is no detailed geometry available. Our obtained ground state geometry of B142− has ladderlike structure with high point symmetry C2h, whereas for the other clusters, the ladderlike structures are less stable than the ground states in energy.
image file: c5ra15940g-f1.tif
Fig. 1 Optimized structures of multicharged B9+, B102−, B11, B13+, B142− and B15 clusters.
3.1.2. Auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3). In this Section, we discuss the structural details of low-lying isomers of auropolyboroenes Au2B9+, Au2B102−, Au2B11, Au2B13+, Au2B142− and Au2B15 which are obtained at TPSS/Au/SDD/B/6-311+G* level. During the geometry optimizations, a large number of optimized isomers for auropolyboroenes are obtained by adding two gold atoms on each possible site of the corresponding pure boron clusters. We only selected the first two low-lying isomers for each type and listed them together with their point symmetries, electron states and some bond lengths in Fig. 2. Results of our final relative energies for the first and second local minima at three levels of theory TPSS, OVGF and CCSD(T) are also summarized in this figure. Some additional higher-energy isomers along with their information are given in Fig. S1–S3 of ESI. To distinguish the different structures, these isomers for each species are followed by a lowercase letter a, b, c, d, … representing their energies from low to high. The calculated vertical detachment energies and adiabatic ionization potentials (VDEs and AIPs) for the ground state structures are gathered in Table 2. In addition, the ground states of the corresponding dihydride polyboroenes H2B9+, H2B102−, H2B11, H2B13+, H2B142− and H2B15 are also obtained at the same theory level and collected in Fig. S4 of ESI as comparison.
image file: c5ra15940g-f2.tif
Fig. 2 Optimized structures of auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3) with their symmetry, electron state, and relative energies at the OVGF/Au/SDD/B/6-311+G* method in curly brackets, those at the CCSD(T)/Au/SDD/B/6-311+G* level in parentheses, and those otherwise at B3LYP/Au/SDD/B/6-311+G*.
Table 2 Adiabatic ionization potentials (AIP) for the ground state of Au2B9+, Au2B13+ and vertical detachment energy (VDE) for the anions Au2B102−, Au2B11, Au2B142−, Au2B15 at TPSS/Au/SDD/B/6-311+G* level
Species
  Au2B9+ Au2B13+ Au2B102− Au2B142− Au2B11 Au2B15
a The energy required to remove one electron from doubly charged Au2B102− and Au2B142− clusters.b The energy required to remove one electron from singly charged Au2B10 and Au2B14.
AIP 7.11 6.26
VDE 0.43a, 3.76b 0.31a, 4.21b 3.82 3.92


As is shown in Fig. 2, our most stable Au2B9+-a structure has Cs symmetry with two Au atoms bind associatively (dimer like) to different binding sites. It lies 0.13, 1.31 and 0.16 eV lower than the next lowest isomer (Au2B9+-b) at TPSS, OVGF and CCSD(T) levels, respectively. The calculated Au–Au bond length (2.72 Å) of Au2B9+-a is longer than the equilibrium Au–Au distance (2.47 Å) in isolated Au2, whereas it is shorter than the Au–Au distance (2.884 Å) in bulk Au. The second low-lying isomer Au2B9+-b has two Au atoms separately connected to the ladder structure of pure boron B9+-b isomer. This structure is similar to the ground state of H2B9+ (Fig. S4) obtained by us at TPSS/6-311+G* level. For the cationic Au2B13+ auropolyboroene, the ground state can be obtained by adding two Au atoms on the apexes of the low-lying isomer of B13+ (ref. 46). In other words, the B13+ moiety retains intact except for minor distortions. It is followed by a ladder structure (Au2B13+-e) with high C2v point symmetry, which is 0.72 eV higher in energy at TPSS level. However, for the four negatively charged auropolyboroenes Au2B102−, Au2B11, Au2B142− and Au2B15, the ladderlike structures were proved to be the ground states at three theory levels TPSS, OVGF, and CCSD(T). Among these ground state structures, there is a very interesting phenomenon. For even numbered boron clusters, the two Au atoms attached terminally to the corresponding pure B102− and B142− dianions in a trans fashion; while for odd numbered ones, the two Au atoms are in the cis position. In addition, the structures of these lowest-energy isomers are nearly identical to the ground state structures of H2B102−, H2B11, H2B142− and H2B15 as shown in Fig. S4. This further provides new examples for Au/H analogy in Au alloy clusters. Thus, comparing with the stable closed-shell H2B102−, H2B9+ species28 and our obtained dihydride polyboroenes at TPSS level, we can conclude that these closed-shell auropolyboroenes considered in present work are all stable. It is worth mentioning that among these four auropolyboroenes, only the ground state of Au2B142− can be obtained by adding two Au atoms to the lowest-energy structure of B142−, while the ground states of Au2B102−, Au2B11 and Au2B15 is obtained by bonding two Au atoms to the second lowest-energy isomer of pure boron clusters. Moreover, for Au2B102−, Au2B11, Au2B142− and Au2B15 clusters, we found that the structures with two Au atoms binding associatively are higher in energy than that with two Au atoms connected to the boron framework separately. This indicates that the Au atoms prefer to bind atomically to negatively charged boron clusters.

3.2. PES spectra, adiabatic ionization potentials and vertical detachment energies

As is well known, well-resolved PES spectra can serve as electronic “fingerprints” of the underlying clusters. To facilitate the experimental PES spectra, the simulated photoelectron spectroscopy spectra for anionic auropolyboroenes are performed and displayed in Fig. S7 of ESI. As can be seen in Fig. S7, the simulated spectra of H2B11 and Au2B11 are compared with Li's28 experimental spectra of H2B11. The numbers of distinct peaks of simulated photoelectron spectra in the low-binding-energy range of 2–4.6 eV and their relative positions overall agree with the experimental spectra. Those increase the confidence in the reliability of the ground-state structures isomers obtained. Unfortunately, only the experimental spectrum of H2B11 has been reported till now.

Furthermore, the adiabatic ionization potential (AIP) for Au2B9+, Au2B13+ and vertical detachment energy (VDE) for the concerned anions Au2B102−, Au2B11, Au2B142− and Au2B15 are calculated. The values are given at TPSS/Au/SDD/B/6-311+G* level, unless mentioned otherwise. The AIP is computed by following formula:

 
AIP = Eoptimized cationEoptimized neutral (1)
where Eoptimized cation and Eoptimized neutral are the total energies of corresponding clusters obtained at their respective ground-state geometries. Our results of adiabatic ionization potential (AIP) for cations Au2B9+ and Au2B13+ are summarized in Table 2. The optimized ground state structures of neutral Au2B9 and Au2B13 are listed in Fig. S5 in ESI. It can be seen that Au2B9+ and Au2B13+ have significantly large ionization potentials (7.11 and 6.26 eV) at TPSS level, respectively. We recalculated the value of AIP for Au2B9+ and Au2B13+ at OVGF and CCSD(T) levels, the obtained results are 5.66 (6.94) eV for Au2B9+ and 5.72 (5.94) eV for Au2B13+, respectively. The value in parentheses is obtained base on CCSD(T).

The vertical detachment energy (VDE) is defined as the energy difference between the anion and its unoptimized neutral counterpart both at the ground state geometry of the anionic cluster. The ith detachment energy (DE) is the energy required to remove the ith electron corresponding to i-valence anion change into (i-1)-valence ion. Then, it might be given by the following definition:

 
DEi = Eni−1Eni (2)

The DEi is vertical if the Eni−1 of corresponding cluster is the total energy of the unoptimized singly charged (neutral) cluster at the equilibrium geometry of the doubly (singly) charged anions. The VDE values for Au2B102−, Au2B11, Au2B142− and Au2B15 clusters are listed in Table 2. From Table 2, we found that the energies required to remove one-electron from doubly charged Au2B102− and Au2B142− are 0.43 and 0.31 eV, respectively. These are much smaller than those (3.76 and 4.21 eV) of removing one electron from singly charged Au2B10 and Au2B14, suggesting that removing one electron from doubly charged anion is easier than from singly charged anion. In other words, the “second extra electron” is only weakly attached on these single charged anions. In addition, we also calculated the energies required to remove two-electrons from doubly charged Au2B102− and Au2B142−. Their respective values are 3.23 and 4.46 eV, which implies that Au2B102− and Au2B142− dianions are more likely to lose one electron. For the singly charged Au2B11 and Au2B15, the calculated VDEs are 3.82 and 3.92 eV, respectively.

3.3. Fragmentation channels

The fragmentation energies of the lowest-energy isomers for Au2B9+, Au2B102−, Au2B11, Au2B13+, Au2B142− and Au2B15, were calculated at TPSS/Au/SDD/B/6-311+G* level. The results were compared with those of the corresponding H2B9+, H2B102−, H2B11, H2B13+, H2B142− and H2B15. The fragmentation energy is defined as the difference between a parent cluster and its daughters.56
 
ΔEnQ = ΔEniQ + XiEnQ, X = Au or H, i = 1 or 2, (3)
where n and Q denote the total number of clusters and different charge, respectively. Here, we mainly investigate two types of fragmentation channels; one is an Au (or H) atom dissociates and the other is that two gold atoms (or H2 molecule) loss. All the fragmentation processes involving the corresponding detail fragmentation products and energies are listed in Table S1 of ESI.

For the sake of brevity, the ΔE of two channels for Au2B9+, Au2B102−, Au2B11, Au2B13+, Au2B142−, Au2B15 and the corresponding dihydride polyboroenes (H2B9+, H2B102−, H2B11, H2B13+, H2B142− and H2B15) as a function of the cluster size are plotted in Fig. 3a and b, respectively. It can be seen that all the auropolyboroenes possess high fragmentation energies (about ≥3 eV) against the two type fragmentation channels. This appears to be comparable with the values of corresponding dihydride polyboroenes at the same theoretical level. The channels Au2BnQ = AuBnQ + Au for auropolyboroenes Au2B9+, Au2B102−, Au2B13+ and Au2B142− have the lower fragmentation energies. However, for both singly charged anions Au2B11 and Au2B15, the fragmentation pathway Au2BnQ = BnQ + Au2 has the lower energies. Thus for n = 11 and 15, the negative charge mainly resides on the B11 and B15 clusters during the dissociation. Additionally, the ΔE exhibits a maximum at n = 14, indicating that Au2B142− is less probable dissociation into B142− and Au2 due to one need to supply relatively high energy 5.40 eV. In the case of dihydride polyboroenes H2B9+, H2B102−, H2B11, H2B13+, H2B142− and H2B15, we found that removing two H atoms of the ground state dihydride polyboroenes H2BnQ to produce an H2 molecule plus a BnQ cluster is proved to have the lowest fragmentation energies. This may be due to the H2 molecule is the most stable form for H atoms. Meanwhile, for all the dihydride boron clusters, the fragmentation energies of the fragmentation channel H2BnQ = HBnQ + H have very close values. Most strikingly, in the case of H2B142−, the fragmentation energies of pathway HB142− + H is higher than that of B142− + H2 by only 0.04 eV (see Table S1).


image file: c5ra15940g-f3.tif
Fig. 3 (a) Dissociation energies for two fragmentation channels of various auropolyboroenes Au2BnQ. (b) Dissociation energies for two fragmentation channels of dihydride polyboroenes H2BnQ.

3.4. Chemical bonding analyses

Next, we performed chemical bonding analyses for the lowest-energy isomers of considered auropolyboroenes using the adaptive natural density partitioning (AdNDP) method.52,53 AdNDP represents the molecular electronic structure in terms of n-center two-electron (nc-2e) bonds, including the familiar lone pairs (1c-2e) and localized 2c-2e bonds or delocalized nc-2e bonds (3 ≤ n ≤ total numbers of atoms in the system). What's more, it is based on the concept that electron pairs are the main elements of the chemical bonds, which has been used successfully to analyze electronic structure and chemical bonding in previous literatures.28,57–59 In order to keep the length of our paper, we only selected one isomer (Au2B9+, Au2B11 and Au2B142−) as example for each type with different charge (cation and singly/doubly charged anion). Results of the analysis are listed in Fig. 4, 5 and 6, respectively, where the blue and purple regions correspond to the different phases of the molecular wave functions for the MOs. Additionally, the calculated values of the bond order for nc-2e are marked by red text.
image file: c5ra15940g-f4.tif
Fig. 4 Chemical bonding analyses for Au2B9+ auropolyboroene using the AdNDP method.

image file: c5ra15940g-f5.tif
Fig. 5 Chemical bonding analyses for Au2B11 auropolyboroene using the AdNDP method.

image file: c5ra15940g-f6.tif
Fig. 6 Chemical bonding analyses for Au2B142− auropolyboroene using the AdNDP method.

According to the AdNDP analysis, the ground state Au2B9+ isomer (see Fig. 4) has ten lone pair (1c-2e) Au atomic 5d-orbitals with occupy number (ON) in the range of 1.81–1.99|e|, The remaining 26 electrons form seven 2c-2e localized B–B σ-bonds (ON = 1.80–1.92|e|), one 3c-2e Au–Au–B σ-bonds (ON = 1.98|e|), three 4c-2e B–B–B–B σ-bonds and two 4c-2e B–B–B–B delocalized over four boron atoms' π-bonds (ON = 1.84–1.88|e|). We found Au2B9+ isomer possesses two 4c-2e π (four π electrons) bonds. However, it can be seen that these 3c-2e and 4c-2e bonds are very weak as revealed by the calculated bond orders, and hence the bonds in this cluster are mainly described as 2c-2e localized B–B σ-bonds and Au d-orbitals. This may be used to explain the reason why the structure of this cluster is nonplanar and different from the ground state of H2B9+. Additionally, this structure may provide some occupation of s- and p-AOs of boron, avoiding the presence of any unoccupied atomic orbital.

For the anionic Au2B11 (see Fig. 5), AdNDP analysis reveals clearly the existence of ten 1c-2e Au 5d lone-pairs (dxy, dxz, dyz, dx2y2 and dz2 on each Au atom), which is similar to the 1c-2e bonding picture of Au2B9+. In addition, there are two 2c-2e Au–B σ-bonds, six 2c-2e B–B localized σ-bonds, two 4c-2e B–B–B–B π-bonds with the adjacent B2 pair in each B4 and one 5c-2e B–B–B–B–B π-bond. Three π-bonds (six π-electrons) in the boron framework is π aromatic according to the Hückel 4n + 2 rule.60 This is similar to the H2B11 presented by Li et al.,28 and three π orbitals should be also compared with those of 1,3,5-hexatriene. However, the elongated shape of Au2B11 is not consistent with π aromaticity. More interestingly, we found that there are five 3c-2e B–B–B and two 5c-2e B–B–B–B–B σ-bonds in this cluster. The latter bonding features give rise to double aromaticity for Au2B11 cluster because each number of π or σ electrons satisfies the Hückel 4n + 2 rule. Thus, this cluster possesses structure stabilized by electron delocalization both in the π or σ framework. By calculating the bond order, we found that three 3c-2e B–B–B σ-bonds in the center of boron framework are much stronger comparing with the others in two sides. Among the three 5c-2e B–B–B–B–B bonds, only one 5c-2e σ-bond has large bond order 0.11, whereas others are nearly zero. Moreover, the bonding patterns of 2c-2e, 3c-2e, 4c-2e and 5c-2e for Au2B11 cluster are very similar to those of H2B11 (see Fig. S6 in ESI). This further supports the Au/H analogy in Au alloy clusters.

As for the Au2B142−, ten 1c-2e Au 5d-orbitals are similar to those observed in both Au2B9+ and Au2B11 clusters. The two rows of boron atoms in the ladder structure are bonded via multicenter σ and π bonds. The bond pattern between the two Au atoms and the boron framework is σ bond. Additionally, eight 2c-2e B–B σ-bonds (ON = 1.82–1.99|e|), six 3c-2e B–B–B σ-bonds (ON = 1.97|e|), and two 4c-2e B–B–B–B σ-bonds (ON = 1.90|e|) are observed in this cluster. It is worth pointing out that the four middle 3c-2e B–B–B σ-bonds are much stronger than the two others with bond order −1.68. Meanwhile, the two 4c-2e B–B–B–B σ-bonds with bond order 0.14 are also significantly strong compared with the two 4c-2e π-bonds. Two 5c-2e B–B–B–B–B π-bonds and two 5c-2e B–B–B–B–B σ-bonds are all very weak according to the calculated bond order. Just like Au2B11, we found the Au2B142− cluster with three delocalized π bonds (six π electrons) is also π aromatic, while Au2B142− is an σ antiaromatic due to it has eighteen σ bonds fulfilling the Hückel 4n rule for antiaromaticity. This may be the reason why Au2B142− gains major stabilization to the ground state. By the comparison between chemical bonding of Au2B142− and H2B142− (Fig. S6), we found that their major bonds are the same especially for the bonds of Au–B and H–B.

In conclusion, the π bonds in auropolyboroenes Au2B9+, Au2B11 and Au2B142− are delocalized only over parts of the boron frameworks. There is on longer a complete ring of 2c-2e B–B σ bonds found in all boron clusters, whereas there are some 3c-2e B–B–B σ bonds exist. Au2B11 and Au2B142− anions can be considered as covalent complexes between π aromatic B11 and B142− core and two monovalent σ-radicals Au, which are connected by two B–Au σ-bonds. Additionally, the bonding patterns of auropolyboroenes bear similarities to those in dihydride polyboroenes which have the similar structures with auropolyboroenes, especially for the Au–B and H–B bondings. This further supports the Au/H analogy in our studied complexes.

4. Conclusions

In summary, we have presented a systematic study on the equilibrium geometries, electronic structure, thermodynamic stability, and chemical bonding analysis of the auropolyboroenes [Au2(B4)xB3], [Au2(B4)xB2]2− and [Au2(B4)xB]+ (x = 2, 3), comparing with the corresponding dihydride polyboroenes. All the results are summarized as follows:

(1) Based on the geometry optimization, the ladderlike structures are proved to be the ground states for Au2B102−, Au2B11, Au2B142− and Au2B15 at TPSS, OVGF and CCSD(T) levels, which are similar to the corresponding dihydride polyboroenes H2B102−, H2B11, H2B142− and H2B15. These auropolyboroenes and dihydride polyboroenes could be potential nanowires.

(2) To facilitate the experimental PES spectra, the simulated photoelectron spectroscopy spectra for anionic clusters were performed. We also calculated the AIPs for cations and VDEs for the concerned anions. The calculated VDEs for the Au2B102− and Au2B142− dianions show that removing one electron from dianions is easier than from singly charged Au2B10 and Au2B14. That is to say, these two dianions are more likely to lose one electron. Meanwhile, the fragmentation channels and products are deeply discussed.

(3) Using the AdNDP method, we take Au2B9+, Au2B11 and Au2B142− as examples to analyze their electronic structure and chemical bonding. The results reveal that the π-bonds are delocalized only over parts of the boron frameworks. Au2B11 and Au2B142− can be regarded as covalent complexes between a π aromatic core (B11 or B142−) and two monovalent σ-radicals Au. Noteworthy is the fact that the chemical bonding patterns of auroboroenes bear similarities to these in dihydride polyboroenes especially for Au–B and H–B bonding.

Acknowledgements

This work was supported by the Shaanxi University of Science & Technology Key Research Grant (No. BJ15-07).

References

  1. N. N. Greenwood and A. Earnshaw, Chemistry of the Elements, Butterworth-Heinemann, Oxford, 2nd edn, 1997 Search PubMed.
  2. F. A. Cotton, G. Wilkinson, C. A. Murrillo and M. Bochmann, Advanced Inorganic Chemistry, John Wiley & Sons, New York, 6th edn, 1999 Search PubMed.
  3. I. Boustani, Int. J. Quantum Chem., 1994, 52, 1081–1111 CrossRef CAS PubMed; I. Boustani, Phys. Rev. B: Solid State, 1997, 55, 16426–16438 CrossRef.
  4. J. I. Aihara, H. Kanno and T. Ishida, J. Am. Chem. Soc., 2005, 127, 13324–13330 CrossRef CAS PubMed; J. O. C. Jimenez-Halla, R. Islas, T. Heine and G. Merino, Angew. Chem., Int. Ed., 2010, 49, 5668–5671 CrossRef PubMed.
  5. J. E. Fowler and J. M. Ugalde, J. Phys. Chem. A, 2000, 104, 397–403 CrossRef CAS.
  6. J. Aihara, J. Phys. Chem. A, 2001, 105, 5486–5489 CrossRef CAS.
  7. B. Kiran, G. G. Kumar, M. T. Nguyen, A. K. Kandalam and P. Jena, Inorg. Chem., 2009, 48, 9965–9967 CrossRef CAS PubMed.
  8. H. J. Zhai, B. Kiran, J. Li and L. S. Wang, Nat. Mater., 2003, 2, 827–833 CrossRef CAS PubMed.
  9. H. J. Zhai, A. N. Alexandrova, K. A. Birch, A. I. Boldyrev and L. S. Wang, Angew. Chem., Int. Ed., 2003, 42, 6004–6008 CrossRef CAS PubMed; A. N. Alexandrova, H. J. Zhai, L. S. Wang and A. I. Boldyrev, Inorg. Chem., 2004, 43, 3552–3554 CrossRef PubMed.
  10. A. P. Sergeeva, B. B. Averkiev, H. J. Zhai, A. I. Boldyrev and L. S. Wang, J. Chem. Phys., 2011, 134, 224304–224315 CrossRef PubMed.
  11. A. P. Sergeeva, D. Yu. Zubarev, H. J. Zhai, A. I. Boldyrev and L. S. Wang, J. Am. Chem. Soc., 2008, 130, 7244–7246 CrossRef CAS PubMed; W. Huang, A. P. Sergeeva, H. J. Zhai, B. B. Averkiev, L. S. Wang and A. I. Boldyrev, Nat. Chem., 2010, 2, 202–206 CrossRef PubMed; Z. A. Piazza, W. L. Li, C. Romanescu, A. P. Sergeeva, L. S. Wang and A. I. Boldyrev, J. Chem. Phys., 2012, 136, 104310–104318 CrossRef PubMed; A. P. Sergeeva, Z. A. Piazza, C. Romanescu, W. L. Li, A. I. Boldyrev and L. S. Wang, J. Am. Chem. Soc., 2012, 134, 18065–18073 CrossRef PubMed.
  12. B. Kiran, S. Bulusu, H. J. Zhai, S. Yoo, X. C. Zeng and L. S. Wang, Proc. Natl. Acad. Sci. U. S. A., 2005, 102, 961–964 CrossRef CAS PubMed.
  13. A. N. Alexandrova, A. I. Boldyrev, H. J. Zhai and L. S. Wang, Coord. Chem. Rev., 2006, 250, 2811–2866 CrossRef CAS PubMed.
  14. D. Yu. Zubarev and A. I. Boldyrev, J. Comput. Chem., 2007, 28, 251–268 CrossRef CAS PubMed.
  15. E. Oger, N. R. M. Crawford, R. Kelting, P. Weis, M. M. Kappes and R. Ahlrichs, Angew. Chem., Int. Ed., 2007, 46, 8503–8506 CrossRef CAS PubMed.
  16. L. Hanley, J. L. Whitten and S. L. Anderson, J. Phys. Chem., 1988, 92, 5803–5814 CrossRef CAS; P. A. Hintz, S. A. Ruatta and S. L. Anderson, J. Chem. Phys., 1990, 92, 292–303 CrossRef PubMed; S. A. Ruatta, P. A. Hintz and S. L. Anderson, J. Chem. Phys., 1991, 94, 2833–2847 CrossRef PubMed; M. B. Sowa-Resat, J. Smolanoff, A. Lapiki and S. L. Anderson, J. Chem. Phys., 1997, 106, 9511–9522 CrossRef PubMed.
  17. W. N. Lipscomb, Boron Hydrides, Benjamin, New York, 1963 CAS; W. N. Lipscomb, Science, 1977, 196, 1047–1055 CAS.
  18. R. P. Bell and H. C. Longuet-Higgins, Nature, 1945, 155, 328–329 CrossRef CAS PubMed.
  19. W. H. Eberhardt, B. Crawford and W. N. Lipscomb, J. Chem. Phys., 1954, 22, 989–1001 CrossRef CAS PubMed.
  20. (a) K. P. Hall and D. M. P. Mingos, Prog. Inorg. Chem., 1984, 32, 237–254 CrossRef CAS PubMed; (b) J. K. Burdett, O. Eisentein and W. B. Schweizer, Inorg. Chem., 1994, 33, 3261–3268 CrossRef CAS.
  21. K. P. Hall and D. M. P. Mingos, Prog. Inorg. Chem., 1984, 32, 237–254 CrossRef CAS PubMed.
  22. H. J. Zhai, L. S. Wang, D. Yu. Zubarev and A. I. Boldyrev, J. Phys. Chem. A, 2006, 110, 1689–1693 CrossRef CAS PubMed.
  23. D. Y. Zu, J. Li, L.-S. Wang and A. I. Boldyev, Inorg. Chem., 2006, 45, 5269–5271 CrossRef PubMed.
  24. H. J. Zhai, C. Q. Miao, S. D. Li and L. S. Wang, J. Phys. Chem. A, 2010, 114, 12155–12161 CrossRef CAS PubMed.
  25. D. Z. Li and S. D. Li, Int. J. Quantum Chem., 2011, 111, 4418–4424 CrossRef CAS PubMed.
  26. W. Z. Yao, D. Z. Li and S. D. Li, J. Comput. Chem., 2011, 32, 218–225 CrossRef CAS PubMed.
  27. Q. Chen, H. J. Zhai, S. D. Li and L. S. Wang, J. Chem. Phys., 2013, 138, 084306–084313 CrossRef PubMed.
  28. W. L. Li, C. Romanescu, T. Jian and L. S. Wang, J. Am. Chem. Soc., 2012, 134, 13228–13231 CrossRef CAS PubMed.
  29. J. Tao, J. P. Perdew, V. N. Staroverov and G. E. Scuseria, Phys. Rev. Lett., 2003, 91, 146401–146405 CrossRef.
  30. M. J. Frisch, G. W. Trucks and H. B. Schlegel, et al., Gaussian 09 (Revision C.0), Gaussian, Inc, Wallingford, CT, 2009 Search PubMed.
  31. P. Schwerdtfeger, M. Dolg, W. H. E. Schwarz, G. A. Bowmaker and P. D. W. Boyd, J. Chem. Phys., 1989, 91, 1762–1774 CrossRef CAS PubMed.
  32. R. Krishnan, J. S. Binkley, R. Seeger and J. A. Pople, J. Chem. Phys., 1980, 72, 650–654 CrossRef CAS PubMed.
  33. A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS PubMed.
  34. J. P. Perdew and Y. Wang, Phys. Rev. B: Condens. Matter, 1992, 45, 13244–13249 CrossRef.
  35. J. P. Perdew, P. Ziesche and H. Eschrig, Electronic Structure of Solids, Akademie Verlag, Berlin, 1991 Search PubMed.
  36. J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865–3868 CrossRef CAS.
  37. A. D. Becke, Phys. Rev. A, 1988, 38, 3098–3100 CrossRef CAS.
  38. Y. Zhao and D. G. Truhlar, J. Phys. Chem. A, 2006, 110, 13126–13130 CrossRef CAS PubMed.
  39. G. D. Purvis and R. J. Bartlett, J. Chem. Phys., 1982, 76, 1910–1918 CrossRef CAS PubMed.
  40. G. E. Scuseria, C. L. Janssen and H. F. Schaefer III, J. Chem. Phys., 1988, 89, 7382–7387 CrossRef CAS PubMed.
  41. G. E. Scuseria and H. F. Schaefer III, J. Chem. Phys., 1989, 90, 3700–3703 CrossRef CAS PubMed.
  42. K. P. Huber and G. Herzberg, Constants of Diatomic Molecules, Van Nostrand Reinhold, New York, 1979 Search PubMed.
  43. D. R. Lide, CRC Handbook of chemistry and physics, CRC Press, New York, 87th edn, 2006 Search PubMed.
  44. B. Rosen, Spectroscopie Data Relative to Diatomic Molecules, Oxford, Pergamon, 1970 Search PubMed.
  45. N. Akman, M. Tas, C. Oözdoğan and I. Boustani, Phys. Rev. B: Condens. Matter Mater. Phys., 2011, 84, 075463–075473 CrossRef.
  46. T. B. Tai, N. M. Tam and M. T. Nguyen, Theor. Chem. Acc., 2012, 131, 1241–1256 CrossRef.
  47. H. J. Zhai, B. Kiran, J. Li and L.-S. Wang, Nature, 2003, 2, 827–833 CrossRef CAS PubMed.
  48. R. J. Bartlett and M. Musial, Rev. Mod. Phys., 2007, 79, 291–352 CrossRef CAS.
  49. L. S. Cederbaum, J. Phys. B: At. Mol. Phys., 1975, 8, 290–303 CrossRef CAS PubMed.
  50. J. V. Ortiz, J. Chem. Phys., 1996, 104, 7599–7605 CrossRef CAS PubMed.
  51. V. G. Zakrzewski and W. von Niessen, J. Comput. Chem., 1993, 14, 13–18 CrossRef CAS PubMed.
  52. D. Y. Zubarev and A. I. Boldyrev, Phys. Chem. Chem. Phys., 2008, 10, 5207–5217 RSC.
  53. D. Y. Zubarev and A. I. Boldyrev, J. Org. Chem., 2008, 73, 9251–9258 CrossRef CAS PubMed.
  54. T. Lu and F. W. Chen, J. Comput. Chem., 2012, 33, 580–592 CrossRef CAS PubMed.
  55. R. Dennington, T. Keith and J. Millam, GaussView, version 5.0.8, Semichem, Inc, Shawnee Mission, KS, 2007 Search PubMed.
  56. B. K. Rao, P. Jena, M. Manninen and R. M. Nieminen, Phys. Rev. Lett., 1987, 58, 1188–1191 CrossRef CAS.
  57. C. Romanescu, T. R. Galeev, W. L. Li, A. Boldyrev and L. S. Wang, Acc. Chem. Res., 2013, 46, 350–358 CrossRef CAS PubMed.
  58. J. K. Olson and A. I. Boldyrev, J. Phys. Chem. A, 2013, 117, 1614–1620 CrossRef CAS PubMed.
  59. H. Bai, H. J. Zhai, S. D. Li and L. S. Wang, Phys. Chem. Chem. Phys., 2013, 15, 9646–9653 RSC.
  60. M. J. Goldstein, J. Am. Chem. Soc., 1967, 89, 6357–6359 CrossRef CAS.

Footnote

Electronic supplementary information (ESI) available. See DOI: 10.1039/c5ra15940g

This journal is © The Royal Society of Chemistry 2015
Click here to see how this site uses Cookies. View our privacy policy here.