Structural and dynamical properties of transition metal clusters

Results of molecular dynamics simulation studies of structural and dynamical properties of 12-, 13-, and 14-atom transition metal clusters are presented. The calculations are carried out using a Gupta-like potential expressed in reduced units. The transformation to absolute units involves two size-dependent parameters which effectively convert the potential into a size-dependent one. The minimum energy geometries of the clusters are obtained through the technique of simulated thermal quenching. A melting-like transition is observed as the energy of the clusters is increased. A novel element of the transition is that it may involve a premelting state.     

Structural and energetic properties of Ni-Cu bimetallic clusters

The lowest-energy structures for all compositions of Ni n Cu m bimetallic clusters with N = n + m up to 20 atoms, N = 23, and N = 38 atoms have been determined using a genetic algorithm for unbiased structure optimization in combination with an embedded-atom method for the calculation of the total energy for a given structure. Comparing bimetallic clusters with homoatomic clusters of the same size, it is shown that the most stable structures for each cluster size are composed entirely of Ni atoms. Among the bimetallic clusters in the size range N = 2-20, the Ni N-1 Cu 1 clusters possess the highest stability. Further, it has been established that most of the bimetallic cluster structures have geometries similar to those of pure Ni clusters. The size N = 38 presents a special case, as the bimetallic clusters undergo a dramatic structural change with increasing atom fraction of Cu. Moreover, we have identified an icosahedron, a double, and a triple icosahedron with one, two, and three Ni atoms at the centers, respectively, as particularly stable structures. We show that in all global-minimum structures Ni atoms tend to occupy mainly high-coordination inner sites, and we confirm the segregation of Cu on the surface of Ni-Cu bimetallic clusters predicted in previous studies. Finally, it is observed that, in contrast to the bulk, the ground-state structures of the 15-, 16-, and 17-atom bimetallic clusters do not experience a smooth transition between the structures of the pure copper and the pure nickel clusters as a function of the relative number of the two types of atoms. For these sizes, the concentration effect on energy is more important than the geometric one.     

Structural and electronic properties of compound metal clusters

The equilibrium geometries, relative stabilities, and vertical ionization potentials of compound clusters involving Li _n , Na, Mg, and Al atoms have been calculated using ab initio self-consistent field linear combination of atomic orbitals — molecular orbital (SCF-LCAO-MO) method. The exchange energies are calculated exactly using the unrestricted Hartree-Fock (UHF) method whereas the correlation correction is included within the framework of configuration interaction involving pair excitations of valence electrons. While the later correction has no significant effect on the equilibrium geometries of clusters, it is essential for the understanding of relative stabilities. Clusters with even numbers of electrons are found to be more stable than those with odd numbers of electrons regardless of their charge state and atomic composition. The equilibrium geometries of homo-nuclear clusters can be significantly altered by replacing one of its constituent atoms with a hetero-nuclear atom. The role of electronic structure on the geometries and stabilities of compound clusters is discussed.     

Structural properties of 3-dimensional carbon clusters

Structural properties of 3d carbon clusters were calculated employing recently developed model potential energy functions for carbon. Primarily, spherical shell structures were included in the present investigation. Configurations corresponding to local energy minima were calculated for various shells of an icosahedron containing different number of C atoms. For C₆₀, the two low-lying isomers, the buckminsterfullerene and truncated dodecahedron, were found to be almost isoenergetic. It was also found that fully relaxed structures of C₉₀ and C₁₂₀ have energies very comparable to that of C₆₀. Furthermore, a systematic analysis carried out in this study for carbon clusters with varying dimensionalities, revealed an interesting relationship between the bond lengths and the distribution of bond angles. In all cases, shorter bond distances were found to be associated with larger bond angles.     

Structural and electronic properties of hydrated MgO nanotube clusters

Hydrated MgO nanotube clusters are constructed and studied by the density functional theory at the B3LYP/6-31G(d) level. A strong exothermicity chemisorption reactivity of MgO nanotube clusters with water, which releases 137.5–171.8 kJ/mol. The averaged charge of Mg ions is steady, and presents a stronger ionic bonding. Mg ions are more sensitive to the coordination number. For the reaction of water onto clusters, electronic properties of hydrated clusters have remarkable change compared with anhydrous clusters.     

Structures and energetics of hydrated oxygen anion clusters

Hydration of the atomic oxygen radical anion is studied with computational electronic structure methods, considering (O(-))(H(2)O)(n) clusters and related proton-transferred (OH(-))(OH)(H(2)O)(n)(-)(1) clusters having n = 1-5. A total of 67 distinct local-minimum structures having various interesting hydrogen bonding motifs are obtained and analyzed. On the basis of the most stable form of each type, (O(-))(H(2)O)(n)) clusters are energetically favored, although for n > or = 3, there is considerable overlap in energy between other members of the (O(-))(H(2)O)(n) family and various members of the (OH(-))(OH)(H(2)O)(n)(-)(1) family. In the lower-energy (O(-))(H(2)O)(n) clusters, the hydrogen bonding arrangement about the oxygen anion center tends to be planar, leaving the oxygen anion p-like orbital containing the unpaired electron uninvolved in hydrogen bonding with any water molecule. In (OH(-))(OH)(H(2)O)(n)(-)(1) clusters, on the other hand, nonplanar arrangements are the rule about the anionic oxygen center that accepts hydrogen bonds. No instances are found of OH(-) acting as a hydrogen bond donor. Those OH bonds that form hydrogen bonds to an anionic O(-) or OH(-) center are significantly stretched from their equilibrium value in isolated water or hydroxyl. A quantitative inverse correlation is established for all hydrogen bonds between the amount of the OH bond stretch and the distance to the other oxygen involved in the hydrogen bond.     

Bonding and energetics in small clusters of gallium and arsenic

Density functional theory is used to study the electronic, geometric, and bonding properties of small Ga_n, As_m, and Ga_nAs_m clusters with up to eight atoms. New ground state structures for many of these clusters are reported. Electronic states, stabilities, and charges are also reported and discussed in terms of the bonding in the clusters. © 2003 Wiley Periodicals, Inc. Heteroatom Chem 14:189–196, 2003; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/hc.10127     

Structural, energetic, and electronic properties of hydrogenated titanium clusters

Hydrogen undergoes dissociative chemisorption on small titanium clusters. How the electronic structure of the cluster changes as a function of the number of adsorbed hydrogen atoms is an important issue in nanocatalysis and hydrogen storage. In this paper, a detailed theoretical investigation of the structural, energetic, and electronic properties of the icosahedral Ti13 cluster is presented as a function of the number of adsorbed hydrogen atoms. The results show that hydrogen loaded Ti13H20 and Ti13H30 clusters are exceptionally stable and are characterized by hydrogen multicenter bonds. In Ti13H20, the dissociated hydrogen atoms are bound to each of the 20 triangular faces of Ti13, while in Ti13H30, they are bound to the 30 Ti-Ti edges of Ti13. Consequently, the chemisorption and desorption energies of the Ti13H20 (1.93 eV, 3.10 eV) are higher than that of Ti13H30 (1.13 eV, 1.95 eV). While increased hydrogen adsorption leads to an elongation of the Ti-Ti bonds, there is a concomitant increase in the electrostatic interaction between the dissociated hydrogen atoms and the Ti13 cluster. This enhanced interaction results from the participation of the subsurface titanium atom at higher hydrogen concentrations. Illustrative results of hydrogen saturation on the larger icosahedral Ti55 cluster are also discussed. The importance of these results on hydrogen saturated titanium clusters in elucidating the mechanism of hydrogen adsorption and desorption in titanium doped complex metal hydrides is discussed.     

Structural properties and thermodynamics of water clusters: a Wang-Landau study

The temperature dependence of structural properties and thermodynamic behavior of water clusters has been studied using Wang-Landau sampling. Four potential models, simple point charge/extended (SPC/E), transferable intermolecular potential 3 point (TIP3P), transferable intermolecular potential 4 point (TIP4P), and Gaussian charge polarizable (GCP), are compared for ground states and properties at finite temperatures. Although the hydrogen bond energy and the distance of the nearest-neighbor oxygen pair are significantly different for TIP4P and GCP models, they approach to similar ground state structures and melting transition temperatures in cluster sizes we considered. Comparing with TIP3P, SPC/E model provides properties closer to that of TIP4P and GCP.     

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.     

Structural,magnetic, and vibrational properties of stoichiometric clusters of CrN

We performed density‐functional‐theoretic calculations to investigate the structural, magnetic and vibrational properties of the stoichiometric clusters (CrN)_n ( ). We show that the building block of the ground‐state structures of these clusters is a square CrNCrN unit; the only exception with n > 2 occurs for (CrN)₃, but this cluster has an isomer not very far in energy from the ground state consisting of a pair of CrNCrN squares sharing a CrN bond. In the smaller CrN, (CrN)₂, and (CrN)₃ clusters the magnetic moments of the N atoms are non‐negligible and antiparallel to those of the Cr atoms, but for the larger species (CrN)₄, (CrN)₅, (CrN)₆, and (CrN)₉ the cluster magnetic moments are almost entirely due to the Cr atoms. Lack of imaginary vibrational frequencies in the predicted ground‐state structures of (CrN)_n ( ) confirms that they are mechanically stable equilibrium states. © 2015 Wiley Periodicals, Inc.     

Electron bubbles in helium clusters. I. Structure and energetics

In this paper we present a theoretical study of the structure, energetics, potential energy surfaces, and energetic stability of excess electron bubbles in ((4)He)(N) (N=6500-10(6)) clusters. The subsystem of the helium atoms was treated by the density functional method. The density profile was specified by a void (i.e., an empty bubble) at the cluster center, a rising profile towards a constant interior value (described by a power exponential), and a decreasing profile near the cluster surface (described in terms of a Gudermannian function). The cluster surface density profile width (approximately 6 A) weakly depends on the bubble radius R(b), while the interior surface profile widths (approximately 4-8 A) increase with increasing R(b). The cluster deformation energy E(d) accompanying the bubble formation originates from the bubble surface energy, the exterior cluster surface energy change, and the energy increase due to intracluster density changes, with the latter term providing the dominant contribution for N=6500-2 x 10(5). The excess electron energy E(e) was calculated at a fixed nuclear configuration using a pseudopotential method, with an effective (nonlocal) potential, which incorporates repulsion and polarization effects. Concurrently, the energy V(0) of the quasi-free-electron within the deformed cluster was calculated. The total electron bubble energies E(t)=E(e)+E(d), which represent the energetic configurational diagrams of E(t) vs R(b) (at fixed N), provide the equilibrium bubble radii R(b) (c) and the corresponding total equilibrium energies E(t) (e), with E(t) (e)(R(e)) decreasing (increasing) with increasing N (i.e., at N=6500, R(e)=13.5 A and E(t) (e)=0.86 eV, while at N=1.8 x 10(5), R(e)=16.6 A and E(t) (e)=0.39 eV). The cluster size dependence of the energy gap (V(0)-E(t) (e)) allows for the estimate of the minimal ((4)He)(N) cluster size of N approximately 5200 for which the electron bubble is energetically stable.     

Structural and electronic properties of small beryllium clusters: a theoretical study

Geometric structures and electronic properties of small beryllium clusters (Be(n), 2< or = n< or =9) are investigated within the gradient-corrected density functional theory. The computations are performed with the Becke exchange and Perdew-Wang correlation functionals. Both low and high multiplicity states are considered. A predominance of higher multiplicity states among the low-energy isomers of the larger clusters is found. An analysis of the variations in the structural and electronic properties with cluster size is presented, and the results are compared with those of earlier studies.     

Structural, electronic, and chemical properties of multiply iodized aluminum clusters

The electronic structure, stability, and reactivity of iodized aluminum clusters, which have been investigated via reactivity studies, are examined by first-principles gradient corrected density functional calculations. The observed behavior of Al13I(x)- and Al14I(x)- clusters is shown to indicate that for x < or = 8, they consist of compact Al13- and Al14++ cores, respectively, demonstrating that they behave as halogen- or alkaline earth-like superatoms. For x > 8, the Al cores assume a cagelike structure associated with the charging of the cores. The observed mass spectra of the reacted clusters reveal that Al13I(x)- species are more stable for even x while Al14I(x)- exhibit enhanced stability for odd x(x > or = 3). It is shown that these observations are linked to the formation and filling of "active sites," demonstrating a novel chemistry of superatoms.     

Structural properties and energetics of Li2FeSiO4 polymorphs and their delithiated products from first-principles

Structural properties, thermodynamic stability and delithiation process for Li(2)FeSiO(4) polymorphs are investigated by using density functional theory (DFT) within the DFT + U framework. Three Li(2)FeSiO(4) polymorphs crystallizing in space group Pmn2(1), P2(1)/n, and Pmnb have been considered. The investigations demonstrate that the strong Si-O bonds remain almost unchanged during the lithiation-delithiation process for all the polymorphs, which contribute significantly to the structural stability. On the other hand, the differences in local environment around FeO(4) tetrahedra will be translated into varying degrees of distortion, which shows a significant influence on the structural stability and average voltages. The average voltages obtained here are in good agreement with the experimental values. Furthermore, the possibility of extracting more than one lithium ions per formula unit from Li(2)FeSiO(4) of P2(1)/n is also discussed.     

Theoretical study on the structure and energetics of alkali halide clusters

The structures of alkali halide clusters Na_nF_n, Li_nF_n and Na_nCl_n, and their metal-excess clusters Na_nF_n−1⁺, Li_nF_n−1⁺ and Na_nCl_n−1⁺ were investigated by the ab initio molecular orbital method for cluster sizes from 1 to 14. The magic numbers for the neutral clusters Na_nF_n, Li_nF_n and Na_nCl_n are 4, 6, and 8. The most stable structure for these cluster sizes is a perfect crystallite for Na_nF_n and Na_nCl_n, and a double ring for Li_nF_n. The magic numbers for the metal-excess clusters are 5 and 8, which are near ideal cuboids with (100) facets.     

Structures, energetics, and spectra of hydrated hydroxide anion clusters

The structures, energetics, electronic properties, and spectra of hydrated hydroxide anions are studied using density functional and high level ab initio calculations. The overall structures and binding energies are similar to the hydrated anion clusters, in particular, to the hydrated fluoride anion clusters except for the tetrahydrated clusters and hexahydrated clusters. In tetrahydrated system, tricoordinated structures and tetracoordinated structures are compatible, while in pentahydrated systems and hexahydrated systems, tetracoordinated structures are stable. The hexahydrated system is similar in structure to the hydrated chloride cluster. The thermodynamic quantities (enthalpies and free energies) of the clusters are in good agreement with the experimental values. The electronic properties induced by hydration are similar to hydrated chloride anions. The charge-transfer-to-solvent energies of these hydrated-hydroxide anions are discussed, and the predicted ir spectra are used to explain the experimental data in terms of the cluster structures. The low-energy barriers between the conformations along potential energy surfaces are reported.     

Metallic ions and clusters: Formation,energetics, and reaction

Several aspects of the properties of metal clusters and metals in cluster ions are discussed. In particular, results from our laboratory on the photoionization and reactions of sodium clusters, and the thermochemistry and formation of clusters containing metallic ions are presented.     

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.     

Structural fluctuation and atom-permutation in transition-metal clusters

The atomic structure and thermodynamic properties of transition-metal clusters containingN atoms are investigated forN=6 and 7 using the method of molecular dynamics, where Gupta's potential taking into account many-body interaction is employed. The caloric curve (total energy — temperature curve) and the structural fluctuations are studied. The “fluctuating state” is found forN=6 in the region of the temperature near below the melting point, where clusters undergo structural transition from one isomer to others without making any topological change. The fluctuating state differs from the coexistence state in that the former involves no atomic diffusion, and goes to a structural phase transition of the bulk whenN is increased. On the other hand, the motion of atom-permutation is found in the low-temperature region of the liquid state, being induced by the cooperative motion of two atoms. It is discussed that such a motion easily occurs along the surface and may be considered to be one of the characteristics of small clusters. The fluctuating state is discussed in relation to the structural fluctuation of gold clusters observed experimentally.