Size-dependent properties of Zn(m)S(n) clusters: a density-functional tight-binding study

We present the results of our theoretical calculations on structural and electronic properties of ligand-free Zn(n)S(n) [with n ranging from 4 to 104 (0.8-2.0-nm diameter)] clusters as a function of size of the clusters. We have optimized the structure whereby our initial structures are spherical parts of either zinc-blende or wurtzite structure. We have also considered some hollow bubblelike structures. The calculations are performed by using a parametrized linear combination of atomic orbitals-density-functional theory-local-density approximation-tight-binding method. We have focused on the variation of radial distribution function, Mulliken populations, electronic energy levels, band gap, and stability as a function of size for both zinc-blende and wurtzite-derived ZnS clusters. We have also reported the results of some nonstoichiometric Zn(m)S(n) (with m+n=47, 99, 177) clusters of zinc-blende modification.     

Structure and properties of ZnS nanoclusters

Earlier studies have shown that the most stable structures for (ZnS)n clusters with n = 10-47 are hollow polyhedral clusters ("bubbles"). We report a detailed study of larger clusters, where n = 50, 60, 70, and 80, for which onionlike or "double bubble" structures are predicted. We report calculations of the vibrational spectra and the electronic structure of bubble and double bubble clusters, which may assist in their experimental identification.     

Molecular dynamics study of the thermal entropy in mixed zinc chalcogenides

Molecular dynamics (MD) calculations were performed to determine the vibrational contribution to the entropy of mixing and its importance for the mixing of ZnO/ZnS and ZnS/Zn3P2. These systems were modeled by cyclic clusters Zn48O48, Zn32S32, and Zn48P32. The mixed cyclic clusters considered were Zn48O47S, Zn32S31O, Zn33S30P2, and Zn47S2P30. For each of the clusters, the entropy was calculated in the range of the experimental temperature of the mixing process. The convergence of the entropy with respect to the number of MD steps was studied. Finally, the thermal part of the entropy of mixing was determined, and its dependence on the number of MD steps was investigated. It was found that the thermal entropy is important for the Gibbs free energy of mixing near the miscibility gaps.     

Interaction behavior study of HCl on (ZnS)n (n = 1–12) clusters and HCl effect on Hg0 adsorption

The interaction behavior of HCl and (ZnS)_n (n = 1–12) clusters and HCl effect on Hg⁰ adsorbed by (ZnS)_n have been studied theoretically. The combined genetic algorithm and density functional theory (GA-DFT) method has been used to obtain the structures of (ZnS)_nHCl and (ZnS)_nHgHCl (n = 1–12) clusters. The structural properties of (ZnS)_nHCl and (ZnS)_nHgHCl have been analyzed. The adsorption energies and interaction energies have been calculated. Bond length and bond order analysis has revealed that S H and Zn Cl bonds form after HCl adsorbed on (ZnS)_n clusters, while Hg⁰ can only weakly bind with (ZnS)_nHCl clusters. According to thermodynamic adsorption analysis, the formation of (ZnS)_nHCl clusters from (ZnS)_n and HCl are spontaneous because of their negative Gibbs free energy changes. The formation of (ZnS)_nHgHCl from (ZnS)_nHCl and Hg are nonspontaneous for n = 1–4 and 9, and the Gibbs free energy changes have small negative values for other sizes. Electron localization function and noncovalent interaction (NCI) analysis of (ZnS)₁₀HgHCl manifest that Hg and its nearest Zn form zinc amalgam. Projected density of state study has been performed to obtain the interaction nature of HCl and (ZnS)_n clusters and Hg⁰ adsorption on (ZnS)_nHCl clusters. Based on our study, HCl is chemical adsorbed by (ZnS)_n clusters except (ZnS)₄ cluster. After (ZnS)_n adsorbs HCl, Hg⁰ can physically adsorb on (ZnS)_nHCl clusters. The strength of Hg⁰ on (ZnS)_nHCl is comparable to that of Hg⁰ on (ZnS)_n, indicating that HCl can hardly affect the adsorption of Hg⁰ on ZnS clusters.     

Chiral II-VI semiconductor nanostructure superlattices based on an amino acid ligand

Reaction of L-cysteine with M(NO3)2 x xH2O (M = Cd, Zn) generates M(L-cysteinate), which feature one-dimensional substructures that can be viewed as fragments of bulk structures of CdS (rock salt high pressure phase) and ZnS (wurtzite) because of the bridging modes accessible to the sulfur atom of L-cysteine. The MS substructures are arranged in a regular and periodic fashion within the crystal via the carboxylate function of L-cysteine. Considering the structural similarities with bulk materials, the optical properties of M(L-cysteinate) were studied and indicate blue shifts of the band gap of 2.59 eV (M = Cd, compared to CdS rock salt) and 1.37 eV (M = Zn, compared to ZnS wurtzite) with respect to the bulk MS structures, due to the low dimensionality of the metal-sulfur arrangement. The chelating nature of the cysteine ligand imposes an unusual mer arrangement of three binding S moieties at Cd with a correspondingly high Cd coordination number in a chalcogenide-based material. Density of states calculations show strong electronic structure similarities with the bulk phases and rationalize the band gap changes.     

Adsorption of molecular hydrogen and hydrogen sulfide on Au clusters

The authors present theoretical results describing the adsorption of H2 and H2S molecules on small neutral and cationic gold clusters (Au(n)((0/+1)), n=1-8) using density functional theory with the generalized gradient approximation. Lowest energy structures of the gold clusters along with their isomers are considered in the optimization process for molecular adsorption. The adsorption energies of H2S molecule on the cationic clusters are generally greater than those on the corresponding neutral clusters. These are also greater than the H2 adsorption energies on the corresponding cationic and neutral clusters. The adsorption energies for cationic clusters decrease with increasing cluster size. This fact is reflected in the elongations of the Au-S and Au-H bonds indicating weak adsorption as the cluster grows. In most cases, the geometry of the lowest energy gold cluster remains planar even after the adsorption. In addition, the adsorbed molecule gets adjusted such that its center of mass lies on the plane of the gold cluster. Study of the orbital charge density of the gold adsorbed H2S molecule reveals that conduction is possible through molecular orbitals other than the lowest unoccupied molecular orbital level. The dissociation of the cationic Au(n)SH2+ cluster into Au(n)S+ and H2 is preferred over the dissociation into Au(m)SH2+ and Au(n-m), where n=2-8 and m=1-(n-1). H2S adsorbed clusters with odd number of gold atoms are more stable than neighboring even n 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).     

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.     

Ab initio study of structural stability of Mo-S clusters and size specific stoichiometries of magic clusters

Using first-principles calculations with ultrasoft pseudopotential formalism and the generalized gradient approximation for the exchange-correlation functional, we study the stability of MonSm (n =1-6 and m ranging from n to 3n) clusters and obtain the optimal stoichiometry for each n corresponding to the magic cluster. It is found that in this size range, the lowest-energy structures favor a core of metal atoms, which is covered by sulfur. In particular, we observe that for Mo6S14 isolated clusters, a 3D structure is significantly lower in energy as compared to platelet structures found recently on Au (111) surface. The composition ratio between S and Mo in the magic clusters is less than 2 for n=3 and greater than 2 for n<3. The structural stability of the magic clusters arises from the optimization of the Mo-Mo and Mo- S bonding as well as the symmetry of the cluster. Addition of a terminal sulfur in a magic cluster generally lowers its binding energy. The presence of partially occupied d-orbitals in Mo atoms contributes to Mo-Mo bonding and for higher S concentration it leads to sulfur-sulfur bond formation. The variation in energy due to a change in the sulfur composition suggests that sulfurization of the magic clusters is generally more favorable than desulfurization.     

Structures and energetics of neutral and ionic silicon-germanium clusters: density functional theory and coupled cluster studies

We have calculated the structural and energetic properties of neutral and ionic (singly charged anionic and cationic) semiconductor binary silicon-germanium clusters Si(m)Ge(n) for s = m + n ≤ 12 using the density functional theory (DFT-B3LYP) and coupled cluster [CCSD(T)] methods with Pople's 6-311++G(3df, 3pd) basis set. Neutral and anionic clusters share similar ground state structures for s = 3-7, independent of the stoichiometry and atom locations, but start to deviate at s = 8. The relative energetic stability of the calculated ground state structures among possible isomers has been analyzed through a bond strength propensity model where the pair interactions of Si-Si, Si-Ge, and Ge-Ge are competing. Electron affinities, ionization potentials, energy gaps between the highest and lowest occupied molecular orbitals (HOMO-LUMO gaps), and cluster mixing energies were calculated and analyzed. Overall, for a fixed s, the vertical ionization potential increases as the number of silicon atoms m increases, while the vertical electron affinity shows a dip at m = 2. As s increases, the ionization potentials increase from s = 2 to s = 3 and then decrease slowly to s = 8. The mixing energies for neutral and ionic clusters are all negative, indicating that the binary clusters are more stable than pure elemental clusters. Except for s = 4 and 8, cationic clusters are more stable than anionic ones and, thus, are more likely to be observed in experiments.     

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.     

ZnS bubble clusters with onion-like structures

Following recent studies which showed that the most stable structures for (ZnS)(n) clusters (n= 10-47) are the so-called "bubble clusters", in which all the atoms are three-coordinated, we have used simulated annealing techniques to find the most stable structure for a larger cluster, (ZnS)(60). We find an onion-like structure, with one small cluster enclosed inside a bigger one. The inner cluster has the structure of a sodalite cage. Bonding between the inner and the outer clusters creates a network of four-coordinated atoms.     

Structures and stability of B-doped Al clusters: AlnB and AlnB2 (n=1-7)

Various structural possibilities for Al(n)B(m) (n=1-7, m=1-2) neutral isomers were investigated using B3LYP6-311G(d) and CCSD(T)6-311G(d) methods. Our calculations predicted the existence of a number of previously unknown isomers. The B atom favors to locate over/inside of all clusters in this series. All structures of the Al(n)B (n=2-7) may be derived from capping/putting a B atom over/inside the Al(n) cluster. All Al(n)B(2) (n=1-5) may be understood as two substitutions of Al atoms by B atoms in the Al(n+2) molecule. The strong B-B bond is a dominant factor in the building-up principle of mixed Al(n)B(2) neutral clusters. The second difference in energy showed that the Al(n)B(m) clusters with even n+m are more stable than those with odd n+m. Our results and analyses revealed that the mixed Al-B clusters exhibit aromatic behaviors.     

Structure and energetics of equiatomic K-Cs and Rb-Cs binary clusters

The basin-hopping algorithm combined with the Gupta many-body potential is used to study the structural and energetic properties of (KCs)(n) and (RbCs)(n) bimetallic clusters with N=2n up to 50 atoms. Each binary structure is compared to those of the pure clusters of the same size. For the cluster size N=28 and for the size range of N=34-50, the introduction of K and Rb atoms in the Cs alkali metal cluster results in new ground state structures different from those of the pure elements. In the size range N>/=38 the binary and pure clusters show not only structural differences, but they also display different magic numbers. Most of the magic Rb-Cs and K-Cs clusters possess highly symmetric structures. They belong to a family of pIh structures, where a fivefold pancake is a dominant structural motif. Such geometries have not been reported for alkali binary clusters so far, but have been found for series of binary transition metal clusters with large size mismatch. Moreover, tendency to phase separation (shell-like segregation) is predicted for both K-Cs and Rb-Cs clusters with up to 1000 atoms. Our finding of a surface segregation in Rb-Cs clusters is different from that of theoretical and experimental studies on bulk Rb-Cs alloys where phase separation does not occur.     

Ultrathin ZnS single crystal nanowires: controlled synthesis and room-temperature ferromagnetism properties

Highly uniform single crystal ultrathin ZnS nanowires (NWs) with 2 nm diameter and up to 10 μm length were fabricated using a catalyst-free colloidal chemistry strategy. The nanowires crystallized in hexagonal phase structure with preferential growth along the direction of the (001) basal plane. The strong polarity of the (001) plane composed of Zn cations or S anions drives the oriented attachment of ZnS nanocrystals (NCs) along this direction via electrostatic (or dipole) interaction. The ultrathin ZnS nanowires show intrinsic ferromagnetism at room temperature and other unusual properties related to its unique nature, such as large anisotropic lattice expansion, large blue-shift of UV-vis absorption band of the excition, and photoluminescence spectrum of the exciton band edge. First-principles DFT computation results show that Zn vacancies can induce intrinsic ferromagnetism in these undoped ZnS NWs. The main source of the magnetic moment arises from the unpaired 3p electrons at S sites surrounding the Zn vacancies carrying the magnetic moment ranging from 0.26 to 0.66 μ(B). Calculated results indicate that the magnetic moment of the ultrathin ZnS NWs can be increased by increasing the Zn vacancy concentration without significant energy cost. The calculated magnetization value (1.96 or 0.40 emu/g for Zn vacancies on the surface of NWs or inside, respectively) by Zn(53)S(54) supercell model is larger than our experimental value (0.12 emu/g at 1.8 K and 0.05 emu/g at 300 K), but the ferromagnetic result is qualitatively in agreement.     

The N-alkyldithiocarbamato complexes [M(S2CNHR)2] (M=Cd(II) Zn(II); R=C2H5, C4H9, C6H13, C12H25); their synthesis, thermal decomposition and use to prepare of nanoparticles and nanorods of CdS

A series of N-alkyldithiocarbamato complexes [M(S2CNHR)2] (M=Cd(II), Zn(II); R=C2H5, C4H9, C6H13, C12H25) have been synthesised and characterized. The decomposition of these complexes to sulfates has been investigated, and a mechanism proposed. The structures of [Zn(S2CNHHex)2], [Cd(SO4)2(NC5H5)4)]n and [Cd(SO4)2(NC5H5)2(H2O)2)]n have been determined by X-ray single crystal method. The cadmium complex [Cd(S2CNHC12H25)2] and zinc complex [Zn(S2CNHC6H13)2] were used as single-source precursors to synthesize CdS and ZnS nanoparticles, respectively. The synthesis of CdS nanoparticles was carried under various thermolysis conditions and changes in the shape of derived nanoparticles were studied by transmission electron microscope (TEM).     

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.     

Geometrical and electronic structures of Au(m)Ag(n) (2 < or = m + n < or = 8)

The structural and electronic properties of Au(m)Ag(n) binary clusters (2 < or = m + n < or = 8) have been investigated by density functional theory with relativistic effective core potentials. The results indicate that Au atoms tend to occupy the surface of Au(m)Ag(n) clusters (n > or = 2 and m > or = 2). As a result, segregation of small or big bimetallic clusters can be explained according to the atomic mass. The binding energies of the most stable Au(m)Ag(n) clusters increase with increasing m+n. The vertical ionization potentials of the most stable Au(m)Ag(n) clusters show odd-even oscillations with changing m+n. The possible dissociation channels of the clusters considered are also discussed.     

Transition-metal-doped aluminum hydrides as building blocks for supramolecular assemblies

Density functional theory calculations were carried out to characterize a series of transition-metal-doped aluminum hydrides, forming TMAl(n)H(2n) and TMAl(n)H(2n+1) (TM = Sc, Ti, V; n = 3,4), in either charged or neutral form. A new electron-counting rule for these clusters was formulated as PSEN (paired skeleton electron number) = 4n, which can characterize both closed-shell and open-shell clusters. On the basis of this electron-counting rule, the superatomic clusters such as TiAl(4)H(9) and TiAl(3)H(6) were identified and can be used to assemble supramolecular structures. Electronic structure analysis showed that three-centered TM-H-Al bonds largely contributed to the structural stability. Also, the spin state of a wide range of clusters in their ground state can be predicted by the electron-counting rule.     

A G3 study of the structure of carbon-nitrogen nanoclusters

Possible structures of the carbon-nitrogen clusters of the form C(m)N(n) (m = 1-4, n = 1-4, m + n = 2-5) were predicted for the neutral, anion, and cation species in the singlet, doublet, and triplet states, whenever appropriate. The calculations were performed at the G3, MP2(fc)/6-311+G*, and B3LYP/6-311+G* levels of theory. Several molecular properties related to the experimental data--such as the electronic energy, equilibrium geometry, binding energy, HOMO-LUMO gap (HLG), and spin contamination --were calculated. In addition the vertical electron attachment, the adiabatic electron affinity, and vertical ionization energy, of the neutral clusters were calculated. Most of the predicted lowest energy structures were linear, whereas bent structures became more stable with the increase of the cluster size and increase of the number of the N atoms. In most of the predicted lowest energy structures, the N atom prefers the terminal position with acetylenic bond. The calculated BE of the predicted clusters increases with the increase of the cluster size for the neutral and cation clusters but decreases with the increase of the cluster size for the anion clusters. The predicted clusters are characterized by high HLG of about 11 eV on the average, with that of the anion clusters is smaller than that for the neutral and cation clusters. It is concluded then that the anion clusters are less stable than the corresponding neutral and cation clusters. Finally, the N(2) loss reaction is treated.