Icosahedral gold cage clusters: M@Au12- (M=V, Nb, and Ta)

We report the observation and characterization of a series of stable bimetallic 18-valence-electron clusters containing a highly symmetric 12-atom icosahedral Au cage with an encapsulated central heteroatom of Group VB transition metals, M@Au(12) (-) (M=V,Nb,Ta). Electronic and structural properties of these clusters were probed by anion photoelectron spectroscopy and theoretical calculations. Characteristics of the M@Au(12) (-) species include their remarkably high binding energies and relatively simple spectral features, which reflect their high symmetry and stability. The adiabatic electronic binding energies of M@Au(12) (-) were measured to be 3.70+/-0.03, 3.77+/-0.03, and 3.76+/-0.03 eV for M=V, Nb, and Ta, respectively. Comparison of density-functional calculations with experimental data established the highly symmetric icosahedral structures for the 18-electron cluster anions, which may be promising building blocks for cluster-assembled nanomaterials in the form of stoichiometric [M@Au(12) (-)]X(+) salts.     

Evolution of size and shape in the colloidal crystallization of gold nanoparticles

The addition of dodecanethiol to a solution of oleylamine-stabilized gold nanoparticles in chloroform leads to aggregation of nanoparticles and formation of colloidal crystals. Based on results from dynamic light scattering and scanning electron microscopy we identify three different growth mechanisms: direct nanoparticle aggregation, cluster aggregation, and heterogeneous aggregation. These mechanisms produce amorphous, single-crystalline, polycrystalline, and core-shell type clusters. In the latter, gold nanoparticles encapsulate an impurity nucleus. All crystalline structures exhibit fcc or icosahedral packing and are terminated by (100) and (111) planes, which leads to truncated tetrahedral, octahedral, and icosahedral shapes. Importantly, most clusters in this system grow by aggregation of 60-80 nm structurally nonrigid clusters that form in the first 60 s of the experiment. The aggregation mechanism is discussed in terms of classical and other nucleation theories.     

Assessment of relative stabilities of positional isomers of polyhedral heteronuclear clusters via a simplified method of bond energy calculations based on tight-binding approach and adjacent matrix method: applications to binary icosahedral clusters

A new and simple method for assessing the relative stabilities of various positional isomers of a given heteronuclear cluster is described. The method is based on a tight-binding approach in conjunction with an adjacent matrix methodology (TBAM). The usefulness of the method is illustrated by bond energy calculations of a number of binary icosahedral clusters, including noncentered icosahedral A(n)B(12)(-)n clusters comprising main-group elements B, C, N, and S as well as B- and A-centered icosahedral A(n)B(13)(-)n clusters that consist of transition metals, Au, Ag, Ni, and Pt atoms. The latter results are compared with the previously reported molecular mechanics calculations based on Lennard-Jones potential and with experimental results, whenever possible. The trends of the total bond energies obtained by the two methods are nearly parallel in all cases, indicating that the relative stabilities predicted by the two methods follow the same order. The TBAM approach provides a simple and efficient way of predicting the relative stabilities of various positional isomers of a given cluster, particularly for clusters where the number of positional isomers is so large that it cannot be handled manually. The total bond energies exhibit a stepwise progression. Each step is characterized by a set of A-A, B-B, and A-B bonds which uniquely determines the total bond energy and, hence, the stability. The step formation implies that positional isomers of a given cluster geometry can be categorized by sets of numbers of A-A, B-B, and A-B bonds, or simply the numbers of the minority (either A-A or B-B) bonds. Three site preference rules, the strong-bond rule, the heterobond rule, and the big-hole rule, were formulated based on these model calculations. These rules are useful in rationalizing and/or predicting the relative stabilities of various positional isomers of a given cluster geometry.     

Theoretical study of the stability of lithium atoms in alpha-rhombohedral boron

The stability of lithium atoms in alpha-rhombohedral boron was investigated by first-principles calculations of total energies and molecular dynamics (MD) simulations. In the case of a low concentration (1.03 at. %), Li at the center of the icosahedral B12 site (the I-site) had a negative binding energy, which suggests Li at the I-site is unstable. However, MD simulations at temperatures below 750 K indicated that Li is still confined in the B12 cage under these conditions, which means Li at the I-site is metastable. Over 800 K, Li began to move away from the B12 site and settled at the tetrahedral site (the T-site) or at the octahedral site (the O-site). Li at the T-site also had a negative binding energy, but MD simulations indicated it was metastable up to 1400 K and did not move to other sites. Li at the O-site was energetically the most favorable, having a positive binding energy. In the case of a high concentration (7.69 at. %), the I-site changed to an unstable saddle point. At this concentration, the T-site was metastable and the O-site became the most stable. In MD simulations at 1400 K, Li atoms at the O-site never jumped to other sites regardless of concentration. Considering these facts, the diffusion coefficient of Li in alpha-rhombohedral boron would have to be very small below 1400 K.     

First principles study of the stability and electronic structure of the icosahedral La13, La(-1) (13), and La(+1) (13) clusters

The structural stability and electronic-structure of icosahedral La(13), La(-1) (13), and La(+1) (13) clusters have been studied by DMOL cluster method based on density-functional theory. The ground state of all-electron with relativity results is shown to be a distorted D(2h) icosahedron by the Jahn-Teller effect. However, the binding energies of D(3d) and D(5d) are very close to that of the D(2h) structure for La(13), La(-1) (13), and La(+1) (13) clusters. The effective core potential results show that the true ground state is D(5d) structure. The clusters have small magnetic moments and the symmetry of cluster is an important factor in determining the magnetic moments of the clusters. The effects of interatomic spacing and coordination on atomic magnetic moment are discussed. Further, 5d electrons dominate the hybrid orbitals below the Fermi level in the neutral cluster and contribute the main spin of clusters.     

Structural transition from icosahedra to decahedra of large Lennard-Jones clusters

The lowest icosahedral and decahedral energies of LJ1001-1610 clusters are obtained using a greedy search method (GSM) based on lattice construction. By comparing the lowest energies of icosahedral and decahedral clusters with the same atoms, the structural transition of LJ clusters is studied. Results show that the critical size from icosahedra to decahedra is located at N = 1034. When the cluster size is larger than 1034, the optimal structures are decahedra except the LJ1367-1422 clusters near the magic number, 1402, of icosahedra. However, the energies of icosahedra near the next magic number, 2044, are higher than that of decahedra, which implies that decahedra will be the optimal structure when the cluster size is larger than 1422, even for those clusters near the magic numbers of icosahedra.     

Chemical identification of icosahedral structure for cobalt and nickel clusters

Reactions of bare and hydrogenated cobalt and nickel clusters with ammonia and with water are used to determine cluster geometrical structure. Saturation measurements determine the total number of ammonia binding sites on cluster surfaces. A pattern of minima in the number of such sites is found to correlate with the sequence of closed shells and subshells expected for icosahedral packing in the 50- to 120-atom size range (50- to 200-atom range for hydrogenated clusters). In many cases there are 12 sites at the minima, the number that would be expected for preferred ammonia binding sites on closed (sub)shells of icosahedral clusters. The equilibrium adsorption of a single water molecule provides a sensitive measure of changes in cluster-water binding energy. A pattern of binding energy maxima is found, once again correlating with icosahedral structure, but for clusters having one metal atom more than the closed (sub)shells. In general, hydrogenation enhances the patterns of minima and maxima. These observations are explained in terms of the expected nature of ammonia and water binding to icosahedral clusters.     

Diffusion processes and growth on aluminum cluster surfaces

Diffusion processes of adatoms on icosahedral and Wulff polyhedral aluminum cluster surfaces have been studied by molecular dynamics simulations using the effective medium theory. Activation energies of diffusion mechanisms along {111} and {100} facets and from one facet to another, including different hopping and exchange processes as well as more exotic events, have been calculated. Exchange diffusion of an adatom by a chain mechanism through a {100} facet between two {111} facets and hopping diffusion across the edge between two {111} facets via a pull of another adatom on the neighbour facet are shown to play an important role. Adatoms on {111} facets are mobile already at very low temperatures, but on {100} facets diffusion starts above the room temperature as well as diffusion from {111} facets to {100} facets. Diffusion from {100} facet to other facets was not observed until at temperatures close to the melting temperatures of clusters. Dynamical simulations at different temperatures confirmed the appearance of diffusion mechanisms predicted by the activation energies.     

Lowest-energy structures of (C60)nX (X=Li+,Na+,K+,Cl-) and (C60)nYCl (Y=Li,Na,K) clusters for n</=13

Basin-hopping global optimization is used to find likely candidates for the lowest minima on the potential energy surface of (C(60))(n)X (X=Li(+),Na(+),K(+),Cl(-)) and (C(60))(n)YCl (Y=Li,Na,K) clusters with n相似文献   

Size and structure dependence of carbon monoxide chemisorption on cobalt clusters

We have carried out a series of ab initio calculations to investigate changes in the structural and magnetic properties of pristine cobalt clusters upon CO chemisorption. Our results show that binding energies of CO to 13-55 atom (0.5-1.5 nm) cobalt nanoparticles and preferred chemisorption sites depend on the cluster structure (whether fcc or icosahedral), size, and surface coverage. In addition, we find a strong influence of CO on the magnetism of the cluster, leading to magnetic moments smaller than in the bulk, at variance with pristine clusters which have magnetic moments larger than the bulk. Our findings suggest important changes in catalytic properties of cobalt at the nanoscale. Our theory suggests that at the nanoscale cluster size and surface coverage might control catalysis.     

Melting and freezing characteristics and structural properties of supported and unsupported gold nanoclusters

Molecular dynamics simulations in conjunction with MEAM potential models have been used to study the melting and freezing behavior and structural properties of both supported and unsupported Au nanoclusters within a size range of 2 to 5 nm. In contrast to results from previous simulations regarding the melting of free Au nanoclusters, we observed a structural transformation from the initial FCC configuration to an icosahedral structure at elevated temperatures followed by a transition to a quasimolten state in the vicinity of the melting point. During the freezing of Au liquid clusters, the quasimolten state reappeared in the vicinity of the freezing point, playing the role of a transitional region between the liquid and solid phases. In essence, the melting and freezing processes involved the same structural changes which may suggest that the formation of icosahedral structures at high temperatures is intrinsic to the thermodynamics of the clusters, rather than reflecting a kinetic phenomenon. When Au nanoclusters were deposited on a silica surface, they transformed into icosahedral structures at high temperatures, slightly deformed due to stress arising from the Au-silica interface. Unlike free Au nanoclusters, an icosahedral solid-liquid coexistence state was found in the vicinity of the melting point, where the cluster consisted of coexisting solid and liquid fractions but retained an icosahedral shape at all times. These results demonstrated that the structural stability in the structures of small Au nanoclusters can be enhanced through interaction with the substrate. Supported Au nanoclusters demonstrated a structural transformation from decahedral to icosahedral motifs during Au island growth, in contrast to the predictions of the minimum-energy growth sequence: icosahedral structures appear first at very small cluster sizes, followed by decahedral structures, and finally FCC structures recovered at very large cluster sizes. The simulations also showed that island shapes are strongly influenced by the substrate, more specifically, the structural characteristic of a Au island is not only a function of size, but also depends on the contact area with the surface, which is controlled by the wetting of the cluster to the substrate.     

A density-functional study of the structures, binding energies and total spins of Ni-Fe clusters using nonlocal norm-conserving pseudopotentials and the generalized gradient approximation

We report ab initio calculations of the structures, binding energies, and total spins of the clusters Ni(13), Ni(19), Ni(23), Ni(26), Ni(12)Fe, Ni(11)Fe(2), Ni(18)Fe, Ni(17)Fe(2), Ni(22)Fe, Ni(20)Fe(3), and Ni(25)Fe using a density-functional method that employs linear combination of atomic orbitals as basis sets, nonlocal norm-conserving pseudopotentials, and the generalized gradient approximation for exchange and correlation. Our results show that the Fe-doped Ni clusters, which have icosahedral or polyicosahedral ground-state structures similar to those of the corresponding pure Ni clusters, are most stable with the Fe atoms occupying internal positions, as has also been inferred from experimental results on the adsorption of molecular nitrogen on the cluster surfaces. We also rule out the possibility that the experimentally observed difference between the (nonpolyicosahedral) configurations of N(2)-saturated Ni(26) and N(2)-saturated Ni(25)Fe be due to the influence of the Fe atom on the energy of the underlying metal cluster.     

Thermal effects on energetics and dynamics in water cluster anions (H2O)n(-)

The electron binding energies and relaxation dynamics of water cluster anions (H(2)O)(n)(-) (11 ≤ n ≤ 80) formed in co-expansions with neon were investigated using one-photon and time-resolved photoelectron imaging. Unlike previous experiments with argon, water cluster anions exhibit only one isomer class, the tightly bound isomer I with approximately the same binding energy as clusters formed in argon. This result, along with a decrease in the internal conversion lifetime of excited (H(2)O)(n)(-) (25 ≤ n ≤ 40), indicates that clusters are vibrationally warmer when formed in neon. Over the ranges studied, the vertical detachment energies and lifetimes appear to converge to previously reported values.     

All-boron analogues of aromatic hydrocarbons: B17- and B18-

We have investigated the structural and electronic properties of the B(17)(-) and B(18)(-) clusters using photoelectron spectroscopy (PES) and ab initio calculations. The adiabatic electron detachment energies of B(17)(-) and B(18)(-) are measured to be 4.23 ± 0.02 and 3.53 ± 0.05 eV, respectively. Calculated electron detachment energies are compared with experimental data, confirming the presence of one planar C(2v) ((1)A(1)) isomer for B(17)(-) and two nearly isoenergetic quasi-planar C(3v) ((2)A(1)) and C(s) ((2)A') isomers for B(18)(-). The stability and planarity/quasi-planarity of B(17)(-) and B(18)(-) are ascribed to σ- and π-aromaticity. Chemical bonding analyses reveal that the nature of π-bonding in B(17)(-) and B(18)(-) is similar to that in the recently elucidated B(16)(2-) and B(19)(-) clusters, respectively. The planar B(17)(-) cluster can be considered as an all-boron analogue of naphthalene, whereas the π-bonding in the quasi-planar B(18)(-) is reminiscent of that in coronene.     

Density functional theory study of BnC clusters

B_nC clusters (n = 3–10) were studied at the density functional theory (DFT) (B3LYP)/6‐311G** level of theory. The calculations predicted that the most stable configurations of the B_nC clusters are the (n + 1)‐membered cyclic structures. For boron–carbon clusters, the configurations containing greater numbers of three‐membered boron rings are more favorable, except for the B₇C and B₉C clusters. Through molecular orbital analysis of these B_nC clusters, we have concluded that π‐electron delocalization plays a crucial role in the stability of n + 1‐membered cyclic structures. In this paper, the relative stability of each cluster is discussed based on their single atomic‐binding energies. The capability of clusters to obtain or lose an electron was also discussed, based on their vertical electron detachment energies (VDEs), adiabatic electron detachment energies (ADEs), vertical electron affinities (VEAs) and adiabatic electron affinities (AEAs). Copyright © 2011 John Wiley & Sons, Ltd.     

The structures of small clusters of C60 molecules

Using a pairwise additive atom-atom intermolecular potential to describe the interaction between C₆₀ molecules, we calculated the lowest-energy structures of (C60)N clusters up to N = 15 and compared the results with predictions derived using Girifalco’s spherical potential. The cluster binding energies calculated on the basis of the former potential are in all cases significantly higher than those obtained from the latter. Moreover, the atom-atom potential predicts that small fullerene clusters have structures based on icosahedral packing, a finding which, for N = 14 and 15, contrasts with the results obtained using Girifalco’s approximation.     

Ground state, growth, and electronic properties of small lanthanum clusters

The DMol cluster method based on density-functional theory has been employed to study the structural stability and electronic structure of La(n) (n=2-14) clusters. The ground states have been found out for lanthanum clusters. The Jahn-Teller effect plays an important role in this process because there are many isomers near the ground state. The magnetism is not sensitive to interatomic spacing when the change of interatomic spacing is in a small range. Lanthanum clusters grow in an icosahedral pattern. The results of the mean binding energy, of the second derivative of binding energy, and of the formation energy show strong odd-even alternation and that 7- and 13-atom clusters are magic. Further, the HOMO-LUMO gap, the mean nearest bond lengths, and the mean magnetic moments suggest that the convergence to bulk is slow and it shows an oscillatory behavior for small lanthanum clusters.     

Structural and electronic properties of Si n, Si n+, and AlSi n-1 (n=2-13) clusters: theoretical investigation based on ab initio molecular orbital theory

The geometric and electronic structures of Si(n), Si(n) (+), and AlSi(n-1) clusters (2< or =n< or =13) have been investigated using the ab initio molecular orbital theory under the density functional theory formalism. The hybrid exchange-correlation energy function (B3LYP) and a standard split-valence basis set with polarization functions [6-31G(d)] were employed for this purpose. Relative stabilities of these clusters have been analyzed based on their binding energies, second difference in energy (Delta (2)E) and fragmentation behavior. The equilibrium geometry of the neutral and charged Si(n) clusters show similar structural growth. However, significant differences have been observed in the electronic structure leading to their different stability pattern. While for neutral clusters, the Si(10) is magic, the extra stability of the Si(11) (+) cluster over the Si(10) (+) and Si(12) (+) bears evidence for the magic behavior of the Si(11) (+) cluster, which is in excellent agreement with the recent experimental observations. Similarly for AlSi(n-1) clusters, which is isoelectronic with Si(n) (+) clusters show extra stability of the AlSi(10) cluster suggesting the influence of the electronic structures for different stabilities between neutral and charged clusters. The ground state geometries of the AlSi(n-1) clusters show that the impurity Al atom prefers to substitute for the Si atom, that has the highest coordination number in the host Si(n) cluster. The fragmentation behavior of all these clusters show that while small clusters prefers to evaporate monomer, the larger ones dissociate into two stable clusters of smaller size.     

用密度泛函理论研究ZrnB(n=1-13)团簇的结构及性质

利用密度泛函理论, 得到了ZrnB(n=1-13)团簇的基态结构, 计算并讨论了团簇能量的二阶差分和离解能. 结果表明, n=2, 5, 12时, 相应团簇较稳定, 特别是Zr5B团簇的稳定性最高. 同时分析了ZrnB团簇的电子性质及磁性, 结果显示能隙随n值的增大出现奇偶振荡趋势, 特别是Zr12B团簇的能隙只有0.015 eV, 表明该团簇已具有金属性. 电荷转移随n值增大, 整体呈增大趋势, 除了二聚体ZrB, 电荷由B原子转移到Zr原子. 利用Mulliken布居分析得到二聚体ZrB(5.000 μB)和团簇Zr4B(3.000 μB)的磁矩较大, ZrnB团簇中总磁矩主要来自Zr原子的4d轨道.     

Structural and electronic properties of boron-doped lithium clusters: ab initio and DFT studies

The lowest-energy structures and electronic properties of the BLi(n) (n = 1-7) clusters are reported using the B3LYP, MP2, and CCSD(T) methods with the aug-cc-pVDZ basis set. Though the results at the B3LYP level agree well with those at the CCSD(T) level, the MP2 method is rather unsatisfactory. The first three-dimensional ground state in the BLi(n) clusters occurs for BLi(4), and the impurity B atom is seen to be trapped in a Li cage from the BLi(6) cluster onwards. The evolution of the binding energies, vertical ionization potentials, and polarizability with size of cluster shows the BLi(5) cluster to be most stable among the BLi(n) clusters. Besides, the BLi(5) cluster is also found to have the largest reaction enthalpy (49.8 kcal/mol) upon losing a Li atom, which is different from the previous prediction. The unique stability of the 8-valence electron BLi(5) can be understood from the cluster electronic shell model (CSM). However, in contradiction to the prediction of the CSM, the 2s level is filled prior to the 1d level in the BLi(n) clusters.