A semi-empirical potential for simulations of transition metal clusters: Minima and isomers of Nin (n=2–13) and their hydrides

First principles insights into the relative stability, electronic and catalytic properties of core-shell, Janus and mixed structural patterns for bimetallic Pd-X nano-alloys (X = Co, Ni, Cu, Rh, Ag, Ir, Pt, Au)

The three well-known orderings of the two constituting atomic species in a bimetallic nano-alloy - core-shell, Janus and mixed structural patterns - may be interconvertible depending on the synthesis conditions. Using first principles electronic structure calculations in the present work, we look for the microscopic origin for such structural transformation considering eight Pd-related bimetallic nano-alloys. Our analysis shows that it is the change in atom-atom covalency that is responsible for such structural transformation. Our study also reveals that the three patterns are distinctly identified in terms of total orbital hybridization. Finally, we have analyzed the trend in the relative catalytic activity for the three structures of each bimetallic nano-alloy using the d-band model. Our analysis indicates that the trend in the catalytic activity for the bimetallic Pd-X nano-alloys seems to be intermediate to those of the pristine Pd and Pt nano-clusters possessing similar structure and equal number of total atoms. Among the studied binary nano-alloys, the bimetallic Pd-Ni nano-alloy appears as the most suitable binary pair to develop a non-Pt catalyst.

[Ni8(CNtBu)12][Cl]: A nickel isocyanide nanocluster with a folded nanosheet structure

The reaction of 1.75 equiv of ^tBuNC with Ni(1,5-COD)₂, followed by crystallization from benzene/pentane, resulted in the isolation of [Ni₈(CN^tBu)₁₂][Cl] (2) in low yields. Similarly, the reaction of Ni(1,5-COD)₂ with 0.6 equiv of [Ni(CN^tBu)₄], followed by addition of 0.08 equiv of I₂, resulted in the formation of [Ni₈(CN^tBu)₁₂][I] (3), which could be isolated in 52% yield after work-up. Both 2 and 3 adopt folded nanosheet structures in the solid state, characterized by two symmetry-related planar Ni₄ arrays, six terminally bound ^tBuNC ligands, and six ^tBuNC ligands that adopt bridging coordination modes. The metrical parameters of the six bridging ^tBuNC ligands suggest that they have been reduced to their [^tBuNC]^2- form. In contrast to the nanosheet structures observed for 2 and 3, gas phase Ni₈ is predicted to feature a compact bisdisphenoid ground state structure. The strikingly different structural outcomes reveal the profound structural changes that can occur upon addition of ligands to bare metal clusters. Ultimately, the characterization of 2 and 3 will enable more accurate structural predictions of ligand-protected nanoclusters in the future.

Influence of Single Impurity Atoms on the Structure, Electronic, and Magnetic Properties of Ni5 Clusters

With a gradient-corrected density functional method, we have studied computationally the influence of single impurity atoms on the structure, electronic, and magnetic properties of Ni5 clusters. The square-pyramidal isomer of bare Ni5 with six unpaired electrons was calculated 23 kJ/mol more stable than the trigonal bipyramid in its lowest-energy electronic configuration with four unpaired electrons. In a previous study on the cluster Ni4, we had obtained only one stable isomer with an O or an H impurity, but we located six minima for ONi5 and five minima for HNi5. In the most stable structures of HNi5, the H atom bridges a Ni-Ni edge at the base or the side of the square pyramid, similarly to the coordination of an H atom at the tetrahedral cluster Ni4. The most stable ONi5 isomers exhibit a trigonal bipyramidal structure of the Ni5 moiety, with the impurity coordinated at a facet, (micro3-O)Ni5, or at an apex edge, (micro-O)Ni5. We located four stable structures for a C impurity at a Ni5 cluster. As for CNi4, the most stable structure of the corresponding Ni5 complex comprises a four-coordinated C atom, (micro4-C)Ni5, and can be considered as insertion of the impurity into a Ni-Ni bond of the bare cluster. All structures with C and five with O impurity have four unpaired electrons, while the number of unpaired electrons in the clusters HNi5 varies between 3 and 7. As a rough trend, the ionization potentials and electron affinities of the clusters with impurity atoms decrease with the coordination number of the impurity. However, the position of the impurity and the shape of the metal moiety also affect the results. Coordination of an impurity atom leads to a partial oxidation of the metal atoms.

Structure, stability, electronic and magnetic properties of Ni4 clusters containing impurity atoms

Using a gradient-corrected density functional method, we studied computationally how single impurity atoms affect the structure and the properties of a Ni4 cluster. H and O atoms coordinate at a Ni-Ni bond, inducing small changes to the structure of bare Ni4 which is essentially a tetrahedron. For a C impurity, we found three stable structures at a Ni4 cluster. In the most stable geometry, the carbon atom cleaves a Ni-Ni bond of Ni4, binding to all Ni atoms. Inclusion of the impurity atom leads to a partial oxidation of the metal atoms and, in the most stable structures, reduces the spin polarization of the cluster compared to bare Ni4. An H impurity interacts mainly with the Ni 4s orbitals, whereas the Ni 3d orbitals participate strongly in the bonding with O and C impurity atoms. For these impurity atoms, Ni 3d contributions dominate the character of the HOMO of the ligated cluster, in contrast to the HOMO of bare Ni4 where Ni 4s orbitals prevail. We also discuss a simple model which relates the effect of a H impurity on the magnetic state of metal clusters to the spin character (minority or majority) of the LUMO or HOMO of the bare metal cluster.

Homonuclear transition-metal trimers

Density-functional theory has been used to determine the ground-state geometries and electronic states for homonuclear transition-metal trimers constrained to equilateral triangle geometries. This represents the first application of consistent theoretical methods to all of the ten 3d block transition-metal trimers, from scandium to zinc. A search of the potential surfaces yields the following electronic ground states and bond lengths: Sc3(2A1',2.83 A), Ti3(7E',2.32 A), V3(2E",2.06 A), Cr3(17E',2.92 A), Mn3(16A2',2.73 A), Fe3(11E",2.24 A), Co3(6E",2.18 A), Ni3(3A2",2.23 A), Cu3(2E',2.37 A), and Zn3(1A1',2.93 A). Vibrational frequencies, several low-lying electronic states, and trends in bond lengths and atomization energies are discussed. The predicted dissociation energies DeltaE(M3-->M2+M) are 49.4 kcal mol(-1)(Sc3), 64.3 kcal mol(-1)(Ti3), 60.7 kcal mol(-1)(V3), 11.5 kcal mol(-1)(Cr3), 32.4 kcal mol(-1)(Mn3), 61.5 kcal mol(-1)(Fe3), 78.0 kcal mol(-1)(Co3), 86.1 kcal mol(-1)(Ni3), 26.8 kcal mol(-1)(Cu3), and 4.5 kcal mol(-1)(Zn3).

Broken-symmetry unrestricted hybrid density functional calculations on nickel dimer and nickel hydride

In the present work we investigate the adequacy of broken-symmetry unrestricted density functional theory for constructing the potential energy curve of nickel dimer and nickel hydride, as a model for larger bare and hydrogenated nickel cluster calculations. We use three hybrid functionals: the popular B3LYP, Becke's newest optimized functional Becke98, and the simple FSLYP functional (50% Hartree-Fock and 50% Slater exchange and LYP gradient-corrected correlation functional) with two basis sets: all-electron (AE) Wachters+f basis set and Stuttgart RSC effective core potential (ECP) and basis set. We find that, overall, the best agreement with experiment, comparable to that of the high-level CASPT2, is obtained with B3LYP/AE, closely followed by Becke98/AE and Becke98/ECP. FSLYP/AE and B3LYP/ECP give slightly worse agreement with experiment, and FSLYP/ECP is the only method among the ones we studied that gives an unacceptably large error, underestimating the dissociation energy of Ni(2) by 28%, and being in the largest disagreement with the experiment and the other theoretical predictions. We also find that for Ni(2), the spin projection for the broken-symmetry unrestricted singlet states changes the ordering of the states, but the splittings are less than 10 meV. All our calculations predict a deltadelta-hole ground state for Ni(2) and delta-hole ground state for NiH. Upon spin projection of the singlet state of Ni(2), almost all of our calculations: Becke98 and FSLYP both AE and ECP and B3LYP/AE predict (1)(d(A)(x(2)-y(2)d(B)(x(2)-y(2)) or (1)(d(A)(xy) (d)(B)(xy)) ground state, which is a mixture of (1)Sigma(g) (+) and (1)Gamma(g). B3LYP/ECP predicts a (3)(d(A)(x(2)-y(2))d(B)(xy) (mixture of (3)Sigma(g) (-) and (3)Gamma(u)) ground state virtually degenerate with the (1)(d(A)(x(2)-y(2)d(B)(x)(2)-y(2)/(1)(d(A)(xy)D(B)(xy) state. The doublet delta-hole ground state of NiH predicted by all our calculations is in agreement with the experimentally predicted (2)Delta ground state. For Ni(2), all our results are consistent with the experimentally predicted ground state of 0(g) (+) (a mixture of (1)Sigma(g) (+) and (3)Sigma(g) (-)) or 0(u) (-) (a mixture of (1)Sigma(u) (-) and (3)Sigma(u) (+)).

On the encapsulation of nickel clusters by molecular nitrogen

The structures and energetic effects of molecular nitrogen adsorbates on nickel clusters are investigated using an extended Huckel model coupled with two models of the adsorbate-nickel interaction. The potential parameters for the adsorbates are chosen to mimic experimental information about the binding strength of nitrogen on both cluster and bulk surface phases of nickel. The first model potential is a simple Lennard-Jones interaction that leads to binding sites in holes defined by sets of near-neighbor nickel atoms. The second model potential has a simple three-body form that forces the model nitrogen adsorbates to bind directly to single nickel atoms. Significant rearrangement of the core nickel structures are found in both models. A disconnectivity graph analysis of the potential energy surfaces implies that the rearrangements arise from low transition state barriers and the small differences between available isomers in the nickel core.

REACTIONS OFTRANSITIONMETALCLUSTERS WITHSMALLMOLECULES

Atoms and small molecules react with transition metal clusters in ways that are analogous to the physisorption and chemisorption reactions observed on the corresponding extended metal surface. However, often underlying these similarities are size-dependent variations in the reaction mechanisms and rates, the interpretation of which requires a detailed understanding of the structures of both the bare metal cluster substrates and the cluster-molecule complexes. Although polyatomic transition metal clusters cannot be characterized by the traditional methods of molecular spectroscopy, the combination of other physical and chemical probes can provide qualitative and semiquantitative structural information. These techniques, when combined with equilibrium geometries calculated using ab initio or semiempirical methods, provide a detailed picture of the structural origin of metal cluster reactivity and its variation with size.