Structure determination of neutral MgO clusters--hexagonal nanotubes and cages

Structural information for neutral magnesium oxide clusters has been obtained by a comparison of their experimental vibrational spectra with predictions from theory. (MgO)(n) clusters with n = 3-16 have been studied in the gas phase with a tunable IR-UV two-color ionization scheme and size-selective infrared spectra have been measured. These IR spectra are compared to the calculated spectra of the global minimum structures predicted by a hybrid ab initio genetic algorithm. The comparison shows clear evidence that clusters of the composition (MgO)(3k) (k = 1-5) form hexagonal tubes, which confirm previous theoretical predictions. For the intermediate sizes (n≠ 3k) cage-like structures containing hexagonal (MgO)(3) rings are identified. Except for the cubic (MgO)(4) no evidence for bulk like structures is found.     

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.     

Cage and tube structures of medium-sized zinc oxide clusters (ZnO)n (n=24, 28, 36, and 48)

Following our recent work which revealed that the lowest-energy structures of (ZnO)n (n=9-18) follow cage and tube structural growth patterns with stacks of small subunits of (ZnO)2 and (ZnO)3 [Wang et al., J. Phys. Chem. C 111, 4956 (2007)], we have extended the search for the most stable structures to some larger clusters, i.e., (ZnO)n (n=24, 28, 36, and 48) by using gradient-corrected density-functional theory (DFT). A number of starting configurations belonging to different structural motifs were generated from handmade constructions with chemical intuition and then optimized via DFT calculations. Within the size range studied, cage and tube structures were found to be the most preferred structural motifs for the (ZnO)n clusters.     

Stable Tin (n = 2-15) clusters and their geometries: DFT calculations

We present a detailed structural analysis for small Tin (n = 2-15) clusters based on ab initio quantum mechanical calculations of their binding energies, frontier orbital gaps, and second energy derivatives. Local density approximation calculations revealed that while the smaller clusters (n < or = 8) prefer hexagonal atomic arrays with bulklike crystal symmetry, the bigger clusters prefer pentagonal atomic arrays. From the stability criteria of the magic number clusters we could identify three magic number clusters Ti7, Ti13, and Ti15. While the most stable configuration of Ti7 is a decahedral bipyramid induced by tetrahedral atomic arrays, the most stable configuration of Ti13 is an icosahedron. The other stable cluster Ti15 takes a closed-shell icosahedron-like configuration with both pentagonal and hexagonal symmetries. The stability of the Tin clusters strongly depends on their geometries and charge states. The HOMO-LUMO gap of the Tin clusters approaches its bulk value for n > 8. While there is not much difference between the HOMO and LUMO isosurface charge distributions for the Ti7 and Ti13 clusters in their most stable configurations, they are very different in the case of Ti15. Such a distinct charge distribution in Ti15 indicates its singular chemical selectivity over the other two magic number clusters.     

Structure of molybdenum and tungsten sulfide M(x)S(y)+ clusters: experiment and DFT calculations

A combination of experiment and density functional theory was used to investigate the energetics of CO adsorption onto several small M(x)S(y)(+) (M = Mo, W; x/y = 2/6, 3/7, 5/7, 6/8) clusters as a probe of their atomic and electronic structure. Experimentally, tandem mass spectrometry was used to measure the relative yields of M(x)S(y)(+)(CO)(n) cluster adducts formed by collisions between a beam of mass-selected M(x)S(y)(+) cluster ions and CO molecules in a high-pressure collision cell (hexapole ion guide). The most probable M(x)S(y)(+)(CO)(n) adducts observed are those with n < or = x, that is, only one CO molecule bound to each metal site. The notable exception is the M(5)S(7)(+) cluster, for which the n = 6 adduct is found to have nearly the same intensity as the n = x = 5 adduct. Density functional calculations were used to search for the lowest energy structures of the bare M(x)S(y)(+) clusters and to obtain their relative stability for sequential CO binding. The calculated trends in CO binding energies were then compared to the experimental adduct distributions for assigning the ground-state structures. In this way, it was possible to distinguish between two nearly isoenergetic ground-state isomers for the M(2)S(6)(+) and M(3)S(7)(+) clusters, as only one isomer gave a calculated CO stabilization energy trend that was consistent with the experimental data. Similar comparisons of predicted and observed CO adsorption trends also provide evidence for assigning the ground-state structures of the M(5)S(7)(+) and M(6)S(8)(+) clusters. The latter contain metallic cores with most of the sulfur atoms bonded along the edges or in the faces of the metal core structure. The n = 6 and 7 adducts of M(5)S(7)(+) are predicted to be more stable than the n = x = 5 adduct, but only the n = 6 adduct is observed experimentally. The DFT calculations show that the n = 7 adduct undergoes internal bond breaking whereas the n = 6 framework is stable, albeit highly distorted. For the M(6)S(8)(+) cluster, the calculations predict that the two lowest energy isomers can bind more than six CO molecules without fragmentation; however, the apparent binding energy drops significantly for adducts with n > 6. In general, the ability of these small M(x)S(y)(+) clusters to bind more CO molecules than the number of metal atoms is a balance between the gain in CO adsorption energy versus the strain introduced into the cluster structure caused by CO crowding, the consequences of which can be fragmentation of the M(x)S(y)(+)(CO)(n) cluster adduct (n > x).     

Geometric and electronic structures of (BeO)(N) (N = 2-12, 16, 20, and 24): rings, double rings, and cages

The structure of (BeO)(N) clusters (N = 2-12, 16, 20, and 24) are investigated using the method combining the genetic algorithm with density function theory. Benchmark calculation indicates that THSSh functional is reliable to predict the structures of (BeO)(N) cluster. The global minimum structures of (BeO)(N) clusters are rings up to N = 5, double rings at N = 6 and 7 and cages at N ≥ 8. Besides, almost all of the structures of (BeO)(N) cluster are aromatic according to the NICS criterion. Adaptive natural density partitioning analysis reveals that C(6), (BN)(3), and (BeO)(3) rings (C(24) and (BeO)(12) fullerenes) are similar in bonding patterns. The building-up principle of (BeO)(N) is different from that of covalent (BN) and ionic (LiF and MgO) clusters.     

Structure, stability, and cluster-cage interactions in nitride clusterfullerenes M3N@C2n (M = Sc, Y; 2n = 68-98): a density functional theory study

Extensive semiempirical calculations of the hexaanions of IPR (isolated pentagon rule) and non-IPR isomers of C(68)-C(88) and IPR isomers of C(90)-C(98) followed by DFT calculations of the lowest energy structures were performed to find the carbon cages that can provide the most stable isomers of M(3)N@C(2n) clusterfullerenes (M = Sc, Y) with Y as a model for rare earth ions. DFT calculations of isomers of M(3)N@C(2n) (M = Sc, Y; 2n = 68-98) based on the most stable C(2n)(6-) cages were also performed. The lowest energy isomers found by this methodology for Sc(3)N@C(68), Sc(3)N@C(78), Sc(3)N@C(80), Y(3)N@C(78), Y(3)N@C(80), Y(3)N@C(84), Y(3)N@C(86), and Y(3)N@C(88) are those that have been shown to exist by single-crystal X-ray studies as Sc(3)N@C(2n) (2n = 68, 78, 80), Dy(3)N@C(80), and Tb(3)N@C(2n) (2n = 80, 84, 86, 88) clusterfullerenes. Reassignment of the carbon cage of Sc(2)@C(76) to the non-IPR Cs: 17490 isomer is also proposed. The stability of nitride clusterfullerenes was found to correlate well with the stability of the empty 6-fold charged cages. However, the dimensions of the cage in terms of its ability to encapsulate M(3)N clusters were also found to be an important factor, especially for the medium size cages and the large Y(3)N cluster. In some cases the most stable structures are based on the different cage isomers for Sc(3)N and Y(3)N clusters. Up to the cage size of C(84), non-IPR isomers of C(2n)(6-) and M(3)N@C(2n) were found to compete with or to be even more stable than IPR isomers. However, the number of adjacent pentagon pairs in the most stable non-IPR isomers decreases as cage size increases: the most stable M(3)N@C(2n) isomers have three such pairs for 2n = 68-72, two pairs for n = 74-80, and only one pair for n = 82, 84. For C(86) and C(88) the lowest energy IPR isomers are much more stable than any non-IPR isomer. The trends in the stability of the fullerene isomers and the cluster-cage binding energies are discussed, and general rules for stability of clusterfullerenes are established. Finally, the high yield of M(3)N@C(80) (Ih) clusterfullerenes for any metal is explained by the exceptional stability of the C(80)(6-) (Ih: 31924) cage, rationalized by the optimum distribution of the pentagons leading to the minimization of the steric strain, and structural similarities of C(80) (Ih: 31924) with the lowest energy non-IPR isomers of C(760(6-), C(78)(6-), C(82)(6-), and C(84)(6-) pointed out.     

Structure and bonding in coinage metal halide clusters M(n)X(n), M = Cu, Ag, Au; X = Br, I; n = 1-6

Ab initio calculations in the framework of density functional theory (DFT) were performed to study the lowest-energy isomers of noble metal halide clusters M(n)Br(n) and M(n)I(n), for M = Cu, Ag, or Au and n = 1-6. For all species, the most stable structures were found to be cyclic arrangements. Calculated bond lengths and infrared frequencies were compared with the available experimental data. The nature of the ionocovalent bonding was characterized. The stability and fragmentation were also investigated. The present work confirms previous observations on the particular stability of the trimer.     

Experimental and theoretical studies on carbon-nitrogen clusters C2nN7-

C(2n)N7(-) cluster ions are produced by laser ablating on the K(3)[Fe(CN)6] sample. DFT calculations have been performed for these cluster anions. Various isomeric structures of these clusters are optimized and their energies are compared to find the most stable isomers. The most stable structure for C8N7(-) is similar to that of adenine by theoretical calculation, which is in agreement with the collision-induced dissociation (CID) experimental results. With the increasing even numbers of C atoms from 8 to 16, the N atoms in the double-ring structure are gradually substituted by C atoms from the six-membered ring to the five-membered ring. All these C(2n)N7(-) (n = 3-9) clusters exhibit planar aromatic characters. The energy difference and incremental binding energy analyses show that C(2n)N7(-) (n = 4-8) clusters are more stable than C6N7(-) and C18N7(-), which are consistent with the observed mass spectrum.     

CnAl (n=2-11)团簇的结构特征与稳定性

采用密度泛函理论(DFT)的B3LYP方法, 在6-311++G**水平上对CnAl(n=2-11)团簇的几何构型和电子结构进行了结构优化和振动频率计算. 结果表明, n=2的CnAl团簇基态结构为Al原子与两个C原子相连形成的环状结构, n=3-11均为Al原子端基配位的线状结构. 通过对基态结构的能量分析, 得到了n为偶数的CnAl团簇比n为奇数团族稳定的结论.     

CnAl2 (n=1-10)团簇的结构特征与稳定性

采用密度泛函理论(DFT)的B3LYP方法, 在6-311G*水平上对CnAl2 (n=1-10)团簇的几何构型和电子结构进行了结构优化和振动频率计算. 结果表明, 富铝的CAl2团簇基态结构为折线型面状结构, 多碳和双聚体的CnAl2团簇基态结构均为Al原子端基配位的线状结构. 通过对基态结构的能量分析, 得到了n为偶数的CnAl2团簇比n为奇数稳定的结论.     

Structures and dynamics of protonated ammonia clusters

The structures and infrared spectra of protonated ammonia clusters NH(4+)(NH3)n, for n < or = 8, are investigated using density functional-theory (DFT) calculations and semiempirical DFT/molecular dynamics simulations. For n < 5 the clusters are found to be mostly stable up to 100 K, while the larger clusters (n > or = 5) isomerize. Temperature effects are taken into account by performing ab initio molecular dynamics simulations with the computationally tractable self-consistent charges density functional tight-binding method. The infrared spectra at 10 K for the most stable isomers for n = 3-8 compare qualitatively with predissociation experiments, and using a common scaling factor almost quantitative agreement is found. For n > or = 6 the notion of multiple isomers present under the experimental conditions is supported. Of the 13 stable structures for n = 8 only three are found to survive at 100 K. All other clusters isomerize. Cluster structures are inferred from the analysis of the cumulative radial distribution function of the ammonia molecules surrounding the NH(4+) core. The infrared spectra are found to be typical for the structure of the clusters, which should help to relate the experimentally measured infrared spectra to the number and identity of the contributing isomers. For clusters that reorganize to a more stable isomer during the dynamics, the infrared spectrum is generally similar to that of the stable isomer itself. The clusters are found to preferably form globular structures, although chain-like arrangements are also among the low-energy configurations.     

Aun(n=2-9)团簇与乙醇分子相互作用的第一性原理研究

利用密度泛函理论研究了Au_n(n=2-9)团簇吸附一个乙醇分子的结构和电子性质. 研究结果表明: Au_n(n=2-9)团簇的最稳定构型为二维平面结构, Au6团簇最稳定; 吸附过程是通过金团簇上一个特定的金原子与乙醇分子中氧原子相互作用完成, 形成了20种稳定构型; 金原子的配位数对吸附作用影响明显; 作为吸附主体的金团簇和被吸附的乙醇分子在吸附前后构型无明显变化, 它们之间为弱相互作用.     

SinP+m(n+m=5)结构和稳定性的理论研究

在实验的基础上 ,利用量子化学方法对 Sin P+ m( n+m=5 )的各种可能构型进行几何构型优化 ,预测各团簇的稳定结构 ,从中得出各个团簇稳定构型之间的基本关系 ,当 n>m时 ,团簇的稳定构型与 Si+ n 相似 ,而当 n相似文献   

B2Cn+(n=1～9)团簇的结构及其稳定性

采用密度泛函理论(DFT)的B3LYP泛函, 在6-311G*水平上对B2Cn+(n=1～9)团簇的几何构型和电子结构进行了优化和振动频率计算. 结果表明, 在B2Cn+(n=1～9)团簇的基态构型中, B2C2+、B2C3+为具有D∞h对称性的线形结构, B2C7+为具有Cs对称性的立体环状结构, 其余均为平面构型; 其成键顺序为C—C成键优于B—C 成键, B—C成键优于B—B成键. 进一步得到了B2Cn+(n=1～9)团簇的总能量(ET)、零点能(EZ)、摩尔热容(Cp)、标准熵(S0)以及原子化能(ΔEn+). 其结果显示, 随着n的递增, ET、EZ、Cp、S0和ΔEn+数值均呈现增大趋势, 其中EZ数值呈现近似等梯度的增加趋势. 通过对B2Cn+(n=1～9)团簇基态结构的垂直电子亲合势的研究发现, n为奇数的B2Cn+团簇比n为偶数的稳定.     

On the structure and chemical bonding of tri-tungsten oxide clusters W3On- and W3On (n=7-10): W3O8 as a potential molecular model for O-deficient defect sites in tungsten oxides

Electronic and structural properties of a series of tri-tungsten oxide clusters, W3On- and W3On (n=7-10), are investigated using photoelectron spectroscopy and density functional theory (DFT) calculations. Both W 5d and O 2p detachment features are observed for n=7-9, whereas only detachment features from O 2p-type orbitals are observed for W3O10- at high electron binding energies (>7 eV). A large energy gap (approximately 3.4 eV) is observed for the stoichiometric W3O9 cluster, which already reaches the bulk value, suggesting that W3O9 can be viewed as the smallest molecular model for bulk WO3. DFT calculations are carried out to locate the most stable structures for both the anion and neutral clusters; time-dependent DFT method is used to predict the vertical detachment energies and to compare with the experimental data. It is shown that W3O9 possesses a D3h structure, in which each W atom is tetrahedrally coordinated with two bridging O atoms and two terminal O atoms. W3O8 and W3O7 can be viewed as removing one and two terminal O atoms from W3O9, respectively, whereas W3O1) can be viewed as replacing a terminal O in W3O9 by a peroxo O2 unit. We show that W3O8 contains a localized W4+ site, which can readily react with O2 to form the W3O10 clusters with a calculated O2 adsorption energy of -78 kcal/mol. It is suggested that the W3O8 cluster can be viewed as a molecular model for O-deficient site in tungsten oxides.     

Unique CO chemisorption properties of gold hexamer: Au6(CO)n- (n = 0-3)

Elucidating the chemisorption properties of CO on gold clusters is essential to understanding the catalytic mechanisms of gold nanoparticles. Gold hexamer Au(6) is a highly stable cluster, known to possess a D(3)(h) triangular ground state structure with an extremely large HOMO-LUMO gap. Here we report a photoelectron spectroscopy (PES) and quasi-relativistic density functional theory (DFT) study of Au(6)-CO complexes, Au(6)(CO)(n)(-) and Au(6)(CO)(n) (n = 0-3). CO chemisorption on Au(6) is observed to be highly unusual. While the electron donor capability of CO is known to decrease the electron binding energies of Au(m)(CO)(n)(-) complexes, CO chemisorption on Au(6) is observed to have very little effect on the electron binding energies of the first PES band of Au(6)(CO)(n)(-) (n = 1-3). Extensive DFT calculations show that the first three CO successively chemisorb to the three apex sites of the D(3)(h) Au(6). It is shown that the LUMO of the Au(6)-CO complexes is located in the inner triangle. Thus CO chemisorption on the apex sites (outer triangle) has little effect on this orbital, resulting in the roughly constant electron binding energies for the first PES band in Au(6)(CO)(n)(-) (n = 0-3). Detailed molecular orbital analyses lead to decisive information about chemisorption interactions between CO and a model Au cluster.     

Theoretical and experimental analysis of ammonia ionic clusters produced by 252Cf fragment impact on an NH3 ice target

Positive and negatively charged ammonia clusters produced by the impact of (252)Cf fission fragments (FF) on an NH(3) ice target have been examined theoretical and experimentally. The ammonia clusters generated by (252)Cf FF show an exponential dependence of the cluster population on its mass, and the desorption yields for the positive (NH(3))(n)NH(4)(+) clusters are 1 order of magnitude higher than those for the negative (NH(3))(n)NH(2)(-) clusters. The experimental population analysis of (NH(3))(n)NH(4)(+) (n = 0-18) and (NH(3))(n)NH(2)(-) (n = 0-8) cluster series show a special stability at n = 4 and 16 and n = 2, 4, and 6, respectively. DFT/B3LYP calculations of the (NH(3))(0)(-)(8)NH(4)(+) clusters show that the structures of the more stable conformers follow a clear pattern: each additional NH(3) group makes a new hydrogen bond with one of the hydrogen atoms of an NH(3) unit already bound to the NH(4)(+) core. For the (NH(3))(0)(-)(8)NH(2)(-) clusters, the DFT/B3LYP calculations show that, within the calculation error, the more stable conformers follow a clear pattern for n = 1-6: each additional NH(3) group makes a new hydrogen bond to the NH(2)(-) core. For n = 7 and 8, the additional NH(3) groups bind to other NH(3) groups, probably because of the saturation of the NH(2)(-) core. Similar results were obtained at the MP2 level of calculation. A stability analysis was performed using the commonly defined stability function E(n)(-)(1) + E(n)(+1) - 2E(n), where E is the total energy of the cluster, including the zero point correction energy (E = E(t) + ZPE). The trend on the relative stability of the clusters presents an excellent agreement with the distribution of experimental cluster abundances. Moreover, the stability analysis predicts that the (NH(3))(4)NH(4)(+) and the even negative clusters [(NH(3))(n)NH(2)(-), n = 2, 4, and 6] should be the most stable ones, in perfect agreement with the experimental results.     

Effect of Ni and Pd on the geometry, electronic properties, and active sites of copper clusters

The geometry, electronic properties, and active sites of copper clusters doped with Ni or Pd atoms, Cu(n)()(-)(1)M (n = 2-6; M = Ni, Pd) have been investigated using first-principles methods. Planar structures are energetically favorable in Cu(n)()(-)(1)Ni (n = 2-6). However, for Pd-doped clusters, three-dimensional structures are competitive in energy, and for n = 6, the most stable structure is not planar. Several properties of doped copper clusters present odd-even oscillations as the number of copper atoms grow. The different atomic ground-state configuration of Ni and Pd determines the bonding and electronic properties of doped copper clusters. The interaction between impurities and copper atoms can modify the chemical hardness and active sites of doped copper clusters markedly inducing directionality in the reactivity. This effect is relevant to the behavior of catalysts as well as in the growth of metallic films.     

Geometries, stabilities, and electronic properties of different-sized ZrSi(n) (n=1-16) clusters: a density-functional investigation

The ZrSi(n) (n=1-16) clusters with different spin configurations have been systematically investigated by using the density-functional approach. The total energies, equilibrium geometries, growth-pattern mechanisms, natural population analysis, etc., are discussed. The equilibrium structures of different-sized ZrSi(n) clusters can be determined by two evolution patterns. Theoretical results indicate that the most stable ZrSi(n) (n=1-7) geometries, except ZrSi3, keep the analogous frameworks as the lowest-energy or the second lowest-energy Si(n+1) clusters. However, for large ZrSi(n) (n=8-16) clusters, Zr atom obviously disturbs the framework of silicon clusters, and the localized position of the transition-metal (TM) Zr atom gradually varies from the surface insertion site to the concave site of the open silicon cage and to the encapsulated site of the sealed silicon cage. It should be mentioned that the lowest-energy sandwich-like ZrSi12 geometry is not a sealed structure and appears irregular as compared with other TM@Si12 (TM = Re,Ni). The growth patterns of ZrSi(n) (n=1-16) clusters are concerned showing the Zr-encapsulated structures as the favorable geometries. In addition, the calculated fragmentation energies of the ZrSi(n) (n=1-16) clusters manifest that the magic numbers of stabilities are 6, 8, 10, 14, and 16, and that the fullerene-like ZrSi16 is the most stable structure, which is in good agreement with the calculated atomic binding energies of ZrSi(n) (n=8-16) and with available experimental and theoretical results. Natural population analysis shows that the natural charge population of Zr atom in the most stable ZrSi(n) (n=1-16) structures exactly varies from positive to negative at the critical-sized ZrSi8 cluster; furthermore, the charge distribution around the Zr atom appears clearly covalent in character for the small- or middle-sized clusters and metallic in character for the large-sized clusters. Finally, the properties of frontier orbitals and polarizabilities of ZrSi(n) are also discussed.