Gold as hydrogen: structural and electronic properties and chemical bonding in Si3Au3(+0-) and comparisons to Si3H3(+0-)

A single Au atom has been shown to behave like H in its bonding to Si in several mono- and disilicon gold clusters. In the current work, we investigate the AuH analogy in trisilicon gold clusters, Si3Au3(+0-). Photoelectron spectroscopy and density functional calculations are combined to examine the geometric and electronic structure of Si3Au3-. We find that there are three isomers competing for the ground state of Si3Au3- as is the case for Si3H3-. Extensive structural searches show that the potential energy surfaces of the trisilicon gold clusters (Si3Au3-, Si3Au3, and Si3Au3+) are similar to those of the corresponding silicon hydrides. The lowest energy isomers for Si3Au3- and Si3Au3 are structurally similar to a Si3Au four-membered ring serving as a common structural motif. For Si3Au3+, the 2pi aromatic cyclotrisilenylium auride ion, analogous to the aromatic cyclotrisilenylium ion (Si3H3+), is the most stable species. Comparison of the structures and chemical bonding between Si3Au3(+0-) and the corresponding silicon hydrides further extends the isolobal analogy between Au and H.     

Photoelectron spectroscopy and density functional calculations of CuSi(n)- (n = 4-18) clusters

We conducted a combined anion photoelectron spectroscopy and density functional theory study on the structural evolution of copper-doped silicon clusters, CuSi(n)(-) (n = 4-18). Based on the comparison between the experiments and theoretical calculations, CuSi(12)(-) is suggested to be the smallest fully endohedral cluster. The low-lying isomers of CuSi(n)(-) with n ≥ 12 are dominated by endohedral structures, those of CuSi(n)(-) with n < 12 are dominated by exohedral structures. The most stable structure of CuSi(12)(-) is a double-chair endohedral structure with the copper atom sandwiched between two chair-style Si(6) rings or, in another word, encapsulated in a distorted Si(12) hexagonal prism cage. CuSi(14)(-) has an interesting C(3h) symmetry structure, in which the Si(14) cage is composed by three four-membered rings and six five-membered rings.     

Structure evolution of gold cluster anions between the planar and cage structures by isoelectronic substitution: Au(n)- (n = 13-15) and MAu(n)- (n = 12-14; M = Ag, Cu)

The structural and electronic effects of isoelectronic substitution by Ag and Cu atoms on gold cluster anions in the size range between 13 and 15 atoms are studied using a combination of photoelectron spectroscopy and first-principles density functional calculations. The most stable structures of the doped clusters are compared with those of the undoped Au clusters in the same size range. The joint experimental and theoretical study reveals a new C(3v) symmetric isomer for Au(13)(-), which is present in the experiment, but has hitherto not been recognized. The global minima of Au(14)(-) and Au(15)(-) are resolved on the basis of comparison between experiment and newly computed photoelectron spectra that include spin-orbit effects. The coexistence of two isomers for Au(15)(-) is firmly established with convincing experimental evidence and theoretical calculations. The overall effect of the isoelectronic substitution is minor on the structures relative to those of the undoped clusters, except that the dopant atoms tend to lower the symmetries of the doped clusters.     

First principles study on the solvation and structure of C2O4(2-)(H2O)n, n = 6-12

The structures and energies of hydrated oxalate clusters, C2O4(2-)(H2O)n, n = 6-12, are obtained by density functional theory (DFT) calculations and compared to SO4(2-)(H2O)n. Although the evolution of the cluster structure with size is similar to that of SO4(2-)(H2O)n, there are a number of important and distinctive futures in C2O4(2-)(H2O)n, including the separation of the two charges due to the C-C bond in C2O4(2-), the lower symmetry around C2O4(2-), and the torsion along the C-C bond, that affect both the structure and the solvation energy. The solvation dynamics for the isomers of C2O4(2-)(H2O)12 are also examined by DFT based ab initio molecular dynamics.     

High-stability hydrogenated silicon-carbon clusters: a full study of Si2C2H2 in comparison to Si2C2, C2B2H4, and other similar species

The structural and electronic characteristics of the Si2C2H2 and Si2C2 clusters are studied by ab initio calculations based on coupled cluster and density functional theory using the hybrid B3LYP functional. In addition, similar species, such as SiC2H2 and Si3C2H2, are also studied for comparison. It is illustrated that the lowest energy structures of all three hydrogenated clusters, which have the general form Si(n)(CH)2, n = 1, 2, 3, are fully analogous to the structures of the corresponding organometallic isovalent carboranes. The most stable structure of Si2C2H2 is obtained by attaching two hydrogens onto the carbon atoms of a higher energy (+1.5 eV) planar trapezoidal structure of Si2C2, followed by geometry optimization which leads to puckering of the planar structure. Furthermore, it is demonstrated that Si2C2H2 and the other two "similar" hydrogenated clusters are much more stable than the corresponding bare nonhydrogenated clusters. Comparison of Si2C2H2 and C2B2H4 shows that their structural and bonding similarity includes also nuclear rearrangement similarity. The two species are isomerizable with an energy difference between their lowest energy puckered 1,2- and 1,3-isomers of about +/-0.3 eV. It is suggested that SiC2H2, Si2C2H2, and Si3C2H2 are special cases of a larger class of stable clusters. It is speculated on the basis of the calculated infrared spectrum that Si2C2H2 and perhaps other members of this class of clusters could be found in appreciable abundance in interstellar space.     

A new class of silicon-carbon clusters: a full study of the hydrogenated SinC2H2, n=3,4,5, clusters in comparison with their isoelectronic carboranes C2BnHn+2

The structural and electronic characteristics of the Si(n)C(2)H(2), n=3,4,5, clusters are studied by ab initio calculations based on coupled cluster and density functional theory using the hybrid B3LYP functional. It is demonstrated that all three clusters are structurally and electronically homologous to the corresponding isoelectronic organometallic carboranes C(2)B(n)H(n+2). This homology, which is in full agreement with the analogy of Si(6) (2-) and B(6)H(6) (2-) demonstrated recently by the author [J. Chem. Phys. 127, 014314 (2007)], includes not only the ground states but also the lower-lying isomers as well. These lowest lying isomers can be obtained by ortho, para, and meta substitutions from the corresponding Si(n) (2-), n=3,4,5, dianions. The energetic ordering of the low-lying isomers is in full agreement with the known valence and topological charge stability rules developed for carboranes. The hydrogenated clusters are much more stable than their nonhydrogenated counterparts. It is suggested that Si(3)C(2)H(2), Si(4)C(2)H(2), and Si(5)C(2)H(2), which can be probably found in interstellar space, are special examples of a general class of silicon-carbon clusters of the form Si(n)C(2)H(2), with analogous properties and similarities to the corresponding carboranes C(2)B(n)H(n+2). It is furthermore illustrated that the lowest energy structures of the Si(n)C(2) clusters can be obtained through a systematic and straightforward procedure from the Si(n)C(2)H(2) clusters. The present results could hopefully make possible the exploitation of the rich borane and carborane chemistry for the design and development of novel silicon and silicon-carbon composite nanomaterials.     

Structure and bonding in cyclic isomers of B2AlH(n)(m) (n = 3-6, m = -2 to +1): a comparative study with B3H(n)(m), BAl2H(n)(m) and Al3H(n)(m)

The structure, bonding and energetics of B(2)AlH(n)(m) (n = 3-6, m = -2 to +1) are compared with corresponding homocyclic boron, aluminum analogues and BAl(2)H(n)(m) using density functional theory (DFT). Divalent to hexacoordinated boron and aluminum atoms are found in these species. The geometrical and bonding pattern in B(2)AlH(4)(-) is similar to that for B(2)SiH(4). Species with lone pairs on the divalent boron and aluminum atoms are found to be minima on the potential energy surface of B(2)AlH(3)(2-). A dramatic structural diversity is observed in going from B(3)H(n)(m) to B(2)AlH(n)(m), BAl(2)H(n)(m) and Al(3)H(n)(m) and this is attributable to the preference of lower coordination on aluminum, higher coordination on boron and the higher multicenter bonding capability of boron. The most stable structures of B(3)H(6)(+), B(2)AlH(5) and BAl(2)H(4)(-) and the trihydrogen bridged structure of Al(3)H(3)(2-) show an isostructural relationship, indicating the isolobal analogy between trivalent boron and divalent aluminum anion.     

Ultraviolet and infrared photodissociation of Si(+)(C6H6)n and Si(+)(C6H6)(n)Ar clusters

Ion-molecule complexes of the form Si(+)(C6H6)n and Si(+)(C6H6)(n)Ar are produced by laser vaporization in a pulsed nozzle cluster source. These clusters are mass-selected and studied with ultraviolet (355 nm) photodissociation and resonance-enhanced infrared photodissociation spectroscopy in the C-H stretch region of benzene. In the UV, Si(+)(C6H6)n clusters (n = 1-5) fragment to produce the Si(+)(C6H6)n mono-ligand species, suggesting that this ion has enhanced relative stability. IR photodissociation of Si(+)(C6H6)n complexes occurs by the elimination of benzene, while Si(+)(C6H6)(n)Ar complexes lose Ar. Resonances reveal C-H vibrational bands in the 2900-3300 cm(-1) region characteristic of the benzene ligand with shifts caused by the silicon cation bonding. The IR spectra confirm that the major component of the Si(+)(C6H6)n ions studied have the pi-complex structure rather than the isomeric insertion products suggested previously.     

Silicon monohydride clusters Si(n)H (n = 4-10) and their anions: structures, thermochemistry, and electron affinities

The molecular structures, electron affinities, and dissociation energies of the Si(n)H/Si(n)H- (n = 4-10) species have been examined via five hybrid and pure density functional theory (DFT) methods. The basis set used in this work is of double-zeta plus polarization quality with additional diffuse s- and p-type functions, denoted DZP++. The geometries are fully optimized with each DFT method independently. The three different types of neutral-anion energy separations presented in this work are the adiabatic electron affinity (EA(ad)), the vertical electron affinity (EA(vert)), and the vertical detachment energy (VDE). The first Si-H dissociation energies, D(e)(Si(n)H --> Si(n) + H) for neutral Si(n)H and D(e)(Si(n)H- --> Si(n)- + H) for anionic Si(n)H- species, have also been reported. The structures of the ground states of these clusters are traditional H-Si single-bond forms. The ground-state geometries of Si5H, Si6H, Si8H, and Si9H predicted by the DFT methods are different from previous calculations, such as those obtained by Car-Parrinello molecular dynamics and nonorthogonal tight-binding molecular dynamics schemes. The most reliable EA(ad) values obtained at the B3LYP level of theory are 2.59 (Si4H), 2.84 (Si5H), 2.86 (Si6H), 3.19 (Si7H), 3.14 (Si8H), 3.36 (Si9H), and 3.56 (Si10H) eV. The first dissociation energies (Si(n)H --> Si(n) + H) predicted by all of these methods are 2.20-2.29 (Si4H), 2.30-2.83 (Si5H), 2.12-2.41 (Si6H), 1.75-2.03 (Si7H), 2.41-2.72 (Si8H), 1.86-2.11 (Si9H), and 1.92-2.27 (Si10H) eV. For the negatively charged ion clusters (Si(n)H- --> Si(n)- + H), the dissociation energies predicted are 2.56-2.69 (Si4H-), 2.80-3.01 (Si5H-), 2.86-3.06 (Si6H-), 2.80-3.03 (Si7H-), 2.69-2.92 (Si8H-), 2.92-3.18 (Si9H-), and 2.89-3.25 (Si10H-) eV.     

Geometries and electronic properties of the neutral and charged rare earth Yb-doped Si(n) (n = 1-6) clusters: a relativistic density functional investigation

The neutral and charged YbSi(n) (n = 1-6) clusters considering different spin configurations have been systematically investigated by using the relativistic density functional theory with generalized gradient approximation. The total bonding energies, equilibrium geometries, Mulliken populations (MP), Hirshfeld charges (HC), fragmentation energies, and highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gaps are calculated and discussed. The optimized geometries indicate that the most stable YbSi(n) (n = 1-6) clusters keep basically the analogous frameworks as the low-lying Si(n)(+1) clusters, while the charged species deviate from their neutral counterparts, and that the doped Yb tends to occupy the substitutional site of the neutral and charged YbSi(n) isomers. The relative stabilities are investigated in terms of the calculated fragmentation energies, exhibiting enhanced stabilities for the remarkably stable neutral and charged YbSi2 and YbSi5 clusters. Furthermore, the calculated MP and HC values show that the charges of the neutral and charged YbSi(n) clusters transfer from the Yb atom to Si(n) atoms and the Yb atom acts as an electron donor, and that the f orbitals of the Yb atom in the neutral and charged YbSi(n) clusters behave as core without involvement in chemical bonding. The calculated HOMO-LUMO gaps indicate that the YbSi2 and YbSi4+ clusters have stronger chemical stabilities. Comparisons of the Yb-doped Si(n) (n = 1-6) with available theoretical results of transition-metal-doped silicon clusters are made. The growth pattern is investigated also.     

Palladium doping of magic gold cluster Au38(SC2H4Ph)24: formation of Pd2Au36(SC2H4Ph)24 with higher stability than Au38(SC2H4Ph)24

A phenylethanethiolate-protected Pd(2)Au(36)(SC(2)H(4)Ph)(24) cluster, which is a two-Pd atom-doped cluster of the well studied magic gold cluster Au(38)(SC(2)H(4)Ph)(24), was synthesized in high purity and its stability was investigated. The results demonstrate that Pd(2)Au(36)(SC(2)H(4)Ph)(24) is more stable than Au(38)(SC(2)H(4)Ph)(24) against degradation in solution and core etching by thiols.     

Electronic structures and thermochemical properties of the small silicon-doped boron clusters B(n)Si (n=1-7) and their anions

We perform a systematic investigation on small silicon-doped boron clusters B(n)Si (n=1-7) in both neutral and anionic states using density functional (DFT) and coupled-cluster (CCSD(T)) theories. The global minima of these B(n)Si(0/-) clusters are characterized together with their growth mechanisms. The planar structures are dominant for small B(n)Si clusters with n≤5. The B(6)Si molecule represents a geometrical transition with a quasi-planar geometry, and the first 3D global minimum is found for the B(7)Si cluster. The small neutral B(n)Si clusters can be formed by substituting the single boron atom of B(n+1) by silicon. The Si atom prefers the external position of the skeleton and tends to form bonds with its two neighboring B atoms. The larger B(7)Si cluster is constructed by doping Si-atoms on the symmetry axis of the B(n) host, which leads to the bonding of the silicon to the ring boron atoms through a number of hyper-coordination. Calculations of the thermochemical properties of B(n)Si(0/-) clusters, such as binding energies (BE), heats of formation at 0 K (ΔH(f)(0)) and 298 K (ΔH(f)([298])), adiabatic (ADE) and vertical (VDE) detachment energies, and dissociation energies (D(e)), are performed using the high accuracy G4 and complete basis-set extrapolation (CCSD(T)/CBS) approaches. The differences of heats of formation (at 0 K) between the G4 and CBS approaches for the B(n)Si clusters vary in the range of 0.0-4.6 kcal mol(-1). The largest difference between two approaches for ADE values is 0.15 eV. Our theoretical predictions also indicate that the species B(2)Si, B(4)Si, B(3)Si(-) and B(7)Si(-) are systems with enhanced stability, exhibiting each a double (σ and π) aromaticity. B(5)Si(-) and B(6)Si are doubly antiaromatic (σ and π) with lower stability.     

Equilibrium geometries, stabilities, and electronic properties of the bimetallic M2-doped Au(n) (M = Ag, Cu; n = 1-10) clusters: comparison with pure gold clusters

The density functional method with relativistic effective core potential has been employed to investigate systematically the geometrical structures, relative stabilities, growth-pattern behaviors, and electronic properties of small bimetallic M(2)Au(n) (M = Ag, Cu; n = 1-10) and pure gold Au(n) (n ≤ 12) clusters. The optimized geometries reveal that M(2) substituted Au(n+2) clusters and one Au atom capped M(2)Au(n-1) structures are dominant growth patterns of the stable alloyed M(2)Au(n) clusters. The calculated averaged atomic binding energies, fragmentation energies, and the second-order difference of energies as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. The analytic results exhibit that the planar structure Ag(2)Au(4) and Cu(2)Au(2) isomers are the most stable geometries of Ag(2)Au(n) and Cu(2)Au(n) clusters, respectively. In addition, the HOMO-LUMO gaps, charge transfers, chemical hardnesses and polarizabilities have been analyzed and compared further.     

Structures and relative stability of medium- and large-sized silicon clusters. VI. Fullerene cage motifs for low-lying clusters Si(39), Si(40), Si(50), Si(60), Si(70), and Si(80)

We performed a constrained search, combined with density-functional theory optimization, of low-energy geometric structures of silicon clusters Si(39), Si(40), Si(50), Si(60), Si(70), and Si(80). We used fullerene cages as structural motifs to construct initial configurations of endohedral fullerene structures. For Si(39), we examined six endohedral fullerene structures using all six homolog C(34) fullerene isomers as cage motifs. We found that the Si(39) constructed based on the C(34)(C(s):2) cage motif results in a new leading candidate for the lowest-energy structure whose energy is appreciably lower than that of the previously reported leading candidate obtained based on unbiased searches (combined with tight-binding optimization). The C(34)(C(s):2) cage motif also leads to a new candidate for the lowest-energy structure of Si(40) whose energy is notably lower than that of the previously reported leading candidate with outer cage homolog to the C(34)(C(1):1). Low-lying structures of larger silicon clusters Si(50) and Si(60) are also obtained on the basis of preconstructed endohedral fullerene structures. For Si(50), Si(60), and Si(80), the obtained low-energy structures are all notably lower in energy than the lowest-energy silicon structures obtained based on an unbiased search with the empirical Stillinger-Weber potential of silicon. Additionally, we found that the binding energy per atom (or cohesive energy) increases typically >10 meV with addition of every ten Si atoms. This result may be used as an empirical criterion (or the minimal requirement) to identify low-lying silicon clusters with size larger than Si(50).     

Photoelectron imaging and theoretical calculations of gold-silver hydrides: comparing the characteristics of Au, Ag and H in small clusters

Structures and electronic properties of the mixed metal hydride anions AuAgH(-), Au(2)AgH(-), AuAg(2)H(-) and their neutrals are studied using anionic photoelectron imaging and theoretical calculations. The three isomers of AuAgH(-) are determined to be linear and those of AuAgH are determined to have C(s) symmetry. The structures of Au(2)AgH(-), AuAg(2)H(-) and their corresponding neutrals are determined to be planar with C(s) or C(2v) symmetries. The vertical detachment energies (VDEs) and adiabatic detachment energies (ADEs) of these anions are reported. Similar to the homonuclear Au(m)(-) and Ag(n)(-) clusters, the metal hydride anions with an even number of valence electrons have higher VDEs than those with an odd number. Variation of the VDEs of these metal hydride anions with interchange of Au, Ag and H (for example Au(m)Ag(n)(-)→ Au(m-1)Ag(n+ 1)(-), or Au(m-1)Ag(n)H(-)) will be shown to be characterized by the electronegativities of Au, Ag and H. The results presented in this study provide important insights into the similar and different characteristics of these three elements in small clusters.     

Structure and stability of oxygen adsorption on Si(n) (n = 5-10) clusters

The structures, binding energies, and electronic properties of one oxygen atom (O) and two oxygen atoms (2O) adsorption on silicon clusters Si(n) with n ranging from 5 to 10 are studied systematically by ab initio calculations. Twelve stable structures are obtained, two of which are in agreement with those reported in previous literature and the others are new structures that have not been proposed before. Further investigations on the fragmentations of Si(n)O and Si(n)O2 (n = 5-10) clusters indicate that the pathways Si(n)O --> Si(n-1) + SiO and Si(n)O2 --> Si(n-2) + Si2O2 are most favorable from thermodynamic viewpoint. Among the studied silicon oxide clusters, Si8O, Si9O, Si5O2 and Si8O2 correspond to large adsorption energies of silicon clusters with respect to O or 2O, while Si8O, with the smallest dissociation energy, has a tendency to separate into Si7 + SiO. Using the recently developed quasi-atomic minimal-basis-orbital method, we have also calculated the unsaturated valences of the neutral Si(n) clusters. Our calculation results show that the Si atoms which have the largest unsaturated valences are more attractive to O atom. Placing O atom right around the Si atoms with the largest unsaturated valences usually leads to stable structures of the silicon oxide clusters.     

Theoretical study of mixed silicon-lithium clusters Si(n)Li(p)(+) (n=1-6, p=1-2)

Theoretical study on the structures of neutral and singly charged Si(n)Li(p)((+)) (n=1-6, p=1-2) clusters have been carried out in the framework of the density functional theory (DFT) with the B3LYP functional. The structures of the neutral Si(n)Li(p) and cationic Si(n)Li(p)(+) clusters are found to keep the frame of the corresponding Si(n), Li species being adsorbed at the surface. The localization of the lithium cation is not the same one as that of the neutral atom. The Li(+) ion is preferentially located on a Si atom, while the Li atom is preferentially attached at a bridge site. A clear parallelism between the structures of Si(n)Na(p) and those of Si(n)Li(p) appears. The population analysis show that the electronic structure of Si(n)Li(p) can be described as Si(n)(p)(-)+pLi(+) for the small sizes considered. Vertical and adiabatic ionization potentials, adsorption energies, as well as electric dipole moments and static dipolar polarizabilities, are calculated for each considered isomer of neutral species.     

Microsolvation of the dicyanamide anion: [N(CN)(2)(-)](H(2)O)n (n = 0-12)

Photoelectron spectroscopy is combined with ab initio calculations to study the microsolvation of the dicyanamide anion, N(CN)(2)(-). Photoelectron spectra of [N(CN)(2)(-)](H2O)n (n = 0-12) have been measured at room temperature and also at low temperature for n = 0-4. Vibrationally resolved photoelectron spectra are obtained for N(CN)(2)(-), allowing the electron affinity of the N(CN)2 radical to be determined accurately as 4.135 +/- 0.010 eV. The electron binding energies and the spectral width of the hydrated clusters are observed to increase with the number of water molecules. The first five waters are observed to provide significant stabilization to the solute, whereas the stabilization becomes weaker for n > 5. The spectral width, which carries information about the solvent reorganization upon electron detachment in [N(CN)(2)(-)](H2O)n, levels off for n > 6. Theoretical calculations reveal several close-lying isomers for n = 1 and 2 due to the fact that the N(CN)(2)(-) anion possesses three almost equivalent hydration sites. In all the hydrated clusters, the most stable structures consist of a water cluster solvating one end of the N(CN)(2)(-) anion.     

Highly luminescent octanuclear Au(I)-Cu(I) clusters adopting two structural motifs: the effect of aliphatic alkynyl ligands

Reactions of the homoleptic (AuC(2)R)(n) precursors with stoichiometric amount of diphosphine ligand PPh(2)C(6)H(4)PPh(2) (P^P) and Cu(+) ions lead to an assembly of a new family of bimetallic clusters [Au(6)Cu(2)(C(2)R)(6)(P^P)(2)](2+) (type I; R=9-fluorenolyl (1), diphenylmethanolyl (2), 2,6-dimethyl-4-heptanolyl (3), 1-cyclohexanolyl (4), Cy (5), tBu (6)). In the case of R=1-cyclohexanolyl, a structurally different complex [Au(6)Cu(2)(C(2)C(6)H(11)O)(6)(P^P)(3)](2+) (7, type II) could be obtained by treatment of 4 with one equivalent of the diphosphine, while for R=isopropanolyl only the latter type of cluster [Au(6)Cu(2)(C(2)C(3)H(7)O)(6)(P^P)(3)](2+) (8) was detected. Steric bulkiness of the alkynyl ligands and O···H-O hydrogen bonding are suggested to play an important role in stabilizing the type I and type II cluster structural motif, respectively. All the complexes exhibit intense photoluminescence in solution with emission parameters that depending on the geometrical arrangement of the octanuclear metal core. The clusters 1-4 and 6 show single emission band in a blue region (469-488 nm) with maximum quantum yield of 94% (4), while structurally different 7 and 8 emit yellow-orange (590 nm) with unity quantum efficiency. The theoretical DFT calculations of the electronic structures have been carried out to demonstrate that the metal-centered triplet emission within the heterometallic core plays a key role for the observed phosphorescence.