Experimental and computational study of the Zn(n)S(n) and Zn(n)S(n)+ clusters

Zinc sulfide clusters produced by direct laser ablation and analyzed in a time-of-flight mass-spectrometer, showed evidence that clusters composed of 3, 6, and 13 monomer units were ultrastable. The geometry and energies of neutral and positively charged Zn(n)S(n) clusters, up to n = 16, were obtained computationally at the B3LYP/6-311+G level of theory with the assistance of an algorithm to generate all possible structures having predefined constraints. Small neutral and positive clusters were found to have planar geometries, neutral three-dimensional clusters have the geometry of closed-cage polyhedra, and cationic three-dimensional clusters have structures with a pair of two-coordinated atoms. Physical properties of the clusters as a function of size are reported. The relative stability of the positive stoichiometric clusters provides a thermodynamic rationale for the experimental results.     

The structures and electronic states of zinc-water clusters Zn(n)(H2O)(m) (n = 1-32 and m = 1-3)

Ab initio and Density Functional Theory (DFT) calculations have been carried out for zinc-water clusters Zn(n)-(H2O)(m) (n = 1-32 and m = 1-3, where n and m are the numbers of zinc atoms and water molecules, respectively) to elucidate the structure and electronic states of the clusters and the interaction of zinc cluster with water molecules. The binding energies of H2O to zinc clusters were small at n = 2-3 (2.3-4.2 kcal mol(-1)), whereas the energy increased significantly in n = 4 (9.0 kcal mol(-1)). Also, the binding nature of H2O was changed at n = 4. The cluster size dependency of the binding energy of H2O accorded well with that of the natural population of electrons in the 4p orbital of the zinc atom. In the larger clusters (n > 20), it was found that the zinc atoms in surface regions of the zinc cluster have a positive charge, whereas those in the interior region have a negative charge with the large electron population in the 4p orbital. The interaction of H2O with the zinc clusters were discussed on the basis of the theoretical results.     

Comparison of the growth patterns of Si(n) and Ge(n) clusters (n = 25-33)

We performed an unbiased search for low-energy structures of medium-sized neutral Si n and Ge n clusters ( n = 25-33) using a genetic algorithm (GA) coupled with tight-binding interatomic potentials. Structural candidates obtained from our GA search were further optimized by first-principles calculations using density functional theory (DFT). Our approach reproduces well the lowest-energy structures of Si n and Ge n clusters of n = 25-29 compared to previous studies, showing the accuracy and reliability of our approach. In the present study, we pay more attention to determine low-lying isomers of Si n and Ge n ( n = 29-33) and study the growth patterns of these clusters. The B3LYP calculations suggest that the growth pattern of Si n ( n = 25-33) clusters undergoes a transition from prolate to cage at n = 31, while this transition appears at n = 26 from the PBE-calculated results. In the size range of 25-33, the corresponding Ge n clusters hold the prolate growth pattern. The relative stabilities and different structural motifs of Si n and Ge n ( n = 25-33) clusters were studied, and the changes of small cluster structures, when acting as building blocks of large clusters, were also discussed.     

Structures and energetics of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters

The structures and energies of Be(n)Si(n) and Be(2n)Si(n) (n = 1-4) clusters have been examined in ab initio theoretical electronic structure calculations. Cluster geometries have been established in B3LYP/6-31G(2df) calculations and accurate relative energies determined by the G3XMP2 method. The two atoms readily bond to each other and to other atoms of their own kind. The result is a great variety of low-energy clusters in a variety of structural types.     

Transition from a Bloch-Wilson to a free-electron density of states in Zn(n)- clusters

We present photoelectron spectroscopy studies on Zn(n) (-) in the size range of n=3-117. We show that zinc clusters exhibit a distinct transition in their electronic structure as a function of size. At small sizes (up to n=18) the clusters follow the Bloch-Wilson picture of the development of a metal from closed-shell atoms, exhibiting a gradual decrease of the gap between the fully occupied s band and the empty p band. For large sizes (n approximately or > 32) the band overlap allows the valence electrons to fully delocalize. This leads to an almost perfect free-electron density of states, as is demonstrated by discussing the spectra in the light of standard free-electron models and by comparison to the results obtained on sodium clusters.     

Theoretical study of the electronic states of Nb4, Nb5 clusters and their anions (Nb4-,Nb5-)

Geometries and energy separations of the various low-lying electronic states of Nb(n) and Nb(n) (-) (n=4,5) clusters with various structural arrangements have been investigated. The complete active space multiconfiguration self-consistent field method followed by multireference singles and doubles configuration interaction (MRSDCI) calculations that included up to 52x10(6) configuration spin functions have been used to compute several electronic states of these clusters. The ground states of both Nb(4) ((1)A('), pyramidal) and Nb(4) (-) ((2)B(3g), rhombus) are low-spin states at the MRSDCI level. The ground state of Nb(5) cluster is a doublet with a distorted trigonal bipyramid (DTB) structure. The anionic cluster of Nb(5) has two competitive ground states with singlet and triplet multiplicities (DTB). The low-lying electronic states of these clusters have been found to be distorted due to Jahn-Teller effect. On the basis of the energy separations of our computed electronic states of Nb(4) and Nb(5), we have assigned the observed photoelectron spectrum of Nb(n) (-) (n=4,5) clusters. We have also compared our MRSDCI results with density functional calculations. The electron affinity, ionization potential, dissociation and atomization energies of Nb(4) and Nb(5) have been calculated and the results have been found to be in excellent agreement with the experiment.     

Theoretical study of Al(n) and Al(n)O (n = 2-10) clusters

The stable structures, energies, and electronic properties of neutral, cationic, and anionic clusters of Al(n) (n = 2-10) are studied systematically at the B3LYP/6-311G(2d) level. We find that our optimized structures of Al5(+), Al9(+), Al9(-), Al10, Al10(+), and Al10(-) clusters are more stable than the corresponding ones proposed in previous literature reports. For the studied neutral aluminum clusters, our results show that the stability has an odd/even alternation phenomenon. We also find that the Al3, Al7, Al7(+), and Al7(-) structures are more stable than their neighbors according to their binding energies. For Al7(+) with a special stability, the nucleus-independent chemical shifts and resonance energies are calculated to evaluate its aromaticity. In addition, we present results on hardness, ionization potential, and electron detachment energy. On the basis of the stable structures of the neutral Al(n) (n = 2-10) clusters, the Al(n)O (n = 2-10) clusters are further investigated at the B3LYP/6-311G(2d), and the lowest-energy structures are searched. The structures show that oxygen tends to either be absorbed at the surface of the aluminum clusters or be inserted between Al atoms to form an Al(n-1)OAl motif, of which the Al(n-1) part retains the stable structure of pure aluminum clusters.     

Size- and shape-dependent polarizabilities of sandwich and rice-ball Co(n)Bz(m) clusters from density functional theory

The dipole polarizabilities of Co(n)Bz(m), (n, m = 1-4, m = n, n + 1) clusters are studied by means of an all-electron gradient-corrected density functional theory and finite field method. The dipole moments are relatively large for most of the clusters, implying their asymmetric structures. The total polarizability increases rapidly as cluster size, whereas the average polarizability shows "odd-even" oscillation with relatively large values at (n, n + 1). The polarizabilities exhibit clear shape-dependent variation, and the sandwich structures have systematically larger polarizability and anisotropy than the rice-ball isomers. The dipole polarizabilities are further analyzed in terms of the highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) gap, ionization potential, and electron delocalization volume. We conclude that the polarizability variations are determined by the interplay between the geometrical and electronic properties of the clusters.     

Density functional calculation of the structure and electronic properties of Cu(n)O(n) (n = 1-8) clusters

Ab initio simulations and calculations were used to study the structures and stabilities of copper oxide clusters, Cu(n)O(n) (n = 1-8). The lowest energy structures of neutral and charged copper oxide clusters were determined using primarily the B3LYP/LANL2DZ model chemistry. For n ≥ 4, the clusters are nonplanar. Selected electronic properties including atomization energies, ionization energies, electron affinities, and Bader charges were calculated and examined as a function of n.     

硅-硫二元团簇[(SiS2)nS]-(n=1～4)的结构和稳定性的量子化学研究

用密度泛函（DFT）方法（B3LYP／6－31＋G^*)研究了硅硫团簇[(SiS2)nS]^-(n=1-4)的可能几何构型,得到各稳定构型的电子结构,并 相应的振动频率,预测了稳定构型的振动光谱,由其稳定构型的比较可在理论上预测团簇的生长规律,并可初步预测团簇的形成机理。     

Structure and stability of (TiO2)n, (SiO2)n, and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] clusters

Structures, energetics, and vibrational spectra are investigated for small pure (TiO(2))(n), (SiO(2))(n), and mixed Ti(m)Si(n-m)O(2n) [n = 2-5, m = 1 to (n - 1)] oxide clusters by density functional theory (DFT). The BP86/ATZP level of theory is employed to obtain constitutional isomers of the oxide clusters. In accordance with previous studies, our calculations show three-dimensional compact structures are preferred for pure (TiO(2))(n) with oxo-stabilized higher hexavalent states, and linear chain structures are favored for pure (SiO(2))(n) with tetravalent states. However, the herein theoretically first reported mixed Ti(m)Si(n-m)O(2n) oxide clusters prefer either three-dimensional compact or linear chain structures depending upon the stoichiometry of the compound. Vibrational analysis of the important modes of some highly stable structures is provided. Coupled-cluster single and double excitation (with triples) [CCSD(T)] computed energy gaps for the TiO(2) dimers compare well with results from previous study. Excitation energies are computed by use of time-dependent (TD) DFT and equation-of-motion coupled-cluster calculations with singles and doubles (EOM-CCSD) for the most stable isomers.     

Lowest-energy structures and electronic properties of Na-Si binary clusters from ab initio global search

The ground state structures of neutral and anionic clusters of Na(n)Si(m) (1 ≤ n ≤ 3, 1 ≤ m ≤ 11) have been determined using genetic algorithm incorporated in first principles total energy code. The size dependence of the structural and electronic properties is discussed in detail. It is found that the lowest-energy structures of Na(n)Si(m) clusters resemble those of the pure Si clusters. Interestingly, Na atoms in neutral Na(n)Si(m) clusters are usually well separated by the Si(m) skeleton, whereas Na atoms can form Na-Na bonds in some anionic clusters. The ionization potentials, adiabatic electron affinities, and photoelectron spectra are also calculated and the results compare well with the experimental data.     

Structure and stability of binary transition-metal clusters (NbCo)n (n < or = 5): a relativistic density-functional study

Equilibrium geometries and electronic properties of binary transition-metal clusters, (NbCo)n (n < or = 5), have been investigated by means of the relativistic density-functional approach. The metal-metal bonding and stability aspects of these clusters have been analyzed on the basis of calculations. Present results show that these clusters exhibit rich structural varieties on the potential-energy surfaces. The most stable structures have a compact conformation in relatively high symmetry, in which the Nb atoms prefer to form an inner core and Co atoms are capped to the facets of the core. Such building features in clustering of the Nb/Co system are related to the order of bond strength: Nb-Nb>Nb-Co>Co-Co. As the binary cluster size increases, the Nb-Co bond may become stronger than the Nb-Nb bond in the inner niobium core, which results in a remarkable increment of the Nb-Nb bond length. Amongst these binary transition-metal clusters, the singlet (NbCo)4 in T(d) symmetry has a striking high stability due to the presence of the spherical aromaticity and electronic shell closure. The size dependence of the bond length and stability of the cluster has been explored.     

Theoretical study on (Al2O3)n (n = 1-10 and 30) fullerenes and H2 adsorption properties

The structures and stabilities of (Al2O3)n (n = 1-10 and 30) clusters were studied by means of first principles calculations. The calculated results reveal that the global minima of small (Al2O3)n (n = 1-5) clusters are cage structures with high symmetries, in which Al and O atoms are three- and two-coordinated, respectively, and are linked to neighbors via single bonds. Beyond (Al2O3)5, we calculated both cage and cage-dimer structures for (Al2O3)n (n = 6-10), and the results show that, at this size range, cage-dimer structures are more stable than cage structures. Furthermore, an onion-like motif for (Al2O3)10 was studied, and it is interesting to find that, at this size, the onion structure is more favorable than cage and cage-dimer structures. For large clusters, a shell structure of Al60O90 is suggested. Electronic properties and calculations on hydrogen adsorption of these aluminum oxide structures are reported, and we discuss their possible use as hydrogen storage materials.     

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).     

Shapeshifting: ligation by 1,4-cyclohexadiene induces a structural change in Ag5(+)

Relatively little is known about structural transformations of very small metal clusters that result from the adsorption of molecules. Here, the ligand-induced structural transformation of Ag(5)(+)(g) by 1,4-cyclohexadiene, which is capable of binding metal clusters as a bidentate ligand, is investigated using equilibrium mass spectrometry experiments and theory. Based on the measured sequential ligand binding free energies of Ag(n)(+)(cyclohexene)(m) and Ag(n)(+)(1,4-cyclohexadiene)(m) (n = 3 and 5; m up to 3), it is found that Ag(5)(+)(1,4-cyclohexadiene) is a particularly stable cluster relative to the other ion-molecule association complexes investigated. These results together with those from electronic structure calculations suggest that upon addition of 1,4-cyclohexadiene to Ag(5)(+), the metal cluster core undergoes a structural transformation from a "bowtie" structure(s), in which two Ag(2) units are bridged side-on by a central Ag atom, into a bidentate Ag(5)(+)(1,4-cyclohexadiene) structure that resembles a "razorback" arrangement of the five Ag atoms. These results raise the prospect of using multidentate ligands to transform larger ionic silver clusters from relatively compact 3D geometries into 2D elongated "razorback" nanowires. However, results from electronic structure calculations for clusters in which the razorback nanowire structural motif is propagated to larger sizes (up to Ag(9)(+)) indicate that the energy required to form such templated structures becomes increasingly unfavourable with increasing size. By calculating the vertical and adiabatic ligand binding energies, the competing effects that contribute to the energy required to form such structures, such as the metal cluster reorganization energy, can be quantified. These results indicate that the tendency for metal clusters to form compact shapes dominates other effects that contribute to the energy for forming templated nanowire structures, and this effect dramatically increases with increasing cluster size.     

Theoretical investigation of the electronic properties, and first and second harmonic generation for cadmium chalcogenide

We present the results of the ab initio theoretical study of the electronic properties, and first and second harmonic generation for CdX compounds with zinc-blende structure performed using the full potential linearized augmented plane wave method. Our calculations show that these compounds have similar structures. The valence band maximum and the conduction band minimum are located at Gamma, resulting in a direct energy gap. The energy gap of these compounds decreases when S is replaced by Se and Se by Te, in agreement with the experimental data and previous theoretical work. This can be attributed to the increase in the bandwidth of the conduction bands. The optical spectra are analyzed and the origin of some of the peaks in the spectra is discussed in terms of the calculated electronic structure. Our calculations for the linear optical properties show excellent agreement with the available experimental data.     

Ab initio calculations for the photoelectron spectra of vanadium clusters

We report ab initio calculations for the electronic and structural properties of V(n), V(n) (-), and V(n) (+) clusters up to n=8. We performed the calculations using a real-space pseudopotential method based on the local spin density approximation for exchange and correlation. This method assumes no explicit basis. Wave functions are evaluated on a uniform grid; only one parameter, the grid spacing, is used to control convergence of the electronic properties. Charged states are easily handled in real space, in contrast to methods based on supercells where Coulombic divergences require special handling. For each size and charge state, we find the lowest energy structure. Our results for the photoelectron spectra, using the optimized structure, agree well with those obtained by experiment. We also obtain satisfactory agreement with the measured ionization potential and electron affinity, and compare our results to calculations using an explicit basis.     

Small, nonstoichiometric zinc sulfide clusters

Recent experiments indicated that the formation of small, nonstoichiometric clusters Zn(n)S(m) and Zn(n)S(m)+ was possible. In this work, the ground states of these clusters, where 1 < or = n, m < or = 4, were studied using density functional theory. Global minima were found to be primarily cyclic structures in which the S-Zn-S preference for large bond angles was preserved. Ionization was shown to lead to structural relaxation and occasionally major changes in conformation. Cohesive energies are reported as a function of cluster composition. Qualitative comparisons were extracted from the energetics resulting from structural optimizations, and such comparisons appear to be consistent with the experiment. The computational data for the ZnS(n) and Zn(n)S(m) (where m > n) clusters indicated that sulfur-sulfur bonding in larger ZnS clusters could be feasible without significant energetic cost and that such structures should at least be considered.     

Theoretical studies on photoelectron and IR spectral properties of Br2.-(H2O)n clusters

We report vertical detachment energy (VDE) and IR spectra of Br2.-.(H2O)n clusters (n=1-8) based on first principles electronic structure calculations. Cluster structures and IR spectra are calculated at Becke's half-and-half hybrid exchange-correlation functional (BHHLYP) with a triple split valence basis function, 6-311++G(d,p). VDE for the hydrated clusters is calculated based on second order Moller-Plesset perturbation (MP2) theory with the same set of basis function. On full geometry optimization, it is observed that conformers having interwater hydrogen bonding among solvent water molecules are more stable than the structures having double or single hydrogen bonded structures between the anionic solute, Br2.-, and solvent water molecules. Moreover, a conformer having cyclic interwater hydrogen bonded network is predicted to be more stable for each size hydrated cluster. It is also noticed that up to four solvent H2O units can reside around the solute in a cyclic interwater hydrogen bonded network. The excess electron in these hydrated clusters is localized over the solute atoms. Weighted average VDE is calculated for each size (n) cluster based on statistical population of the conformers at 150 K. A linear relationship is obtained for VDE versus (n+3)(-1/3) and bulk VDE of Br2.- aqueous solution is calculated as 10.01 eV at MP2 level of theory. BHHLYP density functional is seen to make a systematic overestimation in VDE values by approximately 0.5 eV compared to MP2 data in all the hydrated clusters. It is observed that hydration increases VDE of bromine dimer anion system by approximately 6.4 eV. Calculated IR spectra show that the formation of Br2.--water clusters induces large shifts from the normal O-H stretching bands of isolated water keeping bending modes rather insensitive. Hydrated clusters, Br2.-.(H2O)n, show characteristic sharp features of O-H stretching bands of water in the small size clusters.